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Introduction 



Radiative corrections to jet production in e~^e~ annihilation at order were com- 
puted a long time ago [jl], ^. These calculations were, however, performed for 
massless quarks. 

This is sufficient in most practical applications: in fact, if we consider the heaviest 
quark that can be produced nowadays at LEP, the b quark, we can say that, at 
relatively low energy, the b fraction is strongly suppressed, and, at high energy (i.e. 
on the Z peak and beyond), mass effects can be neglected for many observables, 
because they appear as powers of the ratio of the mass over the total energy. 

Nevertheless there are several reasons why a massive next-to-leading-order calcu- 
lation is desirable. 

- First of all, there are quantities, such as the heavy-flavour momentum correla- 
tion 0, 1^, that, although well defined in the massless limit, cannot be computed 
using the massless result. 

- Secondly, a massive calculation finds application in the area of fragmentation 
functions, where it can be used to improve the resummed cross section. 

- Then, this type of calculation helps to understand the relevance of mass correc- 
tions in the determination of from event shape variables and jets, and may 
be used to measure quark masses P, |^. 

- Finally, in the future accelerators, where high energies will be reached (Next 
Linear Collider), tt pairs will be produced and mass effects are very likely to be 
important in those measurements. 

In this thesis, we describe a Next-to-Leading-Order (NLO) calculation of the 
heavy-flavour production cross section in e~^e~ collisions, including quark mass ef- 
fects, for unoriented quantities. 
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Introduction 



At NLO, we get contributions both from two-, three- and from four-particle final 
states. The amplitudes describing these contributions contain infrared and ultravi- 
olet divergences. Therefore, some regularization procedure is needed. We use the 
dimensional regularization, that consists in changing the space dimensions d from 4 
to 4 — 2e. In this way, soft and ultraviolet divergences appear as poles in e. At the end 
of the calculation, the renormalization procedure takes care of the ultraviolet poles, 
while the infrared poles must cancel for infrared-safe physical quantities. 

Very recently, two calculations have appeared in the literature dealing with the 
same problem 0-0. They both use a slicing method, in order to deal with infrared 
divergences. This method consists to exclude, from numerical integration, the region 
of soft and coUinear singularities. For this reason, a cut-off parameter is introduced 
in the energy region of soft gluons and in the coUinear limit. In the excluded phase 
space regions, the transition amplitudes are replaced by their limiting values and 
the integration is done analytically, in order to obtain the coefficients of the infrared 
poles 1/e and that must cancel the corresponding terms in the virtual diagrams. 
Together with the poles, there is a finite part, which is integrated numerically over 
the remaining three-body phase space. Finally, the full transitions are integrated 
numerically in d = 4 dimensions over the four-body phase space, above the energy 
cut-off. Both results, the semi-analytical and the full-numerical one, depend on the 
choice for the value of the cut-off. Since this parameter was introduced arbitrarily, 
the final result must be independent from this value. 

The smaller this cut-off is, the better result is obtained: in fact, this is an approximate 
method, because it substitutes, in the singular region, the exact transition amplitudes 
with their limiting values. On the other hand, for small cut-off parameters, the 
numerical integration start to fail since the integration region gets closer and closer 
to the singular points. The best solution is found in the region of the cut-off values 
where the computed physical quantities become nearly independent from these values. 

In our work, we preferred to use a subtraction method, since, in this way, we do 
not need to worry about taking the limit for small cut-off parameters. Subtraction 
methods for the calculation of radiative corrections to e~^e~ — s>jets have been used 
in Refs. [T^, and they have also been successfully employed in the calculation of 
hadronic production processes. 

Together with the quoted massless calculations, parts of the massive calculation 
were computed by other groups: in the work of Ref. |TT[], a calculation of the process 
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e^e QQgg is given, but virtual corrections to the process e^e —>■ QQg are not 
included. In the same work, the amplitudes for two QQ pairs in the final state are 



calculated. In Ref. [|12[, the NLO corrections to the production of a heavy-quark 
pair plus a photon are given, including both real and virtual contributions, while, in 
Ref. [0 , the virtual contribution with a three- fermion loop is computed. 

In this thesis we deal with virtual and real contributions to the e^e~ annihilation 
process, that is with terms coming from the interference between one-loop diagrams 
with the tree-level terms and with final states characterized by a heavy couple of 
quarks plus a pair of gluons or light quarks. The computation of these contributions 
is the really hard part of the full calculation. 

This work is organized in a quite reverse order: in fact, we first give two physical 
applications of our calculation, leaving to the other chapters the detailed description 
of the full computation. 

In Chapter |I|, we introduce the heavy-fiavour momentum correlation, its relation- 
ship with the measurement of Rb, and the results we have obtained at leading and 
next-to-leading order. The computation of this quantity was the initial reason that 
conducted us to perform this calculation. 

In Chapter H, we revisit the fragmentation-function formalism for heavy quarks 
and the connection with our fixed order calculation. In fact, we check the validity 
of the initial conditions for the fragmentation functions, the NLO splitting functions 
in the time-like region and we derive an improved cross section, valid at NLO, for 
energies of order of the mass, but valid at Next-to-Leading-Log for energies much 
larger than the mass. 

The complete calculation is given in the following chapters: in Chapter we 
introduce the kinematical variables necessary to describe unoriented quantities, and 
we find the expression of the phase space in ci = 4 — 2e dimensions for the three-body 
final state (QQg) and for the four-body final states {QQgg, QQqq and QQQQ). 

In Chapter ^, we introduce the notation used during the calculation and we present 
the Feynman diagrams and the cut-diagrams which contribute to the differential cross 
section. All the results are in an analytic form and are implemented in a FORTRAN 
program. This program can compute each sort of unoriented infrared-safe shape 
variables and jet-clustering algorithms. This chapter is the hard part of the thesis, 
and some of the calculations involved, together with other details, are given in the 
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appendixes. 

The check of the cancellation of the infrared poles is performed in Chapter ^ 
where we apply the subtraction method. We first give a pedagogical introduction of 
this method before using it on the full cross section. In this chapter, we show how 
we can avoid to compute two-loop contributions and how we obtain the limits of the 
amplitudes in the soft and coUinear regime. At the end, we list some controls we have 
done: both internal consistency and external comparisons with the results of other 
groups. In fact, we found satisfactory agreement between our jet-clustering results 
and those of Ref . . 

In Chapter ^ we present some tables computed with different values of the ratio 
m/E, and we compare these results with the massless case. For this kind of calcu- 
lations, it would be difficult to perform analytical comparisons: for example, at the 
moment, this is impossible, since the other two groups have used a different 

method of computation. For this reason, we have chosen a few shape variables, and 
we have computed some moments, which can be obtained with a great precision. In 
addition, we present the results for the three-jet decay rate, computed according to 
four jet-clustering algorithms, for different values of the cut parameter. 

Finally, we summarize our work in the concluding chapter. 
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Momentum correlation 



In the early of 1996, the discrepancy between the measured and the theoretical value 
of Rb was drawing considerable attention in the physics community. We then began 
investigating on a possible systematic error, of dynamical origin, in the measurement 
of Rb. 

Very soon we realized that our first order result in ag was to be complemented 
by a second order calculation, so that we started the computation of the differential 
cross section for the process e'^e~ — > Z/7 — QQ + X, where Q is a massive quark, X 
is anything else, at order a^. 



1.1 Momentum correlation in Z/'y — > bb 

Rb is defined by the ratio 

R^-^, (1-1) 

J- had 

where F is the decay width of the Z/7 into the specified final state. In several 
experimental techniques for the measurement of Rb, the tagging efficiency is extracted 
from the data by comparing the sample of events in which only one b has been tagged, 
with the one in which both 6's are observed. If the production characteristics of the 
b and the b were completely uncorrelatcd, this method would yield the exact answer, 
without need of corrections. Of course, other correlations of experimental nature 
should be properly accounted for, but their discussion is outside the scope of the 
present theoretical section. Here wc deal with the standard QCD gluon emission, that 
generates a correlation of the quark- antiquark momenta of order a^. Other dynamical 
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effects, like the production of heavy-quark pairs via a gluon-sphtting mechanism, may 
affect the measurement. However, they are, to some extent, better understood: in 
fact, they tend to give soft heavy quarks, and they are, therefore, easily eliminated. 

We begin by considering the simple case in which the efficiency for tagging a B 
meson is a linear function of its momentum. This simplifying assumption allows us to 
make very precise statements about the correlation. Furthermore, it is not extremely 
far from reality, in the sense that the experimental tagging efficiency is often a growing 
function of the B momentum. For this reason, we introduce two kinematical variables 
to describe the momentum carried by the couple quark- antiquark 

_ P b(b) /, „\ 

^1(2) = ^^^^y ' (1-2) 
IP|b 

where p denotes the meson three-momentum. In terms of the usual notation, where 

2^1(2) = , (1-3) 

with E the total centre-of-mass energy and -E'b(b) the energy of the B[B) meson, we 
have 

We then assume that the tagging efficiency e is given by 

ei = Cxi, e2 = Cx2, (1.5) 

and that, in our ideal detector, the detection efficiency is not influenced by the pres- 
ence of another tag. We separate the events when one single B or B is tagged, from 
the events where both B and B are tagged. Denoting with A^^i the number of single 
tags, with N2 the number of double tags and with N the total number of events, we 
have 



N 



Rb (ei + 62) ^RbC ((xi) + {X2)) ^Rb2C {x) 



^ = Rb (ei €2) = Rb C {xi X2) = Rb C {xf X (1 + r) 



:i.6) 



where 

{x) = (xi) = {X2) , (1.7) 
and r is the momentum correlation 



1.1. Momentum correlation in Z/7 ^ 



Solving the system ( |0|) , we obtain 



Rh 



X (1 + r) 



(1.9) 



4NN2 

The quantity r cannot be measured, and therefore one has to compute it in order to 
determine Rf,. It is convenient to rewrite r in the following way 

^^((l-^.)(l-^..))-(l-x>-^ (1.10) 

{x) 

SO that we immediately see that terms, to whatever order in as, with only a bb couple 
in the final state give zero contributions, so that we do not need to compute two- 
loop diagrams at order (see Sec. ^]2| for a more detailed discussion). We can then 
compute this infrared-safe quantity with our program. Introducing the differential 
cross section, normalized to one, we define 



dxi dx2 



da 



dxi dx2 



'l-x) = —a + 0(a^^ 



2n 



fl.lll 



((1-Xi)(l-X2)) 



dxi dx'^ 



da 



a. 



where the integral is extended to the appropriate phase space region. Using eq. ( |1.1CI| ), 
we obtain, for the expansion in of r, 

'•=S''+(S)'(c+2a(,-a=)+0(a= 



;i.i2) 



Our results for r are displayed in Tab. where we have taken the total energy 
E = M^. From this table, we see that the coefficients a and c are both plagued by 



m 


a 


b 


c 


r 


1 GeV 


12.79(1) 


0.6628(1) 


155.5(3) 


0.1055a, + 0.23(1) 


5 GeV 


7.170(2) 


0.6182(1) 


50.94(5) 


0.0984 a, + 0.213(1)^2 


10 GeV 


4.858(1) 


0.5432(1) 


26.23(2) 


0.0865 a, + 0.2004(6)^2 



Table 1.1: Mass dependence of the coefficients a, b, c and r. The number in parenthesis 
(to be taken as zero if not given) represents the accuracy of the last digit. The total 
energy E is taken equal to M^. 



collinear divergences, because of the presence of large logarithms of the ratio E/m, 
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but the coefficients of the expansion of r itself do instead converge in the hmit of 
small mass. This is due to cancellation of these logarithms in the definition of r. 

It is easy to prove that this cancellation must occur to all orders in perturbation 
theory. In fact, according to the factorization theorem, we can write the double 
inclusive cross section for bb production, in the limit of E ^ m, as 

, ^'^ = / dyi dzi dy2 dz2 D{zi) D{z2) 6{yi zi - xi) 5{y2 Z2 - X2) , ^'^ , (1.13) 
axi ax2 J dyi ay2 

where da is the short- distance cross section, which, in the limit E/m^oo, has a 
perturbative expansion in with finite coefficients (observe that the distinction 
between the momentum- or energy-defined Feynman x becomes irrelevant in the limit 
we are considering here), and D{z) is the fragmentation function, that absorbs all the 
divergent terms. We then have 

{x)= fdx,dx2X,-^=([dzzDiz)) [dx,dx2X,-^ 
J dxidx2 \J / J dx\dx2 

(1,14) 

J (iX\(iX2 V-^ / (1X\(1X2 

In the ratio (0:1X2) / (a;)^ (and therefore in r) the integral containing D{z) cancels. 
Thus, the perturbative coefficients of r are finite in the limit E/m 00. Observe 
that, in the derivation of eqs. ( |1.14| ), we have assumed the relation 

D{z)dz = l, (1.15) 

which is appropriate when we can neglect the secondary production of bb pairs via 
gluon splitting. 

We also define the quantity r' 

J _ {xiX2 cut(l,2)) - {xf 

where the function cut(l, 2) assumes value one if the quark and the antiquark are in 
opposite hemispheres with respect to the thrust axis, and zero otherwise. We define 

2 



(1-Xi -X2 + XiX2Cut(l,2)) = ^6'+ (^g) c' + , (1.17) 



and obtain 



r' = ^b'+ (^) (c' + 2ab' - a^) + O (a^) . (1.18) 



271 \27rJ 
The quantities b' , d and r' are given in Tab. |1.2| . 



1.2. Double inclusive cross section at order a 
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m 


b' 


c' 


r' 


1 GeV 


0.3800(7) 


158.4(3) 


0.0605(1) + 0.12(1)^2 


5 GeV 


0.3653(6) 


50.9(1) 


0.0581(1) a, + 0.119(2) «2 


10 GeV 


0.3423(5) 


25.51(6) 


0.0545(1) a, + 0.133(1)^2 



Table 1.2: Mass dependence of the coefficients b', d and r' aX E = M^. 



1.2 Double inclusive cross section at order a. 



Up to now, we have assumed that the efficiency is hnear in the B momentum. Even 
in the more reahstic case in which the efficiency is a more comphcated function of 
the kinematical variables, it is possible to compute the inclusive cross section for the 
production of a bb pair, provided one also knows the B fragmentation function, which 
is, to some extent, measured at LEP. The appropriate formula is given in eq. ( |1.13| ). 
We limit our considerations up to order a^- The expression for the short distance 
cross section is 



da 



dyi dy2 



S{l-y^)5il-y2) + 



2a, 



yl + yl 



37r (l-yi)(l-2/2) 



where the "+" distribution sign specifies the way the singularities at yi 
1/2 = 1 should be treated: for any smooth function of yi and ?/2, we define 



(1.19) 
1 and 



2/1 +2/2 >1 



dyi dy2 



yl + yl 



[\-yi){\-y2. 
y\^yl 



G{yi,y2) 



G'(i/i,i/2)-G'(l,y2)-G(i/i,l) + G(l,l) . (1.20) 



dyi dy2 — r 

2/i+2/2>l (1 - ?/i)(l - 1/2) 

It is easy now to compute the value of r in the massless limit. Using eqs. (|1.8| ) 
and ( |1.14| ), we obtain 

r = — = 0.1061 X a, , (1.211 

3 TT 



that is consistent with the massive results of Tab. |1 . 1| . In addition, eqs. ( |1.13| ) 
and (|1.19|) give 



da 

dxi 



dxo 



da 



b{xi) . 



;i.22) 



dxi dx2 

The above equation fixes the factorization scheme to be the annihilation scheme as 
defined in Ref. W^. The choice of the scheme of eq. (|1.22D defines unambiguously the 
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result, without the need of computing exphcitly the virtual corrections. In fact, the 
most general formula for the short distance cross section at order a, is obtained by 
adding to eq. ( |1.19D terms of the form 



[fiVi) 5(1 - y2) + /(z/2) 5{l -yi)+g 5(1 - y{) 5{l - y^)] 



:i.23) 



where / is a generic distribution and g is a, constant. However, in order for eq. ( |1.22| ) 
to be respected, these terms must be absent. 



From eqs. ( |1.13D , ( |1.19D and ( |1.20| ), we immediately derive the O (a,) formula 



da 



dxidx2 



D(x,)D(x2) + 



2a, 



Svr Jo 



dyi dy2Q{yi + y2 - 1] 



yl + yl 



;i-yi)(i-2/2) 



\y^. 

' Xi\ 



2/2/ 2/12/2 



0(2/1 - xi) 0(?/2 - X2) + b{xi) b{x2) 



-D\-] b{x2) - Q{yi - xi) - b{x,) bi"^]- 0(2/2 - X2) 

.2/1/ 2/1 \2/2/ 2/2 



:i.24) 



where is the Heaviside function. As an illustration, we plot in Fig. |1.1| the double 
inclusive cross section da/dxi dx2 as a function of xi, for several values of X2. We use 
the Peterson parametrization of the fragmentation function 



D(z) 



N' 



z(l-z) 



2 ' 



:i.25) 



where A*"^' is fixed by the condition / dz D{z) = 1. We took the values e = 0.04, which 
gives (x) = 0.70, and a, = 0.12. The positive momentum correlation is quite visible 
in the figure. As X2 increases, the peak of the distribution in xi also moves towards 
larger values. Using the above formula, we can compute r again. The result should 
not depend upon the choice of the fragmentation function. In this case, the quantity 
(x) does not receive corrections at order a,, and we get 



(x) = 0.7036 , (xi X2) = 0.4950 + 0.0525 x a, , 



0.1061 X a. 



;i.26) 



which is consistent with the value of r previously obtained (see eq. (|1.21|) ). 



1.3. Final results 
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Figure 1.1: Double inclusive cross section da/dxi dx2, plotted as a function 
of xi for several values of X2- 



1.3 Final results 



Assuming as{M^) = 0.118 (corresponding to the Particle Data Group average |[T5[| ) 
and m = 5 GeV, we have, at leading order, r = 0.0984 x = 0.0116, and at next- 
to-leading order r = 0.0984 x + 0.213 x = 0.0146. With the assumption of a 
rough geometric growth of the expansion, we can give an estimate of the theoretical 
error due to higher orders: r = 0.0146 ± 0.0007. In addition, we get r' = 0.0064 at 
leading order, and r' = 0.0083 at next-to-leading order. 

As a conclusion, we can say that, as far as its perturbative expansion in powers of 
as is concerned, the average momentum correlation is a quantity that is well behaved 
in perturbation theory, and it is also quite small. Since its effect is typically of the 
order of 1%, one may worry that non-perturbative effects, of order A/E (where A 
is a typical hadronic scale) may compete with the perturbative result. This is a 
very delicate problem, since we know very little about the hadronization mechanism 
in QCD. In Ref. |]TB[, this problem was addressed in the context of the renormalon 
approach to power corrections. It was shown there that corrections to the momentum 
correlation are at least of order (A/i?)^, and thus negligible at LEP energies. Although 
this result cannot be considered as a definitive answer to the problem, it is at least 
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an indication that power corrections to this quantity are small. 



Ckapter 2 



Fragmentation functions 



In this chapter, wc introduce the notion of fragmentation functions for massive quarks. 
Using our O (af) calculation of the differential cross section for the production of 
heavy quarks in e~^e~ annihilation, we verify that the Leading (LL) and Next-to- 
Leading Logarithmic (NLL) terms in this cross section are correctly given by the 
standard NLO fragmentation-function formalism for heavy-quark production, in the 
limit m/E ^ 0. 



2.1 Fragmentation functions for keavy quarks in e^e colksions 

The inclusive heavy-quark production is a calculable process in perturbative QCD, 
since the heavy-quark mass acts as a cut-off for the final state coUinear singularities. 
Thus, the process 

e+e- ^Z/j^Q + X , (2.1) 

where Q is the heavy quark and X is anything else, is calculable. Its cross section 
can be expressed as a power expansion in the strong coupling constant 

J oo 

— {x, E,m)^Yl E, m, //) «:(//) , (2.2) 

n=0 

where E is the centre-of-mass energy, m is the mass of the heavy quark, /i is the 
renormalization scale, and 

«.(/^) = ^ ■ (2.3) 
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As usual we define 

X = — — , (2.4) 

where q and p are the four-momenta of the intermediate virtual boson and of the 
final heavy quark Q. We define the heavy-quark fragmentation function in e^e~ 
annihilation as 

1 dcf 

D{x,E,m) = —{x,E,m). (2.5) 

o"tot dx 

Since each heavy quark in the final state contributes to the fragmentation function, 
its integral with respect to x gives the average multiplicity of these quarks 

{n{E,m)) = dxD{x,E,m) . (2.6) 
Jo 

When E/m is not too large, the truncation of eq. ( |2.2| ) at some fixed order in the 
coupling ds can be used to compute the cross section. On the other hand, if E ^ m, 
since the n^^ order coefficient of the expansion has the form 

a(")~ log- , (2.7) 



m 

we cannot truncate the series. In fact, if 

«,log-^l, (2.8) 
m 

each term of the series (p.2| ) is of the same order of magnitude of the first one, so that 
we cannot trust a fixed order calculation, that computes only a finite number of terms. 
We must then try to resum these large logarithms. The resummation procedure is 
described in Ref. [0, and the resummed differential cross section that is obtained 



with this procedure has the form 



da 

-(x,E,m) 



= J:/3^''\x) (asiE)\og-) +asiE)J2^^-\x) «.(E)log- 

+ al{E)Y.^6^-\x) (a.(E)log-j + . . . + O i^—j , (2.9) 

where terms that are suppressed by powers of m? / E"^ are not included, because they 
can be neglected if compared with the large log (i?/m), and where we have taken 
fi = E, for simplicity of notation. If we compute, with the resummation procedure, 
the first series of eq. ^.9| , we talk of Leading-Order (LO) resummed cross section; if 
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we include the next series, we refer to it as Next-to-Leading-Order (NLO) resummed 
cross section, and so on. 

In the approximation when you can neglect powers of rn? / E"^, it can be observed 
that the inclusive heavy-quark cross section must satisfy the factorization theorem 
formula 

da , „ , ■s-^ /"^ , , X daj 



{x,E,m) = dydz6{x -y z) {z,E,fi) D,{y,fi,m) 

Jo dz 

Af-,/^,m) , (2.10) 



dx 



z dz \z 

■L 

where dai{z, E, jj)/dz are the MS-subtracted partonic cross sections for producing 
the parton i, and Di{y, fi,m) are the MS fragmentation functions for the parton i 
into the heavy quark Q. In order for eq. (|2.10| ) to hold, it is essential that you use 
a renormalization scheme where the heavy flavour is treated as a light one, like the 
pure MS scheme. Thus dai{z, E, fi)/dz has a perturbative expansion in terms of as 
with Uf flavours, where nf includes the heavy one. 

In order to compute the differential short-distance cross section that describes the 
process Z/'y -^i + X, we must introduce a regulator, because we now deal with a 
massless quark (no mass dependence in the expression of dai), so that we have to 
face the presence of coUinear singularities. Choosing the dimensional regularization 
{d = 4 — 2e), the MS prescription amounts to throw away all the terms that contain 
poles in e in the expression of dai. After that, the form of Di is flxed, because its 
convolution with d&i must give the physical, flnite, differential cross section da. 

With the factorization theorem, we succeed in separating the two scales of energy 
E and m into two different terms, through the introduction of a factorization scale 
H, that we take equal to the renormalization scale, for simplicity of notation. The 
scale II should be chosen of the order of E, in order to avoid the appearance of large 
logarithms of E/fi in the partonic cross section. 

The MS fragmentation functions Di obey the Altarelli-Parisi evolution equations 

dDi , , dz „ ( X 

dXogy? J 

that resum correctly all the large logarithms. 

The Altarelli-Parisi splitting functions Pjj have the perturbative expansion 

P,,(x,a,(/i)) =a,(/i)P,f (x) + a2(;x)4')(x) + 0(«^) , (2.12) 



{x,ii,m) =J2 [ — ^ij f->«s(/^)) Dj{z,ii,m), (2.11) 

■ J X Z \Z J 
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where Pj^"^ are given in Ref. and P^j' have been computed in Refs. |2D]-|2^. The 
only missing ingredients for the calculation of the inclusive cross section are the initial 
conditions for the MS fragmentation functions. These were obtained at NLO level in 
Ref. [0 by matching the O (a^) direct calculation of the process (i.e. formula ( p.2| )) 
with the expansion of eq. ( |2.1CI| ) at order a^. They have the form 



(2.13) 



bQ{x,fio,m) = 6{1- x) + asifio) SqIx, fio,'^) + O (af 
bg{x,fio,m) = asifio) df\x,fio,m) + O (afj , 

all the other components being of order a^. Thus, to compute the NLO resummed 
expansion, one takes the initial conditions (|2.13|) , at a value of /iq of order m (so that 
no large logarithms appear), evolves them at the scale (taken to be of order E), and 
then applies formula ( |2.10| ), using a NLO expression for the partonic cross section 



d^\x) + d^^\ 



X, E, jj) asifi) + O [a 



For example, if the parton i is the heavy quark itself, one gets 



dx 



[x, E, fi) = (5(1 — x) + dQ\x, E, /i) asifi) + O {o?^ 



and if the parton i is a gluon 

da, 



dx 



X, E, /i) = d^g\x, E, /i) asifi) + O {afj 



(2.14) 



(2.15) 



(2.16) 



where we have normalized the cross section to one at zeroth order in the strong 
coupling constant. 

The procedure outlined above guarantees that all terms of the form (osL)" (lead- 
ing order) and as(«s-Z^)"' (next-to-leading order), where L is the large logarithm, are 
included correctly in the resummed formula. Notice that, at NLO level, the scale 
that appears in as, in eqs. ( |2.13|) and (|2.14| ), could be changed by factors of order 1, 
since this amounts to a correction of order af. However, one cannot set yUo = in 
eqs. ( |2.13| ) (or /i = m in formula ( p.l4| )), since this amounts to a correction of order 
a^L, and thus it would spoil the validity of the resummation formula at NLO level. 

The validity of this procedure has however been questioned by the authors of 
Ref. |]23|. In their procedure, the heavy-quark short- distance cross section is replaced 
by 

da'n 



dx 



(a:, E^ jj) = (5(1 — x) + as{E) dg^ix, E, jj) + d^Q{x^ /iq, m) 



+ 0{at) , (2.17) 
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and the initial condition by 

D'q{x) = 5{l - x) , (2.18) 

which is to be evolved from the scale /io to the scale /x using the NLO MS evolution 
equations. This procedure differs at the NLO level from the standard procedure 
advocated in Ref. Ijl^. The difference starts to show up in the terms of order afL. 

We have used our calculation of the differential cross section at fixed order a? to 



check the standard formalism and thereby to dismiss the approach of Ref. ||2^. By 
the way, the full consistency of the two results gives support to the correctness of our 
computation. 



2.2 Calculation 

Instead of dealing with the realistic case of Z/7 decay, we perform the calculation for 
a hypothetical vector boson V that couples only to the heavy quark with vectorial 
coupling. 

We introduce the following notation for the Mellin transform of a generic function 

f{N) = I dxx^-^f{x) . (2.19) 
Jo 

We adopt the convention that, when N appears, instead of x, as the argument of 
a function, we are actually referring to the Mellin transform of the function. This 
notation is somewhat improper, but it should not generate confusion in the following, 
since we will work only with Mellin transforms. 

The Mellin transform of the factorization theorem ( p. 101) is given by 



a(iV, E,m) = J2 E, /z) A(iV, /x, m) , (2.20) 

i 

where 

aiN,E,m)= dx x^~^^(x, E,m) , (2.21) 
JO dx 

and a similar one for ai{N, E, n); the Mellin transform of the Altarelli-Parisi evolution 
equations ( p^.H] ) is 

dDi{N, /i, m) 



d log /i^ 



Y^asifi) PS\N) + P,^P{N)a,ifi) + 0(al) D,{N,fi,m). (2.22 
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In order to make a comparison with our fixed order calculation, we need an expression 
for a{N, E, m) valid at the second order in a^. Thus, we solve eq. (p.22| ), with initial 
condition at /i = /io, accurate at order a^. This is easily done by rewriting eq. (|2.22| ) 
as an integral equation 

Di{N,fi,m) = Di{N,^o,m) 

+ rdlogfi'^ a^ifi') MP{N)+P,f{N)a,{fi')] D^{N,fi',m) . (2.23) 



J 



The terms proportional to can be evaluated at any scale (/i or /io), the difference 
being of order a^. Factors involving a single power of can instead be expressed in 
terms of as{f^o) using the renormalization group equation 

as(/i') = a.sifJ'o) - 2 -K bo alifio) log — + C (a^ (/io)) 

^0 (2.24) 

, UCA-^nfTp 

"0 ~ voZ ' 

iZ n 

with Uf the number of fiavours, including the heavy one. Equation ( p.23|) then be- 



comes 



A(iV,/i,m) = A(iV,/io,m) + 5] rc?log/i''«,(/io)P,f (iV)4(iV,/i',m) 
+ Y: o?M P^\n) D,{N, /io, m) log 4 



- 2 TTfeo E «'(/^o) p!:-\N) D,{N, /io, m) ^ log^ ^ ■ (2.25) 

Now, we need to express Dj{N, fi' ,m) on the right-hand side of the above equation 
as a function of the initial condition, with an accuracy of order Og. This is simply 
done by iterating the above equation once, keeping only the first two terms on the 
right-hand side. Our final result is then 

,2 



A(iV, /i, m) = Di{N, /io, m) + J2 as(/io)4°^(iV) A /^o, m) log 



2 



+ E «s(/^o)Pr(iV) Pi-^(iV) /^.(iV, /^o, m) - log^ 



2-° /^^ 



+ E«s(/^o)4'^(iV)4(iV,/io,m) log 4 
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19 



(2.26) 



Since the initial condition is 



Dj{N,fio,m) = 6jQ + asifio) df\N,fio, m) + O [a'i [fiQ 



(2.27) 



eq. (|2.26|) becomes, with the required accuracy, 



Di{N, fi, m) = 5iQ + asifio) 4^\n, /iq, m) + as{no)PiQ{N) log ^ 



>(0), 



1^^ 



+ E «'(/io)P.f (N) df (AT, ^0, m) log 



+ E«.'M4°^(iV)pig(iV)iiog^4 



+ «s'(/^o)Pi^ (iV) log 4 - ^ ^0 «^(/io)P/Q^(iV) log^ 4 
/io /io 



Re-expressing (/xq) in terms of (/i), using eqs. ( p.24|) , we get 



(2.28) 



Di{N,fi,m) = SiQ + as{fi)dl^\N,fio,m) + 2 71 bo ai{fi)d\'''{N, fiQ,m) \og^ 

+ cys{fi)P^SiN) log 4 + E /^o, m) log 4 

+ E n^g(iV) ^ log^ 4 + log 4 



+ 7r6o«,^(/x)P.g^(iV) log^^. 



2 



(2.29) 



The partonic cross sections are given by 



a,{N, E, fi) = 6,Q + % + a,(/i) df\N, E,fi) + [aUfi)) , (2.30) 

where O'^"* vanishes unless i is either Q, Q or (/. Thus, combining eq. ( p.30|) with 
eq. (|2.29|) , according to eq. (|2.2(]| ), we obtain 

ag) (AT, E, fi) + rfg) (AT, fio, m) + P§l{N) log ^ 

f^o. 

{2 2 
E 4'^ (N, E, fi)Pi'^ log ^ + 2 vr 60 /^o, ^) log ^ 



a{N,E,m) = l + as{fi) 
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+ E [41 (^) + (^)] 4^ /^o, m) log 



+ EK2(^) + 4°i(^) 



^S(^)^ log's 



2 /f: 



+ 



PqqW + P^^N)] log < + ^ ^'o P^liN) log' ^ 



(2.31) 



The above formula should accurately describe the terms of order a^L, a^, a^L and 



a 



L'. Terms of order a^, without logarithmic enhancement, are not accurately given 
by the fragmentation formalism at NLO level, and have consistently been neglected. 

The lowest order splitting functions are given by 

1 



P^'^N) 



,(0).,.. ^ \2 + N + N'' 



N{N^ - 1) 
2 + N + N'' 



_N{N +1){N + 2) 
where, restricting ourselves to integer values of N, 

N ^ 

j=l J 



(2.32) 



(2.33) 



The splitting functions Pqq{N) and P-^^{N) are given by 



(1), 



(2.34) 



where the non-singlet components are 

p^m) = psm + psm) + pqW) 



P^^iN) = C'APf{N) + A{N)] 
Pg^iN) = '^Ci.CaPg{N) 
Pqq{N) = nfCpTpPMN) 



P^'oiN) - P^^oW + P^^iN) 



Ca I 



P^Q^N) ^ CIPa{N) 



P^niN) = -\Cj.CaPa{N) 



(2.35) 
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and 



Pf{N), A{N), Pg{N) and Pnf{N) were taken from the appendix of Ref. fT^ and 
Pa{N) is given in eq. (5.39) of Ref. [^. We have obtained our exphcit expression for 
Pqtq{N) using the relation 

r^f . .... .... 

(2.37) 



Pq'q 



pNS 

2nf 



pNS 
QQ 



where Pqq is the singlet component^, calculated in Ref. ||2^. Equation ( |2.37| ) is easily 
seen to follow from eqs. ( |2.34| ) and from eqs. (2.42) of Ref. [0 . 

The expressions for dq^ and dq^ are respectively given in eqs. (A. 12) and (A. 13) 
of Ref. []T^ . The coefficient a^^^ can be obtained by performing the Mellin transform 



of the expression CT,g + CL,g, where CT,g and CL,g are given in eq. (2.16) of Ref. [|14 
Thus 

2 2 1 

(A^- 1)2 ^ iV2 " (A^ + 1)2 



d^^\N,E,fi) = Cf< 



2(2 + N + N^)^ 

Mog — + 4 

/i2 



Ar(Ar2 
2 



A^- 1 



1) 



1) _ —SAN) + 
' N ' N + 



-5i(Ar+ 1) 



(2.38) 



In order to make a more detailed comparison with our fixed order calculation, we 
separate the O [p?^ contributions to o-{N,E,m) according to their colour factors. 
Choosing for simplicity fi = E and /iq = m, and using the notation 



«Q/.(^) = a'Q'/,iN,E,f,)\,=E , d'^'^iN) = d\^>^{N,f,o,m)\ 

we write 

a{N, E,m) = l + a,(E) A{N, E, m) + al{E) B{N, E, m 

with 



A) 



/iO=ni 1 



(2.39) 
(2.40) 



A{N,E,m) 
B{N,E,m) 



d^^\N) + d)^>{N) + Pl,^>{N)log^^ 



(0), 



E' 



(2.41) 



Bc^ {N, E, m) + Bc^ {N, E, m) + B^,. {N, E, m) + Bt^ {N, E, m) 



^We warn the reader that, sometimes, in the Hterature, the notation is used for the "sea" 
component, and Pqq is used for the singlet one. Here we use Pqq for the full QQ splitting function. 
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and 



Bc,{N,E,m) 



[p^qIW [rfg)(iV) + a«(iV)] + P^-{N) + P|-(iV)} loj 
1 



^2 



+ 



J? 



11 



7(1). 



log 



BnAN,E,m) 



'-nfTFd^^\N) 



log — 



^n^r^P^'g(iV)log2^ 

>(0) , o 



^2 



^2 



W PSi + 2 Pi°J (AT) (iV) + 2 P,,, {N)\ log — 



(2.42) 



where the subscripts Cp, Ca, nj and Tp denote the Cp, CfCa, njCpTp and CpTp 
colour components. 

Our fixed order calculation can be used to compute the cross section for the 
production of a heavy-quark pair plus one or two more partons, at order a^. We 
separate contributions in which four heavy quarks are present in the final state, from 
those where a single QQ pair is present together with one or two light partons. These 
last contributions are computed only in a three-jet configuration, and they are singular 
in the two-jet limit, that is to say, when x— >-l. Furthermore, the virtual corrections 
to the two-body process V Q + Q are not included in our calculation. In order to 
remedy for these problems, we proceed in the following way (see Sec. ^.2| for a detailed 
description of the method). The O (af) inclusive cross section ioi V ^Q + Q + 
can be written, symbolically, in the following form 

da 



= a(°)(5(l -x) + asj dY a^^^ (x, Y) + J dY af^ (x, Y) + 2jdYi 



,(2)/ 



(2.43) 

where Y denotes all the other kinematical variables, besides x, upon which the final 
state may depend (see Chapter |^). We assume fi = E, and we do not indicate, for 
ease of notation, the dependence upon E and m of the various quantities. The term 



2.2. Calculation 



23 



a] arises from final states with a single QQ pair plus at most two light partons, 
while a)j arises from final states with two QQ pairs. The factor of 2 in front of the 

(2) 

a\ contribution takes account for the fact that we may detect either one of the two 



heavy quarks, as is illustrated in eq. ( |2.6| ). The moments of the inclusive cross section 
can be written in the following way 

a{N) = J rfxx^-^^ = a + asJdxdY (x^-^ - l) a«(x,F) 

+ al J dxdY [x^-^ - l) al^\x,Y) +2 J dxdV (^x^-^ - a^h\x,Y) 

(2.44) 

where 

a = a(°) + asjdx dV a^^^ (x, Y) + J dx dY af ^ (x, Y) + J dx dY a}f ^ (x, Y) . 

(2.45) 

The expression for a{N) can now be easily computed with our program, since the 
(^x^~^ — 1^ factors regularize the singularities in the two-jet limit, and suppress the 
two-body V ^ Q + Q virtual terms. Furthermore, in the massless limit, 

a = l + 2as{E) + cal{E) +o(^^'j+0{afj , (2.46) 

where c is a constant. In fact, the O (a^) term does not contain any large logarithm, 
as long as ag is the coupling with n/ flavours, including the heavy one. If, instead, 
the cross section formulae are expressed in terms of ai""^ ^\ we have to take account 
of the different number of flavours. From the renormalization group equation (|2.24|) , 
that we rewrite remarking the number rij of flavours, 

= "i-^H/^o) - a'M log^^+O {a^ (/io)) , (2.47) 

6 K ^ ^ 

and the matching condition 

ai"/)(/.o) = ai"/"^n/io), (2.48) 

we derive 

a("/)(/i)-ai'^/-i)(/i) = |Tp«^og4 + (^(«sM) . (2.49) 
In this way, the total cross section of eq. ( p. 461) becomes 



a = l + 2a'^:^f-'\E) + ^Tpa'^{E)\og^ + 0(al) . (2.50) 
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Reintroducing the energy and mass dependence in eq. (|2.44 ) and using eq. ( 2.46 ), we 
have 

a{N,E,m) = 1 + as{E) L{N, E,m) + a^^{E) M{N, E,m) + ca^,{E) 



2' 




Oia'A , (2.51) 



where 

L{N,E,m) = 2 + j dxdY (x^-^ -l)a^^\x,Y) 

M{N, E,m) = J dx dY {x^~^ - l) af ^ (x, Y) 

+ 2 JdxdY (^x^-^ - a^h\x,Y) . (2.52) 

We have calculated L{N, E, m) and M{N, E, m) numerically, using E = 100 GeV 
and m = 8, 4, 3, 2.5, 2, 1.5, 1, 0.6, 0.5, 0.4, 0.2 GeV, for a vector current coupled 
to the heavy quark. We expect that, for small masses, A{N, E, m) of eqs. ( |2.41| ) 
should coincide with L{N, E, m), and M{N, E, m) should differ from B{N, E, m) by 
a mass and energy independent quantity, since such term is actually beyond the 
next-to-leading logarithmic approximation. We find very good agreement between 
A{N,E,m) and L{N, E,m). We present the results for M{N,E,m) separated into 
the different colour components 

M{N, E, m) = {N, E, m) + Mc^ {N, E, m) + (A^, E, m) + Mt^ {N, E, m) . 

(2.53) 

In Fig. ^]l] we have plotted our results for Mc^ (crosses with error bars) and for Bc^ 
(solid lines). An arbitrary (A^-dependent) constant has been added to the curves for 
Bcj^, in order to make them coincide with the numerical result for m/E = 0.015. We 
find, as the mass gets smaller, satisfactory agreement for all moments. Notice that for 



higher moments we need smaller masses to approach the massless limit. In Figs. 



0| and |2.4| we report the analogous results for the remaining colour combinations. 



Again, we find satisfactory agreement. 



If one follows the procedure proposed by the authors of Ref. ||2^, eqs. ( ^.42[ ) are 



modified in the CpCa and in the njCpTp coefficients. More specifically, the terms 
proportional to d^Q all disappear from the expressions of the CfCa and of the npCpTp 
coefficients. In fact, by inspecting formulae ( f^.l7D and ( |2.18| ), and the derivation of 
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Figure 2.1: CfCa component of the coefficient in a{N, E,m), as a 
function of log E'^/m?, for N =2, 5, 8 and 11. The dashed lines correspond 
to the alternative approach proposed by the authors of Ref. J2^. 



eq. ( |2.31| ), we see that the only difference between the two approaches is that the 
term (/iq) dg^ is replaced by as{ij) (ig^ which, using the renormalization group 
equation (|2.24|) , amounts to a difference of —2TxhQa1d^Q log/i^//iQ, precisely what is 
needed to cancel the term of the same form appearing in eq. (p.31|) . The modified 

(dashed lines). It is quite clear that the 



result is also shown in Figs. and |2] 
approach proposed by these authors does not work. 



2.3 Improved cross section 



We can get an improved cross section by merging the fixed order calculation with 
the NLO resummed cross section, to obtain a formula that, for m, is accurate 
to order a^, and, for ^ m, is accurate at NLO level. This can be accomplished 
with the following considerations: we have computed the differential cross section till 



26 



Ckapter 2. Fragmentation functions 




order so that we know the first three coefficients of the series ( |2.2| ), that is 
da 



a^^\x,E,m) + a^^\x,E,m)a,{E) + a^'^\x,E,m)al{E) , (2.54) 



fix 



where we have taken ^ = E, for simplicity of notation. In this equation, we have not 
computed the contribution to a^^-* in the two-jet region. This is not a hmiting point 
since an approximate expression for a'-^^ till order O (^{rn? / E"^)^^ is now available 

The NLO resummed cross section is given by (see eq. 



da 

^{x,E,m] 



n=0 



Y.f3^''\x) a,(E)log- +a,{E)Y.j^-\x) a,(E)log- 



n=0 



(2.55) 



so that, the improved formula reads 

a^°\x,E,m) + a'^^\x,E,m) as{E) + a^^\x, E,m) al{E) 

E' 



da 

^{x,E,m^ 



imp 



+ ^/?(")(a^) «s(i?)l0g 



n=3 



m 



°° / E\ 

+ asiE)J2^^-\x) a,(i?)log- 
^2 V mJ 



(2.56) 
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Figure 2.3: Same as in Fig. 2J, for the Cp component. 



The actual way this improved formula is obtained is a bit more complicated, because 
we do not have eq. ( p.55|) in this form, but we have a numerical function of cxs- We 
start from the following considerations 



lim a^^Hx.E, 



m— i>0 



lim a^^\x, E, m) 

m — ^0 



lim a^'^\x, E, m) 

m — ^0 



X] 



E 



pW^x) log- + 7(°)(x) 



m 

(3^ (a;) iog2 - + ^(1) (x) log - + (x) 



(2.57) 



which allow us to write 



da 



fix 



\ da 

hm < — (x, E, m) 
m^o \ dx 



fix> 



da 

+ —ix,E,m) 
dx 



a'^^\x,E,m) + a^^\x,E,m) as{E) + a^^\x, E,m) al{E) 

CO , F \ " 

+ «.(E)log- 



+ a,iE)J2c^-\x) (a,iE)\og-) - 5^'\x) a'^{E) . (2.5. 

n=2 ^ 
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2.4: 


Same as in Fig. 2.1, for the 


CpTp component. 



By comparing eq. ( p.56| ) with eq. ( p.58| ) we obtain 



da 
dx 



[x, E, m] 



imp 



da 
dx 



[x, E, m] 



fix 



{da 

— lim < — (x, E, m] 

m^O dx 



fix> 



da 

+ —[x,E,m) 
dx 



+ 6^'\x) al{E) , 
where S^^^ (x) is defined by the hmiting procedure of eq. (|2.57| ) 



(2.59) 
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Kinematics 



3.1 Kinematics and four-tody pkase space witk two massless 
particles 





Figure 3.1: Four-body kinematics. 



The four-body processes we are considering are illustrated in Fig. and sum- 
marized by 

e+(p'J+e-(pe) ^ Zhiq) ^ Q(p) + Q(p') + g(k) + g{l) 
e+ (p'J + e- (pe) ^ (q) ^ Qip) + W) + Qik) + , 

where Q is the massive quark, q is the massless quark and the momenta satisfy 



P = k^ = 



2/2 2 

p = p = m , 



9^ = (P'e + Pe 



(3.1) 



We hope that no confusion arises between the total momentum q and the light quark 
q. In the centre-of-mass system of the two massless particles, we can express the 
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(i- mo merit a in the following way 

I = Iq {1, . . . ,sm9 sm(j),sm9 cos(f),cos9) 

k = fco (1, . . . , — sin 6' sin 0, — sin6' COS0, — cos( 



P 



P 



Po 1,...,0, 



\ 



Po 



j2 sm a, 



\ 



Po 



cos a 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



where the dots indicate d — 3 equal and opposite components in the expression for / 
and k, and d — 3 zeros in the expression for p and p'. 

To describe the unoriented four-body phase space, we need five independent vari- 
ables, which we choose to be 



Xi 



2q ■ p 



X2 



2q-p' 



{k + lf 

y = — 5 — ' 



9, 



q^ q^ 
In the centre-of-mass system of the light particles 



lo — ko 



(3.6) 



(3.7) 



where we have used the definition of y of eq. (|3.6| ). From momentum conservation, 
written in the following way. 



(g-p)2 = (p' + fc + /)2 

iq-p'f = ip + k + lf 



we have 



and from 



Po 



l-X2-y 



Po 



{p + p +k + lf = q^ 



1-Xi-y 



(3.8) 

(3.9) 
(3.10) 



we obtain 

Po + p'o + \/(fVy 



Pq — j sin^ a — [ ^Pq — + ^^Jpq — cos a 

We can solve this expression to give 

y{p-xi- X2) + (1 - a;i)(l - X2) - y'^ 



j2 



COS a 



l-xi- yf -Apy J{l-X2-yf - Apy 



(3.11) 
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where 

p = -r- (3-12) 

Starting from the four-particle phase space in d dimensions 

f d^-^p d'^-'p' d'^-H d^-^k / , M 

y 2po(27r)°'-i 2p'o(27r)^-i 2/o(27r)^-i 2fco(27r)'^-i ^ ; 

(3.13) 

we introduce the following two identities 

/ -l-k) = l, / ^ 2" - = 1 , (3.14) 

which allow us to integrate eq. (|3.13|) in to, and to obtain 

2 fdy , d'^-H d'^-^p d'^-^p' ,^ 

. ' 

X -Ar—^ _ {271^ 5\t -k-l), (3.15) 



27r 



2/o(27r)^-i 2ko{27Ty 



d*2 



where 



to = V|t|' + (3-16) 
and where Q{y) is the Heaviside function. In this way we succeed in dividing the four- 
body phase space into two simpler Lorentz-invariant scalars, that we can evaluate in 
the most appropriate reference system. 

We compute c?$2 in the centre-of-mass system of the two massless particles. In 
this system 

t = k + l = 0, to = ko + lo = J^y/^ , (3.17) 



where we have used eq. ( p.7|) . Integrating in k, we get 



4/2(27r)^-2^' J All{2TiY- 

where we have introduced the integration over the solid angle f2 in c? — 3 dimensions: 

n'^= I dn'^ = . (3.18) 
J r(f) 
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We perform the following change of variables 

ly = Iq sin 9 cos (p 

h = /ocos^ (3.19) 



1^1 = ^/|l|2-/2_;2 = /psin^sin0, (3.20) 



with Jacobian 

l2 



dly dlz d\l±\ = Iq sin 6 dO dcj) dlo , (3-21) 
and considering that |1| = Zq; we can write 



d-3 



r dO d(h din , r, ^ , ^ .xW-a 2 7r 2 

The integration over Iq is straightforward and, with the last change of variable 

= i(l -cos^) , (3.23) 



we have, in = 4 — 2e dimensions 



d^2 = ^ [q^y]'' fdv [v{l - v)r^ r#(sin0)-^^ , (3.24) 

STT i (1 — ej ^ ^ JO JS^ Jo 



where 



= I # (sin0)-2^ = 4% ^^|^_^^'| (3.25) 



and where we have used eq. ( ^.17| ). 

We can immediately integrate in t the three-body phase space of eq. (p.l5|) to 
obtain 

f 271 d'^-^p d'^-^p' ^, , , , , 

^^^ = 7 27; 2M^ 2p'q(2.Y'^ - to - Po - Po) (3.26) 

with the condition 

q = t + p + p' . (3.27) 

We evaluate this integral in the laboratory frame, where 

g=(go,0), t = -p-p'. (3.28) 

We orient our reference axes in such a way that p is along the z axis and p' belongs 
to the yz plane, forming an angle (3 with p. In this system 

d'^-^p = |p|^-2^|p|l]'^-i 

rf'^-y = dp', d'^-^p'^ = dp', \p'^\'^-^d\p'^\Q'^-^ = Ip'I d\p'\ rf/3|pi|^-3 Q'^-^ 

= Ip'f-^dlp'l {sin f3f-'' d(3n''-\ (3.29) 
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where we have used 



Pz 



Ip'j 



I p' I cos [3 
I p' I sin /3 , 



so that eq. ( ^.261 ) can be written 



2n |p|'^"2^|p|r]'^-i |p'|'^-2^|p'|^]'^-2(sin/5)'^"^rfcos/5 



2to 2po(2vr)'^-i 2p'o(27r)^-i 
From eqs. ( 3.16|) and (p.28|) we have 



^(go -to-Po -Po) ■ 
(3.30) 



to = V IpP + IpP + 2 IpIIp'I cos/3 + q'^y 
and we can integrate in f3 to obtain 



dcosP 



TT 



pI'^-^ ci|p| n^-^ |p7-3 d\p'\ n^^-^ (sin/3) 



iPllP 



v(i-4 



2po(27r)^-i 2p^(27r)'>'-i 
with /? determined by the 5-function argument 



qo -po~p'q = vIpP + IpP + 2 IpIIp'I cos/3 + q^y . 

If we evaluate Xi and X2 of eq. (|3.6| ) in the laboratory system, we get 



Xi 



2po 



Qo 



X2 



and considering that 



2 2 1(2 \ 



J 1 2 „/2 



/2 2 y / 2 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



we can write (i$3 as 



il'^Y ^''J dxi dx2 (sin/3) (^x^ — (^Xg — 



r(2 - 2e) 2 (47r)3 



r(2 - 2e) 2 (4 
From eq. ( |3.33|) , we have 



cos /3 



{xl-4y) -{xl- p) - (xl- p) 



(3.36) 



(3.37) 
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where we define 

= 2 — xi — X2 , (3.38) 

and the integration range for xi and X2 is determined by the reality condition for 

cos P 

-l<cos/3<l. (3.39) 
Inserting the expressions of eqs. (|3.24|) and ( p.36|) into eq. ( |3.15|) , we obtain 

" (47r)2r(l-6) {^J Hjdx,dx2dyy 'e{y) 

X {4 {xj - p) {xl - p) - [{xl -4y)-{xl-p)- {xl - p) 

X f dv[v{l-v)]-^^ rd<j> {sin 

Jo 1^6 Jo 



(3.40) 



with 



H = 



IGtt 



2e 



(3.41) 



r(2 - 2e) 2 (47r)3 \ ) ' 

The last step to perform is the determination of the integration region. First of 
all, the physical boundaries for the integration variables are 

v/p<xi,X2<l, Q<y<{l-^f . (3.42) 

If we choose to integrate first in we solve eq. ( |3.39| ) with respect to y, and we get 

y-<y<y+. (3.43) 

with 



y± = 



4 



x„ 



{xl - p) - {xl - p) ± 2^x1 - p^Jxl - p] . (3.44) 

The Heaviside function Q{y) of eq. (p.40|) imposes the condition 

y+>0, (3.45) 

which forces the allowed regions for xi and X2 integration to be "region I" and "region 
II" of Fig. p.2| . In these regions, the relation 

y<{i- Vp? (3-46) 
is always satisfied. In addition, we have 

y- > if Xi,X2 G region II . (3.47) 
Following this procedure, we succeed in dividing the phase space into two regions: 



3.1. Kinematics and four-body pkase space witk two massless particles 



35 



Vp 













region I 










region II 







X, , 1 



Figure 3.2: The two different areas in the xi-X2 plane correspond to the 



region I and to the region II of eq. (3.5i). 



1. region I: y can reach zero, where colhnear and soft divergences arise. The 
integration range is (see Fig. p.2|) 



1 r-x2+ ry+ 

dxi I dx2 I dy , 

^ Jx2- Jo 



(3.48) 



where 



X2± 



A{l-xi)+p 



2 (1 - xi) (2 - xi) + p (2 - xi) ± 2 (1 - xi) ^Jxl-p 

(3.49) 



2. region II: y cannot reach and the region is free from infrared divergences. 
The integration range is 



xi+ j-x2- ry+ 

dxi / dx2 / dy , 

Vp J y/p Jy- 



where 



xi+ = 2 



2-p 



(3.50) 



(3.51) 
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We can then summarize the full four-body phase space 

.2 /4^\^ 



(47r)2r(l-e) U'j ^ 

{/•I r^2+ ry+ rxi+ rx2- ry+ ^ 

/ dxi / dx2 / dy + dxi / dx2 / dy>y 
J %/p Jx2- Jo Jy/p Jy- J 

" V ' V ' 

region I region II 

X U [xj - p) [xl - p) - [{xl - %) - (xj - p) - (xj - p) 



X f'dv[vii-v)r^ 

Jo ly^ Jo 



smi 



-2e 



(3.52) 



A statistical factor 1/2! must be supplied if we consider the final state with the two 
identical gluons. 

Sometimes we will need an analogous set of final-state variables, in which the 
role of p and p' is interchanged. The variable y remains the same, xi and X2 are 
exchanged, and the other two variables, denoted by v' and 0', are related to v and 
by the equations 



^a-cos0'^ 



2 ' '2 

1 — cos a — 2 (v — v' cos a) 



1 — cos a — [2Jv{l — v) cos sin a — 2 1> cos a 



(3.53) 



cos ( 



2 sin a \ v'(l — v' 



These relations can be obtained considering the definition of the involved angles 



p' ■ k = |p'||k| (sin cos sin a + cos ^ cos a) = |p'||k| cos^' 
p ■ k = |p||k| cos^^ = |p||k| (sin a sin ^' cos 0' + cos 6'' cos a) 



(3.54) 



where the same scalar products have been computed in the reference system of 
eqs. ( p.2|) -( p75|) and in the reference system where p and p' have been interchanged. 

Exchanging the roles of I and k brings about the following transformations 



V ^1 — V , 



•7r + 



l-v' 



■ TT -t- 0' . 



(3.55) 
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Figure 3.3: Three-body kinematics. 



3.2 Kinematics and tkree-body pkase space witk one massless 
particle 



We consider the three-body process depicted in Fig. |3^ 



e+ (p'J + e- (pe) ^ Z/^ (q) ^ Q{p) + Q{p') + g{k) 



(3.56) 



where the momenta of the particles satisfy 



2/2 2 

p = p = m , 



(3.57) 



For unoriented shape variables, we can express the three-body phase space in terms of 
two variables, that we choose to be Xi and X2 of eq. ( p.6|) . Using eqs. (|3.36| ) and (|3.371 ) 
we have 



H I dxi 



dx-) 



X U (xl - p) (xl - p) - [x] - (xl -p)- {xl - p) 



(3.58) 



where we put y = 0, because the mass of the light system is zero (final gluon on- 
shell). The integration range can be found by imposing the reality condition (|3.39|) , 
with ?/ = 0. 
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Figure 3.4: Four massive body kinematics. 



3.3 Kinematics an J four-body pkase space witk four massive 
particles 



The four-body process describing the production of four massive quarks is illustrated 
in Fig. |3.4| and summarized by 



e+ (p'J + e (pe) ^ Z/^ (q) ^ Q{p) + Q{p') + Q{r) + Q{r') 



(3.59) 



where 



2/2 2/2 2 

r = r = p = p = m 



The four-body phase space is obtained with a procedure similar to the one given 
in Sec. with the simplification that now the entire cross section has no soft or 
coUinear divergences, so that we can put ourselves directly in d = 4 dimensions and 
we do not need to divide the allowed phase-space into two different regions. In the 
centre-of-mass frame of one heavy quark-antiquark pair we have 

r = (ro, |r| sin6'sin0, |r| sin6'cos0, |r| cos^) 
r' = (ro, — |r| sin6'sin0, — |r| sin6'cos0, — |r| cos^) 



P 



P„ I 1,0.0.^1 



p' = Pi 1.0. 



1 2 a, 

P'o 



\ 



1 2 cos a 

p'o 



(3.60) 



where 



(r -|- r 
q2 



/\2 
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and po, p'q and cos a are given by ( p.9|) and (|3.11| ), while 



Tq — m? 



(3.61) 



Starting from eq. (|3.15|) , we need to re-compute the two-body phase space in c? = 4 
dimensions for massive particles 



2ro(27r)3 2r^(27r)3 
|rp(i|r| sin^fi^ 



{27:)^ 6\t - r - r') 



6 (to - 2ro 



1 1 



4rg (27r)2 

1 r2TT 



5 (to - 2ro) 



1 - - / dv r d(t), (3.62) 
y Jo Jo 



4r2(27r)2 " 4(27r)2' 

where we have used eqs. ( [3. 611 ) and ( |3.23| ). The four-body phase space, according to 
eqs. ( p.l5| ), ( p.36| ) and ( |3.62| ) becomes 

^4 



(i$4 = ^ f, / dy dxidx2jl — - dv 

Utt) J V y Jo Jo 



The physical boundaries for the kinematical variables are 

P<y<il-y/pf, y/p < Xi,X2 <l-2p 

and from the reality condition ( p.39|) , we have 

X2- < X2 < X2+ , 

where 

X2± 



(3.63) 

(3.64) 
(3.65) 



1 



4(1 -a;i) +p 



d^A 



(2 -x,)i2 + p-2y- 2xi) ± 2^{xl - p) [{xi -l + y^-py] 

(3.66) 

(3.67) 

(3.68) 

(3.69) 



According to eq. ( p.64| ), X2+ must satisfy 
which implies 

xi<l-y- y/py , 

so we have 



y+ 



dyjl 



p 



y J VP 



X2+ 



dxi / dx2 / dv 



X2- 



1 r2TT 



Jo 



where 



(3.70) 



xi+= i-y - y/py ■ 

A statistical factor l/(2!2!) = 1/4 must be supplied to eq. ( p.69|) , because of the 
presence of two pairs of identical particles in the final state. 
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Next-to-leading-order amplitudes 



4.1 Introduction and notation 



In writing the amplitude for the process we are investigating, we disregard, at present, 
the contributions coming from the decay of the Zj 7 boson into a hght couple of quarks. 
We postpone all the comments on this subject to Sees. ^]5| and |4.6.2| . We can then 
write the amplitude, up to an irrelevant phase, as 



A 



(4.1) 



where |/) refers to states with four- momentum q and Jy(0), J4(0) are the currents 
that describe the decay of the vector boson into the final state |/). In this equation 
we see the contributions coming from the 7 propagator, second term, and the Z 
propagator, first term, where we have neglected terms proportional to q^q^ because 
of leptonic current conservation, always verified for the vectorial part but verified in 
the massless limit for the axial part 



The notation used in eq. ( 



with 
II) is 

9z = 



J- V 



7^ r:^ 



9 



7^7^ and rrie 







(4.2) 



Vi 



2 sin 6*-^ cos 6'w 
= Tq,- — 2c,- sin^ 9 



'-3i 
T3i 



w 



(4.3) 
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where g is the electromagnetic couphng, T^i is the third component of the (left) 
isospin of fermion Cj is its electric charge in units of the positron charge and is 
the Weinberg angle. and are the Z mass and total decay width. 

We are interested in describing only unoriented events: for this reason, we can 
neglect the axial- vector interference term in the square of the amplitude. In fact, for 



the three-parton final state of Fig. p.3|, there are not enough momenta to construct 



an invariant with an e symbol. For the four-parton final state of Fig. |3.1| , one could 
in principle build such an invariant, but the cross section must be symmetric in the 
light parton momenta, so that such an invariant cannot survive. This is strictly true 
for the QQgg final state, and for the state QQqq, if the weak current is coupled to 



the heavy quark. We refer to Sec. [4.6.2| for further details. 

In addition, we have no problems between the use of the dimensional regularization 
procedure and the axial coupling. In fact, we can circumvent the presence of 75 
considering the case of a generic vector current coupled to two fermions with different 
masses mi and 1712. One can easily convince oneself that the case of the axial coupling 
can be obtained by setting mi = m and m2 = — m, since one can turn — m into m by 
a chiral rotation. This procedure is bound to work if there are no anomalies involved 
in the calculation, and this is certainly the case for the graphs we have chosen to 
compute (see Sees. ^^4.(j| ). 

When squaring the amplitude to compute the differential cross section, we have 
two different tensorial structures which describe the leptonic and hadronic part of the 
process. The leptonic tensor, obtained by averaging over the initial polarization, is 
given by 



(4.4) 



and, if we further average over the incoming electron beam direction, we obtain 



3 



9'"' + 



1^, 



(4.5) 



We can then write the differential cross section as 

4:71 



da 



3g2 



P2(g )(w, +aJt;Q + c,CQ-2pi(g )veV^CeC, 



+ dTA p2{q'^) (vl + al) al 



4.2. QQ cross section 



43 



where a = g'^/{47c) is the electromagnetic couphng constant, Nc = 3 is the number of 
colours and 



2n 



1 g2 ^^2 _ ^2 



4 sin^ cos2 (M| - q^y + M| Tl 

A 



P2(g ) = ' 



4 sin^ cos2 / (M^ - q^f + VI 



We have also defined 



n 

-Mi^i = l|^f-^?.. + ^)(0|J(^M(0)|/„)(/„|J{;M(0)|0), (4.6) 

where represents the ri-body phase space, and represents an n-body final 
state. The q^(f term in the projector in eq. (^4.6| ) is, of course, irrelevant for the vector 
current component, but it should be kept for the axial current when the quark mass 
is non-zero. 

In the following we will be interested in strong corrections up to the second order, 
and in the final states: QQ^ QQ9-, QQqQi QQqq and QQQQ. We will use the following 
simplified notation: 

(2) — 

- J^v/A for the QQ tree-level term 

- J^yj^ to indicate the three-body Born QQg, O (a^) term 

- -MyJa to indicate the three-body virtual QQg, O {a1) term 

- Myj\ or Mgg foT thc four-body QQgg, O {al) term 

- M^^\ or Mqq for the four-body QQqq, O (a^) term, 

and equivalent ones for the dTy/A terms. 

We will drop the V/A suffix when not referring specifically to the axial or vector 
contribution. 
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"(Tifirff-tP 



Figure 4.1: Tree (T) and Born (B) Feynman diagrams. 



4.2 cross section 



In the higher part of Fig. [4.1| we have depicted the Feynman diagram representing 
the tree-level amplitude 

^^/^ = w(p)r^/^t;(p'), (4.7) 

that, once squared, gives, for the two-body contribution to eq. ( |4.6|) , at zeroth order 
in Us, 



where p is defined in eq. (|3.12| ) and 



P 



(4.9) 



Multiplying eq. ([4.8|) by the 2-body phase space j3/{8n), we get the zeroth-order 
contributions to the total cross section 



(2) 



T 



(2) 



(4.10) 



This justifies our choice for the normalization factor in eq. (|4.6|) : in the massless limit. 



rp{2) _ rp(2) 

— 



1. 



4.3. QQg cross section at order as 45 



4.3 QQg cross section at order q;< 



The amplitude for the Born term represented in the lower part of Fig. is 



u{p) 



{p — qY — w? 



v{p') , (4.11) 



where t"^ are the generators of SU(Ai"c) gauge symmetry, i and j refer to the colour 
indexes of quarks, while a to the colour index of the external gluon. We remind here 
that 



a=l k=l 



where 



We define 



N?-l 4 



2Nr. 



z,j = 1...N, 



for Nr = 3 . 



(4.12) 
(4.13) 



9s 

where the sum refers to the spin and colour of the quarks and to the colour of the 
gluon in the final state. In the Feynman gauge, the sum over the polarization of the 
final gluon gives 

Y^e''{k)t^\k) = -g^^' (4.15) 

pol 

so that we can define 

Mv/^ = -wM^7;. (4.16) 

We need the expressions for My/ a in = 4 — 2e dimensions. Computing the trace in 
eq. (CT) , 



Ml 



V 



X2 I ^2 
1 2 



+ 



16 



m 



(1 - xi)(l - X2) (1 - a;i)2(l - \ g2 
- 3 (x? + xl) - 8 (1 - xi)(l - X2) + 2 
16e 



2XiX2(Xi + X2) 

32 / 



1 — Xi)^(l — X2Y \ J ^ 



(1 - Xi)(l - X2) 



Xi + X2 + {I - Xi){l ~ X2) + Xg - 1 



m 



x„ 



(1 - a;i)(l - X2) 



+ 



x„ 



{l-X,){l-X2) ^' 



(4.17) 
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and 



Ma 



X2 I ^2 
1 2 



+ 



16 



m 



[1 - - X2) (1 - - X2)2 V g2 ^ 

- llxiX2(xi + X2) + 8(x^ + xl) + x'^(x2 — 1) + xl{xi — 1) 



12(xi + X2 — 2x1X2] 



+ 2xixi + 4: (1-e) + 



64 



m 



:i-xi)2(l-x2)2 \ g2 



+ 



(1 - Xi)(l - X2) 



X 



3 ' 



(4.18) 



where Xi, X2 and are defined by ( |3.6| ) and (p.38|) . According to eq. ( [4.6| ) we have 



27r 



(4.19) 



where is the mass parameter of dimensional regularization, introduced in order to 
keep Qs dimensionless. 

For a future use, we introduce now a unit three- vector j, belonging to the event 
plane (i.e. the plane defined by p, p' and k) and perpendicular to k. In the centre- 
of-mass system of the process, we can then write 



f = 1 



k-j = 



j = ap + 6k . 



(4.20) 



The generalization of these equations is easily obtained. From the fact that j is a 
unit, purely space-like vector, we have 



■2 -2 -2 

3 = Jo - J 



Jo = . 



Since q= (go, 0) 



and 



9 ■ J = goJo - q ■ j = 
k- 3 = ^0 Jo - k ■ j = 



(4.21) 



(4.22) 



j = ap + bk + cq . (4.23) 
We can then determine the coefficients a, b and c by imposing the validity of eqs. (|4.22| ) 

q ■ j = =^ a {p ■ q) + b {k ■ q) + cq'^ = 

k-j = =^ a{p-k) + c{q-k) = , (4.24) 
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that, once solved, give 



p ■ k 



N 



q ■ k 



b = N 



p-q 



q ■ k p ■ k 



q ■ k 



(4.25) 



with determined by the normahzation condition ( |4.21| ). From the expression of j 
we can compute 



1 



P-J 



2 X n 



xj (4xi - p- A) +xl (4x2 - p - 4) - 2xiX2(p + 6) 



+ 4(p + 2) (xi +X2) -4(p + 1) 



(4.26) 



We are now in a position to give a decomposition of the tensor M°"'^' of eq. ( |4.14] ) 
(we disregard from now on the suffix V/A) according to the direction of j: in fact, 
the Lorentz indexes of this tensor take origin from the composition of the three 
independent vectors q, p and k. Using eq. ( |4.23| ), we substitute j instead of p so that 
we can write M°"'^' in a quite general form 



M"" = M^gY + M^ff + terms involving g or A; 
where (see eq. (p.25|) ) 



9T = 9"" 



1] ■ k 



with 7] = (/co, — k), so that 



r]^ = 



r]-j = 



q.gT = 0. 



(4.27) 



(4.28) 



(4.29) 



Contracting eq. ( [4. 271 ) with g±a(7' and with jo- jo-' , we obtain, m d = A dimensions 

M^=M; + 2M'^'^j.Jo'. (4.30) 

This expression can easily be computed, starting from the consideration of the pre- 
vious paragraph about the dependence of M°"'^' from g, p, k, and from eqs. (HT 
and (CT) 



2 c 



V/A 



4xiX2(Xi + X2) - p{xi + X2f - 4(xi + x\ 



:i-Xi)2(l-X2)2 

- 12x1X2 + 4(p + 2)(xi + X2) - 4(p + 1 



(4.31) 
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where 

Cv = p + 2, c^ = _2(p-l). (4.32) 

We also define, consistently with our previous notation 

•^?M ' = ^-^C, gl^^^ M(^;; , M^^^l^ = ^ gl^^^ M^/i . (4.33) 

In the cases when the V j A suffix needs not be specified, we will simply write M.l'^'^ 
Mi and 



4.4 Virtual contributions 

Corrections to the three-jet decay rate to order come from the interference of the 
one-loop graphs with the tree-level Born ones. In Fig. [4.2| we have depicted these 
contributions. 

The algebra to compute these terms in c? = 4 — 2e dimensions has been carried out 
in a straightforward way, using a MACSYMA program, which reduces the original 
Feynman graphs to a linear combination of scalar, one-loop integrals. The scalar 
integrals have been computed analytically and their values are listed in Appendix |B[ 

Loop corrections to on-shell external lines (not illustrated in Fig. |4.2| ) require 
particular attention. 

We start from the self-energy corrections to heavy-fiavour external lines. As de- 
scribed in Appendix |C.1| , the effect of the fermion self-energy correction to an external 



line, including the mass counterterm, is equivalent to multiply the external propa- 
gator by the factor Zq of eq. ( |C.13| ), so that we have a contribution to M.^ equal 
to 

/3 



z^^M^ = -N^Cpgl I^^J (^^ + 4j X . (4.34) 
We have to consider also the diagrams obtained with the mass counterterm in- 



sertion in internal fermion lines. In fact, according to eq. (|C.11|) , we have to add 
a counterterm in order to keep m as the pole mass, after radiative corrections. 
These diagrams are depicted in Fig. ^]3|. 

Similar considerations apply to the self-energy corrections to external gluon lines. 
As remarked in Appendix |U.2| , gluon, ghost and light-fermion self-energy corrections 
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to external gluon lines vanish in dimensional regularization. Only the correction 
coming from a heavy-flavour loop needs be considered, and it gives a contribution at 
order equal to (see eq. (|C.22|) ) 



(4.35) 



After that, charge renormalization is all that is needed, since we are computing 
a physical cross section. Charge renormalization in the mixed scheme of Ref. 



IS 



described in Appendix |C.3[ From eq. (|C.24|) we can see that the last correction to 
our virtual term is equal to 



Tpuii - —Ca +1 — 1 —Tf 



3e 



3e 



3e 



X Mb . 



(4.36) 



We can now summarize the combined effect of external line corrections and renor- 
malization to be included with Ai^ 



m? I \e J \3e 3e / \ 



4 



11 



X Mb ■ (4.37) 



The factor of 2 in front of the fermion external line corrections is to account for the 
two fermion lines. 
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Figure 4.2: Radiative corrections to the process Zj^-^QQg. The interfer- 
ence of these diagrams (left hand side of the figure) with the tree-level Born 
term ( right hand side ) gives rise to contributions to order a1 . 



4.5. Ha giwara contritutions 
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Figure 4.3: Contributions to order coming from the interference between 
diagrams with the mass counterterm insertion and the tree-level Born term. 



4.5 Hag iwara contributions 

At order a^, we have also contributions coming from the interference between terms 
in which the weak current is coupled to the heavy quarks and to quarks of different 
flavours. 

These diagrams are represented in Fig. the higher part of the figure describes 
the contributions coming from light-quark loops, while the lower part describes the 
contributions coming from heavy-quark loops. 

By C-invariance (Furry's theorem), these diagrams vanish for vector currents. For 
axial currents, they cancel in pairs of up-type and down-type quarks, because they 
have opposite axial couphng (see eqs. ( |4.3| )), as long as the loop of different flavour 
may be regarded as massless. Thus, the up-quark contribution cancels with the down- 
quark, and, if the charm mass is neglected, the charm contribution cancels with the 
strange. Only the graph with a top quark loop remains. 

Paired to this last type of diagrams are the contributions where the massive quark 
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Figure 4.4: Radiative corrections to the process Z/j ^ QQg (Hagiwara con- 
tributions). The interference of these diagrams (left hand side of the figure) 
with the tree-level Born term (right hand side) gives rise to contributions to 
order . 



in the loop is the same as the heavy quark in the final state. 

We have not included these diagrams in our calculation, because they have been 



computed in Ref. |T3], where it was shown that their contribution is of order of 1%. 



4.6 Real contributions 

The square of the Feynman diagrams with four particles in the final state gives rise 
to the real contributions. These amplitudes are easily obtained, in an analytic form, 
with a little "Diracology" in c? = 4 dimensions. The behavior of these amplitudes 
in the soft and collinear limit is obtained in the next chapter, where the asymptotic 
forms of the amplitudes in derived in ci = 4 — 2e dimensions. 
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4.6.1 Real contributions to tke QQqQ cross section 

The diagrams contributing to the process 

Z/-f-^QQgg 



e e 



(4.38) 



are depicted in Fig. [4.5| . From the square of these eight diagrams, we can obtain thirty- 
six terms, but most of them are related by interchange of the momentum labels, so 
that only thirteen amplitudes need be considered. We follow the notation used in 
Ref. and indicate with Bij the interference of diagram Bi with diagram Bj (i> j). 

The different contributions to the cross section can be classified into three classes, 
according to their colour and spatial structure. We will always factorize out the colour 
factor common to the Born term, equal to CpNc, so that we have: 

1. Cp class: planar QED-type diagrams 

2. Cp — \Ca class: non-planar QED-type graphs 

3. Ca class: QCD graphs, involving the three-gluon vertex 



where Ca = Nc. In Tab. |4.1| we collect the thirty-six terms and the label interchanges, 
needed to compute all of them from the thirteen we have calculated (first row). 



and that are represented in Fig. |4.6| . In the first row of this figure, there are the 
contributions belonging to the first class, in the second, the contributions to the 
second class and in the last row, the contributions to the third class. One final 



permutation 


Cp class 


Cp — \Ca class 


Ca class 


(1^2) 
(3-4) 
(1 ^2) (3^4) 


Bll B21 B22 B32 

B64 B66 
B44 B54 B55 B65 

B31 B33 


B41 B42 B53 B52 
B61 

B51 B62 
B43 B63 


B71 B72 B82 B77 B87 
B84 B86 B76 B88 
B74 B75 B85 
B81 B83 B73 



Table 4.1: Label interchanges needed to compute the thirty-six terms contributing to 
the process e"^e~ ^ Z/7 ^ QQqQ- 



remark is needed, in order to sum over the final gluon polarizations. In fact, due to 
non-conservation of gluon current, we can use eq. ( [4.15| ) to sum over polarization only 
if we include B7- and B8-like diagrams with "external" ghost. 
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Figure 4.6: The thirteen amplitudes necessary to complete the computation 
of the cross section for the process Z/^ ^ QQgg. Cut propagators refer to 
on-shell particles, while numbers label the different momenta. 



4.6.2 Real contributions to tke QQqQ cross section 

The diagrams contributing to the process 

e+e- ^Z/^^QQqq (4.39) 
are illustrated in Fig. [4.7|. From the square of these four diagrams, we generate ten 



terms. In Fig. ^]8| we have depicted nine of them, because the tenth, A41, can be 
obtained with the massive loop of diagram A31 and the massless loop of diagram 
A42. 



56 



Ckapter 4. Next-to-leaJing-orJer amplitudes 




Although the colour coefficient is the same for all the ten diagrams (Tp), their 



infrared structure is different: in fact, only the diagrams in the ffist row of Fig. ^ 



contain infrared divergences, due to the massless- quark coUinear region, all the other 
diagrams being finite. 

For this reason, it is mandatory to include the ffist three diagrams, to check 
infrared cancellation, and it is custom to assign the other diagrams to the light chan- 
nel: in fact, they are, in general, characterized by a large invariant mass for the light 
quarks and a small invariant mass for the heavy-quark couple. We have not included 
these last diagrams in our calculation. 



4.6.3 Real contributions to tke QQQQ cross section 

The eight diagrams contributing to the process 

e+e-^Z/^^QQQQ (4.40) 

are illustrated in Fig. [4.9| . From the square of these diagrams we have thirty-six 
infrared-safe amplitudes. We need to compute only twelve of them, that we have 
represented in Fig. [4.1U| , because the others can be obtained by permuting the quark 
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All 



A21 



A22 






A33 



A43 



A44 





A31 



A32 



A42 



Figure 4.8: The nine amplitudes necessary to complete the computation 
of the cross section for the process Z/j ^ QQqq. Cut propagators refer to 
on-shell particles. 



momenta. In Tab. ^.2| we give the twelve amplitudes we have computed (first row) 
and the label interchanges necessary to build all the others. As previously done for 
the QQgg process, we divide these contributions according to their colour factors and 
according to their structure: 



1. 1®* Tp class: diagrams illustrated in the first row of Fig. |4.1U| 

2. Cp — \Ca class: diagrams from D81 to D62 

3. 2"*^ Tp class: diagrams composed by two different massive loops, coupled to 
the weak current, that can be split apart by cutting the two joining gluon lines 
(singlet contributions). 



58 



Ckapter 4. Next-to-leaJing-or Jer amplituJes 



permutation 


1^* Tp class 


Cp — \Ca class 


2"^^ Tp class 


(1^3) 
(2-4) 
(1 ^ 3) (2 ^ 4) 


DU D21 D22 
D55 D65 D66 
D77 D87 D88 
D33 D43 D44 


D81 D52 D51 D82 D71 D62 
D54 D61 D64 D53 
D72 D83 D73 D84 
D63 D74 


D31 D32 D42 
D75 D76 D86 

D85 

D41 



Table 4.2: Label interchanges needed to compute the thirty-six terms contributing to 
the process e+e" Z/^^ QQQQ. 




Figure 4.9: Feynman diagrams contributing to the process Z/j ^ QQQQ. 
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Figure 4.10: The twelve amplitudes necessary to complete the computation 
of the cross section for the process Z/^^QQQQ. Cut propagators refer to 
on-shell particles, while numbers label the different momenta. 
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Infrared cancellation 



5.1 Tke subtraction metkocl 

In this section we want to introduce the method we have used in order to deal with 
the infrared divergences. In fact, as we have seen from the previous sections, infrared 
divergences arise both in the virtual and in the real terms: 



1. in the virtual terms of Sec. divergences arise from loop integration, as you 
can see from Appendix Since we have chosen to regularize our integrals with 
a dimensional method, we obtain poles in e; 



2. in the real terms of Sees. 4.6.1 and 4.6.2| , divergences in the differential cross 



section appear when we integrate over a particular region of the phase space. 
In the diagrams depicted in these sections, we have propagators proportional to 
1/y (according to eq. (p.6|), y is proportional to the mass of the light system), 
so that we have divergences when integrating over the phase space region of 
small y. In this region, a four-particle final state masquerades as a three-particle 
event, because one of the emitted particles becomes soft, or two particles become 
coUinear. Notice that the Born and the virtual graphs always have y = 



According to B loch- Nor dsiek ||26| and Kinoshita-Lee-Nauenberg theorems, in 



frared divergences must cancel for sufficiently inclusive physical quantities. 

Now, we show how we have checked this cancellation. If we compute the O (a^) 
differential cross section, in d = 4 — 2e dimensions, we get a result of the form 

da = da^'^ + ( ^\ da^'^ + ( da^'^ , (5.1) 

\27T/ \27l/ 
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with 

= + ^^$3 + (5.2) 

(i$2 (i$4 

where we sum up the contributions coming from two-, three- and four-particle final 
state. The two-particle final-state contribution at this order arises from the inter- 
ference of two-loop diagrams with the tree-level term. We have not computed this 
differential cross section term, because we want to calculate only three-jet-related 
quantities, that give zero contribution in the two-jet limit. We give a detailed de- 
scription of this fact in Sec. pl2| , so that we neglect in eq. ( ^.2|) . 

Using eqs. (|3.58|) and (|3.4CI|) , we can rewrite the three- and four-body phase space 



as 

d^3 = dXidx2 J3{xi,X2) 

(i$4 = dxidx2dyd9d(f)J4^{xi,X2,y,9,(j)), (5.3) 

where J3 and J4 represent all the other factors, whose exact structure we do not need 
to know in this section. 

In order to implement the cancellation of the soft and coUinear singularities, we 
now imagine to compute some physical quantity G {G stands for "Generic"), which 
depends on the final-state variables. G may be a combination of theta functions that 
characterize a histogram bin for some infrared-safe shape variable, like the thrust, 
the c parameter and the heavy-jet mass. In general the definition of G is specified 
for any number of particles in the final state. Since we are dealing with three- and 
four-parton final states, G is characterized by only two functions, G3{xi,X2) and 
G^i^Xi, X2, y, 0, (p). Soft and coUinear finiteness of G requires that 

lim G4(a;i,a;2,?/,6',(/)) = G'3(xi,a;2) . (5.4) 

When computing the contributions to G at second order in a^, we have 

/C dcT^'^^ 
da^'^^G = J dxidx2J3{xi,X2)-j^G3{xi,X2) 

f da? 
+ dxi dx2dy d6 dcj) J4{xi, X2,y,0, (p) ——— G4{xi, X2,y,0, (f)) , (5.5) 
J d^A 



where each term on the right-hand side contains soft and collinear divergences that 
cancel in the sum. The integration of the complete differential cross section is too 



5.1. Tke subtraction metkocl 



63 



difficult to be performed analytically, but it is not a limiting point. In fact, we can 
add and subtract a suitable quantity, to obtain 

.(2) . ^^(2) 



J da^'^^G = J dxidx2G3{xi,X2)<^-^^J3{xi,X2) + Jdyd6d(j) J4{xi,X2,y,6,(, 
+ j dxidx2dyded(j)J4{xi,X2,y,6,(f))<-^G4{xi,X2,y,6,(f))--j^G3{xi,X2) 



(5.6) 

where al is chosen in such a way that it has the same soft and collinear singular 

(2) 

behaviour of , or, mathematically 



da^ 



(2) 



jim ^ = 1 . (5.7) 



dn 



d^A 



(2) 

The aim of this procedure is to obtain an approximation for dai that can be in- 
tegrated in dyd9d(f), so that the divergent parts appear as poles in e. In this way, 
the ffist term of eq. ( |5.6| ) can be computed analytically: the single and double poles 
present in da^ / d^^ all cancel with the poles arising from the dy d9 dcp integration of 
dal /d^4^, and thus this term is finite. 



The second term in eq. (|5.6| ), because of eqs. (|5.4|) and ( ^.7|) , has no soft or collinear 
singularities, and thus can be evaluated directly in four dimensions^. 

To implement numerically the computation of this term, we first generate a four- 
body configuration xi,X2,y,9,(j), in the corresponding four-body phase space, with 
weight J4. We associate to this configuration two events: one four-body event, with 
kinematics Xi,X2, y,0,(f) and weight da'^^^/d^^, and one three-body event, with kine- 

(2^ 

matics Xi, X2, y = 0, and weight —da^ /d^4. The computation of a shape variable 
using the above scheme reproduces exactly the second term of eq. ( |5.6| ). 



^ Observe that both eq. (5^) and eq. (5/7) must be satisfied in d dimensions in order for this 
argument to apply. 
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5.2 Two-loop diagrams 

We want now to show why, in at least two cases, we do not need to compute the 
two-body contribution at order a^, to the differential cross section. 

1. If we compute the average value of some physical quantity, we have 

G ^ ^ = i 1,0,...) + (^) / x„ 0, . . .) 

+ (^) ' / da^''>G{xr, X2, y, ^, 0) + O (a^) | (5.8) 
where a is obtained by integrating da of eq. (|5.1|) 



and where we no longer label G according to the number of particles in the final 
state, which is completely specified by the arguments of G itself: 

- for the two-particle final state, xi = X2 = 1 and y = 0, while 6 and are 
undefined, so that G2 = G{1, 1, 0, . . .) 

- for the three-particle final state, we must give Xi and X2, but y = 0, while 
9 and are undefined, so that G3{xi,X2) = G{xi, X2, 0, . . .) 

- for the four-particle final state, we have to specify all the five arguments, 
so that G^ixi, X2, y, 0, 0) = G{xi,X2, y, 0, 0). 

Adding and subtracting the two-body contributions in the integrals, we obtain 
G = ^|aG(l,l,0,...)+(^)y'rfaW[G'(xi,a;2,0,...)-G(l,l,0,...)] 

+ (^)V [ciorf + c^af + daf ] [Gix,,X2,y,e,<P)-Gil,l,0,...)] 
+ (5.10) 
where we have used eq. ( ^.2| ) to expand da^'^\ Since 

da2 = 6 {xi — 1) 6 {x2 — 1) 5 (y) (T2 dxi dx2 dy , (5-11) 

(2) 

the contribution of da2 to the second integral is zero. The same thing happens 
for the two-body contribution, at order ag, da2^^: in fact, this term too is 
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proportional to the product of the three (5's, and gives zero contribution to the 
first integraL 

The dependence from a2^^ and a2'^ is in the expression of a. We can then expand 
the denominator till order a^, to obtain 



G = G(l,l,0,...)+(^)|^[G(a;i,a;2,0,...)-G(l,l,0,...)] 

[G(xi,X2,2/,^,0)-G(l,l,O,...)] 



a 



^(0) / iw °' • • •) - 1' 0' • • •)] I + («^) . (5.12) 

(2) 

We can see that no dependence from a2 is left, confirming that we do not 
need to compute two-loop diagrams. We need, however, the total two-body 
contribution at first order a^2^ , but the analytic expression of this term is known 



from long time (see, for example, Ref. |28 



2. The second case, where we do not need to compute the two-loop diagrams, arises 
if we deal with physical quantities that assume zero value in the two-jet region. 
In fact, starting from eq. (|5.10| ), not normalized to the total cross section 

(G) = J daG (5.13) 

we see that if (7(1, 1, 0, . . .) = 0, we accomplish our goal. This is what we have 
done in our program: we have computed quantities like 

- 1 —t, t being the thrust of the process 

- the c parameter 

- M| — m^, M| being the heavy-jet mass 

which give zero contribution in the limit xi—>-l,X2—>-l,y—>- 0. 



5.3 Soft and collinear limit of tke QQgg cross section 



In order to apply the subtraction method described in Sec. ^.1| , we need to derive 



an expression for the singular part of the four-body cross section, valid in both the 
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coUinear and the soft limit. These hmits are both characterized by y ^ 0, except that, 
in the soft hmit, at the same time, v -^0 {I soft) or — > 1 {k soft). In fact, from the 
definition of y given in eq. ( p.6|) , we see that, if k becomes colhnear with / or if A; or 
/ become soft, we have y—>-0. In this limit, from eqs. ( |3.2| )-( pTTlD , we have 

cos a 1 — ay 
sin a ~ ^J2a^ 



where 



Po V Po 



\R1 = 

P'o \ 



1 



Po 



i-fiy 



p'o \ 



1 2 cos a ~ 1 — 6 y 

p'o 



p'o \ 



1 TT sin a 

Po 



2a^/y 



h 

/2 



(1-Xi)(l-X2) 

2 

(l-a;i)(l-a;2) 

2m2 1 

(1 — XiY 

2m2 1 



Xi + X2-I 



m 

X1+X2-I TT 



2 + 



2 + 



1 — X2 1 — Xi 



1 — Xi I — X2 
1 — Xi 



I - X2 



We can then write, always in the y — > limit. 



p ■ I 



p ■ I 



p ■ k 



9^(1 - X2) 
4 

- xi) 



- X2) 



f2{l-2v) 



2v 



1 - a;2 



2v + y 

2v — 2y/2a cos0 y'f (1 — f)^^ + V 

2(1 -t;) 

1 - X2 



6(1 -2t;) 



2v 



1 — Xi 



2(l-t;)-|/ 



/2(l-2^;) + 



(5.14) 
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p' .k^ ^ '^'' '(2(1 - v) + 2\/2a cos^Jv{l-v)y/y 



1 — Si 



p.p — 



-(1 - xi - X2) 5- + 2/ 



where we have made use of eq. ( p.23| ). We can see that 



1 soft 



and 



(p + lf -m^ = 2p-l^0 
ij)' + lf -m^ = 2p -1^0 



y — i> and f — s> 



p ■ I 



1 



m 



2(1-X2) 

1 



1/ + (1 - X2fv 



m 



p' ■ I 



2(1-X2) 

4 



[y + hv] 



2V2a , ^ ^ 2 
y — cos (p vvVy + ^ ^ 

y — c cos \fv\fy + gv 



- k soft 



and 



{p + ky - m= 2p ■ k ^0 
{p' + kf -m^ = 2p' -k^O 



y -^0 and v ^1 



p ■ k ~ 



1 



m 



2(1-X2) 

1 

2(1-X2) 



l/ + ^(l-a;2)'(l-i;) 



[y + h{l- v)] 



p' ■ k ~ q 



2 6(1 - xi) 



2V2a , , ^ 2 , 

y H 7 — cos VI - t'Vy + 7 (1 - f j 



6(1 -Xi) 

4 



?/ + c cos (j) y/l — + g {1 — 
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where 



h= -^{1- X2f 



2V2a 
b 

2 

b ' 



(5.15) 

(5.16) 
(5.17) 



5.3.1 Soft contribution 

We begin with the soft singularities of Aigg (since the same formulae apply irrespective 
of the vector or axial case, we will always drop the V/A suffix). They are given by 
eq. (p.45|) , which we now rewrite 



M 



soft 
99 



p ■ k 



+ 



p' ■ k 



{p-l){k-l) {p'-l){k-l) 



+ 2 C 



p ■ p 



2 J {p-l){p'-l) 



m 



+ 



m 



+ {k ^ l) \ X Mb . 



(5.18) 



_{p-if {p'-iy_ 

Exploiting the formulae ( p.l4| ) for the approximated scalar products in the soft limit, 
we can write, in the limit of / soft, 

p-k 2h 1 



{p ■ l){k ■ I) y[y + h v] 



p ■ p 



K 



1 



{p ■ I) {p' ■ I) m? y + hv y — c cos (py/y^/v + gv 
^2 4h 1 



= Ep^p';l{Xl,X2,y,V,, 



m 



= Ep^p.i{xi,X2,y,v) 



(5.19) 



{p ■ ly [y + hvY 

where the first two indexes in E refer to the numerator, while the third index refers 
to the soft momentum, and 

2m2' 



K 



X2 4 



1 — xi b 



Xi+ X2-1 - 



(5.20) 



We will also need analogous formulae in which the roles of p and p' are interchanged. 
With the help of eqs. ( ^.531 ), we have, in the limit of / soft, 

p' ■ k 



{p'-l){k-l) 



E'p>,k-Axi,X2,y,v') = Ep^k;i{x2,xi,y,v') 
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p ■ p 



{p-l){p'-l) 



^p',p;«(a;i, 2:2, y, v', 0') = Ep,p>.i{x2, xi, v', 0') 



m 



Ep,p,.i{xi,X2,y,v') = Ep^p.i{x2,xi,y,v') . 



ip' ■ ly 

The contributions for the case when k is soft are instead obtained from the above 
using eqs. (|3.55|) . For example 

Ep^i.k{xi,X2, y, v) = Ep^k;iixi, X2, y,l-v) 

X2, y, V, 4>) = Ep^p>.i{xi,X2, y,l - v,4> + Tc) . 
We can now write down our approximate soft cross section: 



(5.21) 



M 



soft 
99 



9s f^^) Ca Ep^k-i + Ep, f..i + Ep^i-k + -E'p/^/.fc 



+ {Cf ^) Ep^pi-i + E'p,p.i + Ep^pi.f, + E'pip.^ 



Ep,p;l + Epipi.i + Ep^p-k + Eplpl.j. 



X Ml, 



(5.22) 



The soft cross section written in this way is symmetric under the interchange of k 
and I, and of p and p' . 



5.3.2 Collinear contribution 



The colhnear part of the cross section that receives contributions from the QQgg final 



state can be written, according to eq. (p.29|) 
4Ca f 



M 



coU+soft 
99 



2 2e 
Qsl^ 2 



q y 



z 1 — z 



-2z{l- z){l -e) 



gaa' 
g±acT' 



X M 



<7(T 
h 5 



(5.23) 



kl 2 - 2e_ 

where z is the momentum fraction of / versus I + k m the collinear limit. It can be 
chosen to be equal to v or to v' . Notice that, as it is explained in Appendix p.l| , 
parts of the collinear singularities are already contained in the soft-limit expression. 
In fact, for y — at f fixed, we have 

2 „ „, 2 



q'^yv 



Ep,i,k ~ E'„, 



p',l;k 



q^y (1 - v) 



(5.24) 
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Thus, the Xj z and 1/(1 — z) terms in the colhnear hmit formula ( 5.23| ) should not be 
included, since they are already present in the soft term, and we get 



M 



coll 
99 



2 2e 
9s o 



ry 



[-2 + v{l-v)]g, 



2v{l~v){l~e) 



9±o 



ki 2-2.. f^'^^' 

According to eq. ( p.ll| ), the perpendicular direction refers to a direction orthogonal 
to / + A; in the centre-of-mass system and in the colhnear limit. For this reason, 
lies in the same perpendicular plane as the vector j defined in eq. ( |4.20 ), and forming 
the angle with j, so that 

(k±-j)^ 2, ^ 



cos 



— COS 



(5.26) 



Using eq. ( |4.27|) , the azimuth-dependent term of eq. ( |5.25D becomes 



M 



coll 
99 



-2v{l - v){l - e) 



Mt + Mi 



-v(l - v) Ml 



{k±-jy 



kl 

2(l-e) + l 
1 



2 - 2e 



(5.27) 



= -v{l- v) Mi [-2(1 - e) cos^ q 

Considering now the integration over the allowed phase space of eq. ( p.52| ), we obtain, 
for the d) contribution, 



smi 



\-2e 



M 



coll 
99 



smi 



-2e 



-2(1 -e) cos^0 + l 



. 



(5.28) 



Despite of this fact, this term must be present in the program, because, for e = 0, the 
corresponding phase space integral becomes indefinite, since the y integration is di- 
vergent, and the azimuthal term (that gives zero) preserves the numerical integration 
to fail. 

We thus arrive to the following expression for the colhnear term to be added to 
the soft term 



M 



coll 
99 



2 2e4C^A 

9sy^ ^— 

q y 



M, 



v{l-v)+v'{l-v') 



+ 



Mi 



v(l-v) 2(1 -e) cos' 



1) +v'(l-v') (2(1 -e) cos 



-,2 A' 



(5.29) 



where we have symmetrized the expression in v, v' and 
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5.4 Collinear limit of tke QQqq cross section 



With a procedure analogous to that one used in Sec. |5.3.^ for the collinear part of the 
QQgg cross section, we can obtain the colhnear part of Aiqq. From eq. ( p.35| ), we 
have 



M 



coll 



2 2. 4nifTp J 



1 



Mi 



4(1-6) 

vil — v)[2 cos^ & 



1-e 



1-e 



. (5.30) 



The expressions A^^f (eq. (p:22|)), A^^f (eq. dOOQ ) and A^^f (eq. (|O0| )) depend 
upon xi and X2 via M.^ and Al^. These expressions are meaningful only if xi and X2 
belong to the domain of the three-body phase space. We thus define 



Al,, = (Ai^f + Ai^f ) Q,{x,,x,) 
Mgg = M':f e3(xi,x2), 



(5.31) 



where the G3 function is precisely defined to be zero when xi and X2 are outside 
the three-body phase-space region. More specifically, using the integration limits of 
eq. (|3.58D, we have 



QsiXl, X2) = 6(1 - Xi) Q{Xi - ^/p) Q{X2+ - X2) 0(X2 - X2-] 



(5.32) 



We are now in a position to specify the subtraction procedure outlined in Sec. p?T| . 
Our expression for the second-order contribution to a soft- and collinear-safe quantity 
G is given by 



- j Mgg{xi, X2, y, V, 0) ^(xi, X2, y, v, 0) 

+ J d^4Mqq{Xi,X2,y,V,(f))G{Xi,X2,y,V,(f)) + J d^3My{Xi,X2)G{Xi,X2) , 

where all quantities are computed in d = 4 — 2e dimensions. The factor 1/2 in front 
of the gg contribution accounts for the two identical gluons in the final state. We 
rewrite the above expression adding and subtracting the same expression, that is the 
soft and collinear limit of the cross sections, to obtain 

- j d(^4 (^Mgg{xi, X2, y, V, 0)G(xi, X2, y, v, 0) - Mgg{xi,X2, y, v, 0)G(xi, X2)) 
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+ j (i$4 (^Mqq{Xi,X2,y,V,(f))G{Xi,X2,y,V,(j)) - Mqq{Xi,X2,y,V,(f))G{Xi,X2)) 

+ J d^3 (^M^{xi,X2) +Mi{xi,X2)) G{xi,X2) , (5.33) 



where we have defined 



Mi{xi, 2^2) = 2 y ^*4/3 Mgg{xi,X2, y, f , 0) + J d^4/3 Mqq{xi,X2, y, V, 0) , (5.34) 
and (i$4/3 is defined by 

d$4 03(a;i, X2) = rf$4/3 . (5.35) 

The exphcit expression for (^$4/3 can be obtained from eqs. ( p.52|) and (|3.58|) . We 
first notice that the four-body phase space is almost proportional to the three-body 
phase space, except for the ratio 



4 {xl - p) {xl - p) - \(xl - Ay) - {xl - p)- {xl - p) 



\ 4 (a;? - p) {xl - p)- [xl - {xl - p) - {xl - p) 



l + 0{ye) 



(5.36) 

On the other hand, terms of order ye can be neglected, since they cannot generate 
infrared singularities, because of the y factor, and therefore they can only produce 
terms of order e. Thus we can write 



ty+ /"i 1 

d^4/3 = N,R.q'q-'' dyy-' dv [v{l - v)]-' — , 

JO Jo i\<h Jo 



sm( 



-2e 



(5.37) 



or the analogous one in the v', 0' variables. The normalization factor A^^^ is defined as 

1 



while 



-iN{e) 
1 



167r2 



(47r)^r(l + e) 



(5.38) 



(5.39) 



r(l + e)r(l-e) 6 
Since we are free to choose the set of variables we prefer in the (i$4/3 integration, it 
is easy to see that the J^i{xi,X2) term reduces to 

2 q^y 



Mi{xi,X2) = gip^'Jd^^/3 
2 nif Tp 1 



+ 



q^y 1 - e 



v^ + {l-vf-e 



+ 2 



_Ca\ 



2 / 



Ep^p'-i — 4 CpEpp. 



X Ml 
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where the term proportional to has been dropped, since it vanishes in c? = 4 — 2e 
dimensions, after the azimuthal integration, as shown in eq. ( p.2S ). 

We can now integrate over (i$4/3 to obtain the analytic expression of A4i{xi, X2). 
We define, for the coUinear term. 



rcoU 
99 



dyy " dv [v{l — v)] 



[1 - elog (?/+)] 



11 67 



sm( 



-2e 



1 



[v{l-v)-2] 



6 



e +0 e) 



(5.40) 



and 



rcoU 



y+ 



dyy " dv [v{l — v) 



N, 



smi 



-2e 



1 v^ + {l-vf 
y 1 - e 



1 



10 



-[l-elog(y+)] ^- + -ej+0{e) . 
For the integration of the soft term, we define 



^p,k;l — Q 



dyy " dv [v{l — v)]" 



N^Jo 



[smi 



-2t 



p,k;l ; 



(5.41) 



(5.42) 



and the analogous ones for Ip^p-i and Ip,p'-i- With this notation, we have 

„2 



Ip,p;l = 4/l h 



Ip,p';l — K Is 



(5.43) 



where the values of Ji, I2 and Is are collected in Appendix 0, and are given precisely 
by eqs. ( ^121) , (gJ^ and ( W^ - 

Our final expression for Aii{xi, X2) is therefore 

,2' 



Mi{xi,X2) = N,R,gU^ 



coll 



+ 



2 CAlp,k;l + 2 ( Q 



Ip,p';l ~ 2 Cplp^p-l 



5.5 Ckecks of tke calculation 

We have performed several checks to control the correctness of our results, both 
internal and external, by comparing our results with the known ones. 
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The divergences coming from UV and IR poles all cancel. This is surely one 
of the most important analytical check, that covers the virtual terms and the 
integrals of the soft and coUinear limits of the four-body final states. 



The full calculation, as m— i>0, agrees with the massless results of Ref. [T^. In 



the tables of Chapter ^ this is easily seen. In fact, in the first column, we report 
these massless results, and we can see that our massive calculation reaches the 
massless limit as the ratio m/ E goes to zero. 

Our four- dimensional matrix elements for the processes e^e~ ^ Zj^i ^QQgg 



and e e~ ^ Z/'^ ^ QQQQ agree with Ref. ||Tl| . We have performed a numerical 
comparison between their results and ours. Furthermore, the soft and collinear 
limits of the four-body matrix elements for the process Z/7 — > QQ plus two light 



partons are correctly given by formulae (|5.31 ). 



4. Near the production threshold, we recover the Coulomb singularity. If (3 is the 
velocity of the two massive quarks in the fermion centre-of-mass system, then 
(see Ref. 0) 

o?(t{^/^(xi, X2) ^ [Cf - ^) da':^]^{xi,X2) . (5.44) 

By evaluating {p+p'Y centre of mass of the two massive quarks, for small 

/3, we get 

2 



{p + p 



l\2 



(q-ky = q'{xi + X2-l) . (5.45) 



Choosing for example Xi = X2 we have 

xi=X2 = ^ (1 + P + P/5^ 

By letting /3 get smaller and smaller we have checked that the behaviour of the 
virtual differential cross section is in agreement with eq. ( ^.44| ). 



5. The last check we have made, much more involved than the previous ones, is 



described in Chapter and in Ref. [30|. We give here only a brief sketch of it. 



We have used the fact that the semi-inclusive differential cross section for the 
production of a heavy quark is calculable in perturbative QCD: in fact, the mass 
of the final quark acts as a cut-off for the collinear divergences and logarithms 
of the ratio q^ /m? appears in the final result. 
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On the other hand, using the factorization theorem and the AltareUi-Parisi 
evolution equations, we can obtain an expression of the differential cross section 
in which the large logarithms are correctly resummed, while powers of the ratio 
/q^ are completely neglected. 

We have made an expansion in ag of the resummed expression till order a^, and 
we have checked that some moments of the coefficients of the large logarithms 
are correctly given by our program, in the limit of small masses. 
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Numerical results 



We implemented our analytical result in a FORTRAN program, which behaves like 



a "partonic" Monte Carlo generator, analogous to the program EVENT |jTO|. We 
collect here some results obtained with our code. Since for this kind of calculations it 
would be difficult to perform analytical comparisons, the only possible alternative is 
to choose a few shape variables, and compare numerical results, in the spirit of what 
has been done in Ref. [BTl, for the case of the massless calculation. 



We include in these results only the contributions from cut graphs of Sees. ^]4| 
4.6. 1| and |4.6.2|, in which the weak current couples to the same heavy-flavour loop 



and there is a single QQ pair in the final state, which is the really hard part of the 
calculation. 

For the contributions involving two heavy-quark pairs in the final state, it is easier 
to compare directly the value of the matrix elements squared (this part of our program 
was, in fact, checked in this way with the program of Ref. [pT|). 

We have chosen a set of shape variables for which it should be easy to obtain quite 
accurate numerical results. We have fixed the centre-of-mass energy to be 100 GeV, 
and the mass of the heavy quark has been taken to be equal to 1, 10 and 30 GeV. 
We present separately the results for a hypothetical vector boson with purely axial or 
purely vector couplings, normalized to the massless total cross section at zeroth order 
in Us- We have chosen the following shape variables: the thrust t, the c parameter, 
the mass of the heavy jet squared (according to the thrust axis), the energy- 
energy correlation EEC, the three-jet fractions according to the E, EM [^, JADE, 
and DURHAM schemes. 

For some shape variables, the presence of massive particles in the final state may 
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introduce ambiguities in the definition, owing to the fact that, in the massless case, 
energy and momentum can be interchanged. We thus refer to the exact definitions 
given in Ref. for t, c, M| and in Ref. 0] for the EEC. We collect here these 
definitions. Thrust is defined as 

t = max ^ , (6.1) 

" Z^IPil 

i 

where Pi denotes the three-momentum of the i^^ particle in the centre-of-mass system, 
and the sum extends over all final-state particles. The direction of n that maximizes 
the above quantity is called thrust axis. 

The c parameter is derived from the eigenvalues of the infrared-safe momentum 
tensor 



Ciab 



(6.2) 



where is the a*^ component of the three- momentum pj. If we denote with Ai, A2 
and A3 the eigenvalues, we define 

c = 3(AiA2 + A2A3 + AiA3) . (6.3) 

The definition of the heavy-jet mass M|, according to the thrust axis, is obtained 
with the following procedure: it is the maximum value between the invariant masses 
of the particles belonging to the two different hemispheres separated by the plane 
orthogonal to the thrust axis. 

For t, c, M| and EEC we present moments, instead of distributions, because 
they can be obtained with higher precision. For thrust, for example, we compute, 
according to the notation of Sec. 



a«\ , . fa 



dTy/A (1 - tr = ( K^Ain) + [-^ ) Bi^^M ■ (6.4) 



2 



We further decompose 



^V/A — Bv/A,Ca + Bv/a,Cf + Bv/a,Tf ' (6-5) 
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where the Ca, Cf and Tp subscripts denote the CfCa, Cp and nfCpTp colour com- 
ponents. For the other quantities, moments are defined as 

where the sum runs over all the final particles, and 6ij is the angle between the 
corresponding three-momenta. 

Jet clustering algorithms are defined giving two ingredients: 

- the rule to compute the resolution parameter yij for each pair of particles ij in 
the final state 

- the recombination rule for the two particles. 

Before starting the algorithm, you fix the value of a resolution parameter Hcut to be 
used as a discriminant condition. Then, according to the given rule, you compute the 
resolution parameter j/j^ for each couple of final-state particles. If the minimum values 
of yij 's is less than ycut, then the two particles, for which this value was computed, are 
recombined into one pseudo-particle, and you start again, by calculating a new set of 
yij's. If, instead, no value of y^^'s satisfies the condition yij < i/cut, then the algorithm 
is over, and the final number of pseudo-particles obtained gives the number of jets. 

We have considered four different jet-clustering algorithms: E, EM, JADE and 
DURHAM. Their resolution parameter is defined by 

E 
EM 
JADE 
DURHAM 



Pi ■ Pj 
Vij ~ 2 ^ , 

EE- 

yij = 2 -^{l - cos%) , 

y,, = 2minf^,^j(l-cos%) , (6.6) 
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while their recombination rule is the same for all of them 

Pij = Pi +Pj ■ (6.7) 

Observe that the E scheme is not infrared-safe if ycut < Tn^/q^- In fact, in this 
range, the configuration made up of two heavy quarks plus a soft gluon cannot be 
reduced to two pseudo-particles, since the recombination parameter will fail the cut. 
The cancellation of soft divergences cannot therefore work for these values of the cut 
parameter. 

For the jet clustering algorithms, we have computed 

j dTv/ASN^iy^^,),3 = (g) ^y/A(l/cni) + ( ^ ) ^ ^y/A ( J , (6.8) 

where X stands for one of the jet-clustering algorithms, and NxiVcut) is the number 
of pseudo-particles in the final state after the clustering procedure is over. 

We have chosen the renormalization scale fi = E, and riif = 5. The results are 



given in Tabs. |6.1| to |6.9| . The first column of each table contains the massless limit, 

in 



obtained with the program that has generated the results of Refs. and ^1 
order to allow a comparison with our massive calculation. 

Further results for m/E = 0.2 GeV can be found in Ref. . 
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n 


m/S = 


m/E = 0.01 


rn/£; = 0.1 


m/E = 0.3 






1 


71.56 ±0.026 


71.53 ±0.045 




57.14 ±0.024 


21.83 ±0.006 


2 


5.249 ± 9 10"^ 


5.303 ± 0.005 




4.47 ±0.0034 


1.524 ±8 10-4 


3 


0.7956 ± 1.710--* 


0.8056 ±0.0012 




0.6887 ±8 10-4 


0.2127 ±2 10-4 


4 


0.1582 ±4 10-5 


0.1604 ±3.4 10-"* 




0.1381 ±2.4 10-4 


0.04003 ±6 10-5 


5 


0.03614 ± 110-5 


0.03668 ± 1 10-"^ 




0.03174 ± 710-5 


0.008845 ±2 10-5 






1 


-4.019 ± 0.0023 


-4.34 ±0.06 




-0.704 ±0.026 


3.929 ±0.003 


2 


2.566 ± 3 10-^ 


2.48 ±0.007 




1.78 ± 0.004 


0.699 ± 5 10-4 


3 


0.5486 ± 8 10-5 


0.5348 ± 0.002 




0.3883 ± 0.001 


0.1293 ± 1.210-4 


4 


0.1269 ± 2.510-5 


0.1241 ±6 10-" 




0.09007 ±2.7 10-4 


0.02743 ±3.4 10-5 


5 


0.03186 ± 8 lO-'^ 


0.03116 ± 1.810-4 




0.02255 ±8 10-5 


0.006392 ± 110-5 






1 


-22.52 ±0.0016 


-22.37 ±0.004 




-18.48 ±0.005 


-7.767 ±0.002 


2 


-1.555 ± 1.210-^ 


-1.552 ±6 10-4 




-1.38 ±710-4 


-0.5844 ±3 10-4 


3 


-0.2155 ±2.3 10-5 


-0.2153 ± 1.3 10-4 




-0.1958 ± 1.710-4 


-0.08272 ±7 10-5 


4 


-0.03868 ±6 10-6 


-0.03864 ±3.5 10-5 




-0.03559 ±4 10-5 


-0.01515 ± 1.810-5 


5 


-0.007909 ± 1.810-6 


-0.007897 ± 110-5 




-0.007344 ± 1.310-5 


-0.003174 ±5 10-6 






B\ 






1 


71.56 ±0.026 


71.5 ± 0.045 




54.27 ±0.023 


12.46 ±0.0035 


2 


5.249 ± 9 10-^ 


5.301 ± 0.005 




4.284 ± 0.003 


0.9336 ± 5 10-4 


3 


0.7956 ± 1.710-^ 


0.8054 ±0.0012 




0.6637 ±8 10-4 


0.138 ± 1.210-4 


4 


0.1582 ±4 10-5 


0.1603 ±3.4 10-4 




0.1336 ±2.3 10-4 


0.02708 ±3.6 10-5 


5 


0.03614 ± 110-5 


0.03667 ±110-4 




0.03077 ± 710-5 


0.006162 ±1.1 10-5 






B'a 






1 


-4.019 ±0.0023 


-4.3 ±0.06 




0.5 ±0.025 


3.722 ±0.002 


2 


2.566 ± 3 10-^ 


2.481 ± 0.007 




1.801 ± 0.0036 


0.5694 ±3 10-4 


3 


0.5486 ± 8 10-5 


0.535 ± 0.002 




0.3871 ± 9 10-4 


0.1041 ±8 10-5 


4 


0.1269 ±2.5 10-5 


0.1241 ±610-4 




0.0893 ±2.6 10-4 


0.02206 ± 2 10-5 


5 


0.03186 ±8 10-6 


0.03116 ± 1.810-4 




0.02228 ±8 10-5 


0.005136 ±7 10-6 






B'a 






1 


-22.52 ±0.0016 


-22.36 ±0.004 




-17.5 ±0.005 


-4.457 ±0.0013 


2 


-1.555 ± 1.210-4 


-1.552 ±6 10-4 




-1.311 ±710-4 


-0.3533 ± 1.710-4 


3 


-0.2155 ±2.3 10-5 


-0.2152 ± 1.3 10-4 




-0.1861 ± 1.610-4 


-0.05156 ±4 10-5 


4 


-0.03868 ±6 10-6 


-0.03862 ±3.5 10-5 




-0.03382 ±4 10-5 


-0.009594 ± 110-5 


5 


-0.007909 ± 1.810-6 


-0.007893 ± 110-5 




-0.00697 ± 1.210-5 


-0.002024 ±3 10-6 



Table 6.1: The thrust t. 
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n 


m/E = 


m/E = 0.01 


m/E = 0.1 


m/E = 0.3 






^V,Ca 


in) 




1 


68.14 ±0.026 


68.32 ±0.05 




56.71 ± 0.03 


19.46 ± 0.008 


2 


4.542 ±0.001 


4.591 ± 0.006 




4.014 ±0.004 


1.029 ± 0.0024 


3 


0.6258 ± 2 lO^'' 


0.6336 ±0.0013 




0.5549 ±0.001 


0.1072 ±0.001 


4 


0.1141 ± 5 10"^ 


0.1156 ±3.6 10-'' 




0.1002 ±3 10-" 


0.0167 ± 4.5 IQ-" 


5 


0.0241 ± 1.4 10"^ 


0.02447 ± llO-'' 




0.02084 ±110-" 


0.0037 ± 2 10-" 








(n) 




1 


-21.54 ±0.003 


-21.39 ±0.05 




-10.55 ±0.026 


3.166 ± 0.004 


2 


0.2522 ±4 10"^ 


0.218 ±0.007 




0.4605 ± 0.004 


0.736 ±0.0011 


3 


0.0826 ± 1 lO"'' 


0.0784 ±0.002 




0.119 ±0.001 


0.1834 ±5 10-" 


4 


0.01719 ± 310"^ 


0.0165 ±6 10-'' 




0.02626 ±3 10-" 


0.055 ±2 10-" 


5 


0.003789 ± 1 10"^ 


0.00365 ± 1.810-'' 




0.00619 ±110-" 


0.01829 ± 9 10-5 






V,Tf 


(n) 




1 


-23.28 ±0.0017 


-23.2 ±0.005 




-20.04 ±0.006 


-7.809 ±0.0025 


2 


-1.705 ± 1.310"^ 


-1.704 ±6 10-^ 




-1.563 ±8 10-" 


-0.5961 ±5 10-" 


3 


-0.2507 ±2.5 10-5 


-0.2505 ± 1.4 10-" 




-0.2308 ±2 10-" 


-0.1014 ±2 10-" 


4 


-0.04763 ± 6 10"'^ 


-0.04758 ±3.6 10-5 




-0.04375 ± 5 10-5 


-0.02601 ± 8 10-5 


5 


-0.01032 ± 2 10"'^ 


-0.0103 ± 110-5 




0.009445 ± 1.510-5 


-0.0082 ±3 10-5 








(n) 




1 


68.14 ±0.026 


68.29 ±0.05 




53.84 ±0.03 


10.99 ±0.006 


2 


4.542 ± 0.001 


4.59 ± 0.006 




3.841 ± 0.004 


0.6034 ±0.0018 


3 


0.6258 ± 2 lO"'' 


0.6334 ±0.0013 




0.533 ±0.001 


0.0638 ± 8 10-" 


4 


0.1141 ± 5 lO^-"" 


0.1156 ±3.6 10--* 




0.0965 ±3 10-" 


0.0102 ±3.410-" 


5 


0.0241 ± 1.4 10"^ 


0.02447 ± 110-" 




0.02013 ±110-" 


0.00246 ± 1.4 10-" 








(n) 




1 


-21.54 ±0.003 


-21.35 ±0.05 




-8.804 ±0.025 


3.617 ±0.0026 


2 


0.2522 ±4 10"^ 


0.22 ±0.007 




0.546 ± 0.004 


0.7156 ±7 10-" 


3 


0.0826 ± 1 10"" 


0.0787 ±0.002 




0.1304 ±0.001 


0.1817 ±3 10-" 


4 


0.01719 ±3 10-5 


0.0166 ±6 10-" 




0.02845 ±3 10-" 


0.05504 ± 1.2 10-" 


5 


0.003789 ±110-5 


0.00366 ± 1.810-" 




0.0067 ±9 10-5 


0.01825 ±5 10-5 






^A,Tf 


(n) 




1 


-23.28 ±0.0017 


-23.19 ±0.005 




-19.01 ±0.005 


-4.617 ±0.0015 


2 


-1.705 ± 1.310-^ 


-1.704 ±6 10-" 




-1.491 ±8 10-" 


-0.4038 ±3.6 10-" 


3 


-0.2507 ±2.5 10-5 


-0.2504 ± 1.4 10-" 




-0.2207 ± 1.7 10-" 


-0.07702 ± 1.4 10-" 


4 


-0.04763 ±6 10-6 


-0.04756 ±3.6 10-5 




-0.04188 ±4.5 10-5 


-0.02099 ±5 10-5 


5 


-0.01032 ±2 10-6 


-0.0103 ± 110-5 




0.009048 ± 1.310-5 


-0.00675 ±2 10-5 



Table 6.2: The mass of the heavy jet squared M^. 
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n 


m/i; = 


m/E = 0.01 


m/E = 0.1 


m/E = 0.3 






1 


304.3 ±0.1 


303.3 ±0.18 


233.9 ±0.09 


84.7 ±0.023 


2 


70.37 ±0.011 


70.97 ±0.06 


58.25 ± 0.03 


19.18 ±0.007 


3 


30.21 ± 0.006 


30.55 ± 0.03 


25.56 ±0.02 


7.789 ± 0.004 


4 


16.16 ±0.0034 


16.36 ±0.02 


13.85 ±0.012 


4.008 ±0.003 


5 


9.658 ±0.0022 


9.784 ±0.012 


8.341 ±0.008 


2.335 ± 0.002 






1 


-35.42 ±0.01 


-36.2 ±0.2 


-13.58 ±0.1 


13.41 ±0.013 


2 


31.69 ±0.003 


30.42 ± 0.07 


20.29 ± 0.036 


7.845 ± 0.004 


3 


20.7 ±0.002 


20.1 ± 0.04 


13.69 ±0.02 


4.23 ± 0.003 


4 


13.52 ±0.0014 


13.19 ±0.025 


9.04 ±0.013 


2.479 ±0.002 


5 


9.228 ±0.001 


9.024 ±0.017 


6.185 ±0.009 


1.543 ±0.0013 






1 


-95.95 ±0.007 


-95.12 ±0.016 


-75.92 ±0.02 


-30.07 ±0.009 


2 


-20.85 ±0.0016 


-20.81 ± 0.006 


-18.02 ±0.007 


-7.206 ± 0.003 


3 


-8.125 ±8 10"* 


-8.117 ±0.003 


-7.227 ±0.004 


-2.915 ±0.0018 


4 


-3.862 ±5 10-'* 


-3.859 ±0.002 


-3.495 ±0.0027 


-1.435 ±0.0012 


5 


-2.009 ± 3 lO""* 


-2.007 ±0.0013 


-1.842 ±0.0018 


-0.7771 ± 8 lO-'' 






B% 






1 


304.3 ±0.1 


303.2 ±0.18 


222 ±0.09 


48.14 ±0.014 


2 


70.37 ±0.011 


70.95 ± 0.06 


55.73 ± 0.03 


11.58 ±0.0045 


3 


30.21 ± 0.006 


30.54 ±0.03 


24.59 ±0.018 


4.942 ± 0.0027 


4 


16.16 ±0.0034 


16.36 ±0.02 


13.36 ±0.01 


2.643 ±0.0018 


5 


9.658 ±0.0022 


9.782 ±0.012 


8.071 ± 0.007 


1.587 ±0.0013 






Bl 






1 


-35.42 ±0.01 


-36.04 ±0.2 


-8.24 ±0.1 


13.16 ±0.008 


2 


31.69 ±0.003 


30.44 ±0.07 


20.63 ± 0.034 


6.425 ±0.0027 


3 


20.7 ±0.002 


20.11 ±0.04 


13.64 ±0.02 


3.394 ±0.0017 


4 


13.52 ±0.0014 


13.19 ±0.025 


8.94 ±0.013 


1.975 ±0.001 


5 


9.228 ±0.001 


9.024 ±0.017 


6.091 ± 0.009 


1.223 ±810"'' 




bXtM 


1 


-95.95 ±0.007 


-95.07 ±0.016 


-71.88 ±0.02 


-17.2 ±0.005 


2 


-20.85 ± 0.0016 


-20.8 ±0.006 


-17.11 ±0.007 


-4.312 ±0.0017 


3 


-8.125 ±8 10"* 


-8.113 ±0.003 


-6.865 ± 0.004 


-1.788 ±0.001 


4 


-3.862 ±5 10"'* 


-3. 857 ±0.002 


-3.317 ±0.0025 


-0.8903 ±7 10-4 


5 


-2.009 ± 3 10-"* 


-2.006 ±0.0013 


-1.745 ±0.0017 


-0.4836 ±4.5 10-4 



Table 6.3: The c parameter. 
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Ckapter 6. Numerical results 



n 


m/E = 


m/E = Om 


m/E = Q.l 


m/E = 0.3 






pEEC 

-'-^ .(• 






A 

u 


202.8 ±0.07 


202.3± 0.12 


154.4 ±0.06 


46.81 ±0.013 


1 
1 


142.6 ±0.06 


143.1 ±0.1 


117 ±0.05 


38.2 ±0.01 


z 


115.3 ±0.05 


115.9 ±0.09 


97.55 ±0.04 


33.55 ±0.009 


q 
O 


99.28 ± 0.05 


99.79 ± 0.08 


85.27 ±0.04 


30.38 ± 0.008 


A 


88.43 ± 0.05 


88.92 ±0.08 


76.65 ± 0.04 


28.01 ± 0.007 


5 


80.48 ± 0.05 


80.95 ±0.07 


70.18 ±0.04 


26.13 ±0.007 






A 
U 


-23.62 ±0.006 


-24.13 ±0.14 


-8.83 ±0.07 


8.224 ±0.007 


1 
1 


-9.022 ± 0.006 


-9.82 ±0.12 


-3.46 ±0.05 


6.861 ±0.006 




z 


-5.268 ± 0.006 


-5.96 ±0.1 


-1.79 ±0.05 


5.99 ±0.005 


q 
O 


-3.717 ±0.006 


-4.33 ±0.1 


-1.06 ±0.04 


5.387 ±0.005 


A 


-2.899 ±0.005 


-3.46 ±0.1 


-0.664 ±0.04 


4.937 ±0.0045 


5 


-2.399 ±0.005 


-2.92 ±0.1 


-0.43 ±0.04 


4.585 ±0.004 






A 
U 


-63.97 ±0.005 


-63.43 ±0.01 


-50.32 ±0.014 


-17.12 ±0.005 


1 
i 


-45.4 ± 0.004 


-45.25 ±0.009 


-38.32 ±0.01 


-13.98 ±0.004 




z 


-36.92 ±0.0034 


-36.83 ± 0.008 


-32.05 ±0.01 


-12.27 ±0.0035 


q 
o 


-31.87 ±0.003 


-31.8 ±0.008 


-28.07 ±0.009 


-11.11 ±0.003 


A 


-28.45 ±0.003 


-28.39 ±0.007 


-25.26 ±0.008 


-10.23 ±0.003 


5 


-25.93 ±0.003 


-25.88 ±0.007 


-23.15 ±0.008 


-9.544 ±0.0027 




Bf^SM^) 


A 
U 


202.8 ±0.07 


202.2 ±0.12 


146.6 ± 0.06 


26.66 ±0.008 


1 
1 


142.6 ± 0.06 


143 ±0.1 


111 ±0.05 


21.72 ±0.006 




z 


115.3 ±0.05 


115.8 ±0.09 


92.62 ± 0.04 


19.05 ± 0.006 


q 
o 


99.28 ±0.05 


99.75 ±0.08 


80.96 ±0.04 


17.23 ±0.005 


A_ 


88.43 ±0.05 


88.87 ±0.08 


72.77 ±0.04 


15.87 ±0.005 


5 


80.48 ±0.05 


80.91 ±0.07 


66.62 ±0.035 


14.79 ±0.005 




Bl%{n,Q) 


n 
u 


-23.62 ±0.006 


-24.04 ±0.14 


-5.3 ±0.06 


8.116 ±0.0045 


1 
1 


-9.022 ± 0.006 


-9.76 ±0.12 


-0.89 ±0.05 


6.709 ± 0.004 





-5.268 ±0.006 


-5.9 ±0.1 


0.32 ±0.045 


5.85 ±0.003 


q 
O 


-3.717 ±0.006 


-4.28 ±0.1 


0.78 ±0.04 


5.259 ±0.003 


4 


-2.899 ±0.005 


-3.41 ±0.1 


0.98 ±0.04 


4.82 ±0.003 


5 


-2.399 ±0.005 


-2.88 ±0.1 


1.07 ±0.04 


4.476 ± 0.003 









-63.97 ±0.005 


-63.4 ±0.01 


-47.66 ±0.013 


-9.928 ±0.003 


1 


-45.4 ±0.004 


-45.22 ±0.009 


-36.29 ±0.01 


-8.091 ±0.0024 


2 


-36.92 ±0.0034 


-36.81 ±0.008 


-30.35 ±0.009 


-7.085 ±0.002 


3 


-31.87 ±0.003 


-31.79 ±0.008 


-26.57 ±0.008 


-6.401 ± 0.002 


4 


-28.45 ± 0.003 


-28.38 ± 0.007 


-23.91 ± 0.008 


-5.889 ±0.0017 


5 


-25.93 ± 0.003 


-25.87 ±0.007 


-21.91 ±0.007 


-5.486 ± 0.0016 



TaKlp 4- TVip PTiprcrv— PTiprcrv pnrrplaf.inn P.P.r^ h — 
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n 


m/i; = 


m/£; = 0.01 


m/E = 0.1 


m/E = 0.2, 






A 
U 


-25.89 ±0.006 


-26.06 ±0.025 


-21.45 ±0.014 


-7.006 ±0.003 


1 
1 


-9.426 ±0.004 


-9.56 ±0.016 


-8. 967 ±0.01 


-3.152 ± 0.0023 


Z 


-4.75 ±0.004 


-4.827 ±0.013 


-4.892 ±0.008 


-1.898 ±0.002 


o 
o 


-2.844 ± 0.004 


-2.894 ±0.011 


-3.086 ±0.008 


-1.316 ±0.0017 


A 


-1.894 ±0.0036 


-1.93 ±0.01 


-2.133 ±0.007 


-0.9876 ±0.0015 


5 


-1.354 ±0.0035 


-1.38 ±0.01 


-1.57 ±0.007 


-0.779 ±0.0014 






A 
U 


8.722 ±0.0013 


8.73 ±0.03 


3.607 ±0.015 


-2.075 ±0.002 


1 
1 


0.427 ±0.0012 


0.593 ±0.02 


0.272 ±0.01 


-1.345 ±0.0014 




z 


-0.6778 ±0.001 


-0.588 ±0.02 


-0.327 ±0.009 


-0.9457 ±0.0012 


q 
O 


-0.81 ±0.001 


-0.76 ±0.018 


-0.441 ± 0.008 


-0.7102 ±0.001 


A 


-0.7574 ±0.001 


-0.73 ±0.017 


-0.437 ±0.007 


-0.5588 ±9 10-4 


5 


-0.6728 ±0.001 


-0.66 ±0.016 


-0.403 ±0.007 


-0.4546 ±8 10-4 






A 
U 


10.77 ± 7 lO-'^ 


10.75 ± 0.0026 


9.205 ±0.003 


3.381 ±0.0013 


1 
1 


4.462 ± 5 lO-'^ 


4.47 ± 0.0018 


4.256 ± 0.0023 


1.718 ±9 10-4 




z 


2.494 ± 4 10"^ 


2.5 ± 0.0014 


2.506 ± 0.002 


1.108 ±7 10-4 


q 
o 


1.625 ±3.6 10-'* 


1.629 ±0.0012 


1.679 ±0.0016 


0.7997 ±6 10-4 


A 


1.16 ±3.4 10-'' 


1.162 ±0.001 


1.218 ±0.0014 


0.616 ±5 10-4 


5 


0.8785 ± 3 lO-"* 


0.8805 ±0.001 


0.933 ±0.0013 


0.4948 ±5 10-4 




Bf^SMA) 


A 
U 


-25.89 ±0.006 


-26.05 ±0.025 


-20.4 ± 0.014 


-4.12 ±0.002 


1 
1 


-9.426 ±0.004 


-9.555 ±0.016 


-8.558 ±0.01 


-1.907 ±0.0014 


o 
z 


-4.75 ± 0.004 


-4.825 ±0.013 


-4.68 ± 0.008 


-1.165 ±0.0012 


q 
o 


-2.844 ±0.004 


-2.894 ±0.011 


-2.957 ±0.007 


-0.8135 ±0.001 


A 


-1.894 ±0.0036 


-1.93 ±0.01 


-2.047 ±0.007 


-0.6129 ±9 10-4 


5 


-1.354 ±0.0035 


-1.38 ±0.01 


-1.508 ±0.006 


-0.4846 ±9 10-4 




Bl%{n,l) 


A 
U 


8.722 ±0.0013 


8.71 ± 0.03 


2.855 ±0.014 


-1.98 ±0.0012 


1 
1 


0.427 ±0.0012 


0.586 ± 0.02 


-0.025 ±0.01 


-1.228 ±9 10-4 





-0.6778 ±0.001 


-0.59 ±0.02 


-0.486 ±0.008 


-0.853 ±7 10-4 


q 
O 


-0.81 ±0.001 


-0.763 ±0.018 


-0.541 ±0.007 


-0.6373 ±6 10-4 


4 


-0.7574 ±0.001 


-0.732 ±0.017 


-0.506 ±0.007 


-0.4999 ±6 10-4 


5 


-0.6728 ±0.001 


-0.66 ±0.016 


-0.454 ± 0.007 


-0.4059 ±5 10-4 




Bf%{n,l) 





10.77 ±7 10"^ 


10.75 ±0.0026 


8.765 ± 0.003 


2.054 ±7 10-4 


1 


4.462 ± 5 10-4 


4.468 ±0.0018 


4.06 ±0.002 


1.07 ± 510-4 


2 


2.494 ±410-'' 


2.499 ±0.0014 


2.393 ±0.0018 


0.6961 ±4 10-4 


3 


1.625 ±3.6 10-4 


1.628 ±0.0012 


1.603 ±0.0015 


0.5043 ±3.6 10-4 


4 


1.16 ±3.4 10-4 


1.162 ±0.001 


1.164 ±0.0013 


0.3888 ±3 10-4 


5 


0.8785 ± 3 10-4 


0.88 ±0.001 


0.8916 ±0.0012 


0.3123 ±3 10-4 
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Ckapter 6. Numerical results 



Vcut 


m/E = 


m/E = 0.01 


m/E = 0.1 


m/E = 0.3 




Bv,Ca iVcut) 


0.01 


1343 ± 1.5 


1357 ± 1.5 






0.05 


336.5 ± 0.34 


341.3 ±0.5 


374 ± 0.3 




0.10 


131.7 ±0.15 


133.7 ±0.3 


138.6 ±0.17 


459.5 ±0.17 


0.15 


60.06 ±0.07 


60.75 ±0.18 


63.1 ±0.13 


79.06 ±0.04 


0.20 


26.66 ± 0.04 


27.13±0.16 


29 ±0.12 


24.54 ±0.027 




Bv,Cp iVcut) 


0.01 


-225.5 ± 1 


-88 ±4 






0.05 


129.8 ±0.15 


129.8 ± 1.3 


249.1 ±0.36 




0.10 


73.59 ±0.06 


71.8 ±0.7 


120.1 ±0.2 


47.95 ±0.13 


0.15 


37.61 ±0.035 


37.5 ±0.5 


53.8 ±0.13 


41.25 ±0.025 


0.20 


17.97 ±0.02 


18.1 ±0.34 


24.26 ±0.13 


21.58 ±0.014 




Bv,TFiycut) 


0.01 


-450 ±0.1 


-452.3 ±0.13 






0.05 


-103.7 ±0.03 


-103.8 ±0.04 


-121.3 ±0.05 




0.10 


-38.29 ±0.012 


-38.34 ±0.02 


-42.14 ±0.03 


-154.4 ±0.07 


0.15 


-16.54 ±0.007 


-16.56 ±0.015 


-17.79 ±0.02 


-28. 17± 0.012 


0.20 


-6.869 ±0.004 


-6.871 ±0.009 


-7.424 ±0.013 


-9.268 ±0.008 






0.01 


1343 ± 1.5 


1356 ± 1.5 






0.05 


336.5 ± 0.34 


341.2 ±0.5 


355.5 ±0.3 




0.10 


131.7 ±0.15 


133.6 ±0.3 


132.8 ±0.17 


254.5 ±0.09 


0.15 


60.06 ±0.07 


60.9 ±0.24 


60.83 ±0.13 


46.55 ±0.03 


0.20 


26.66 ± 0.04 


27.1 ±0.14 


27.83 ±0.08 


15.4 ±0.016 




^Xcf iVcut) 


0.01 


-225.5 ± 1 


-88 ±3.6 






0.05 


129.8 ±0.15 


130 ± 1.3 


243.1 ±0.3 




0.10 


73.59 ± 0.06 


71.8 ±0.7 


117.2 ±0.2 


49.67 ±0.07 


0.15 


37.61 ±0.035 


37.45 ±0.5 


52.54±0.12 


31.7±0.016 


0.20 


17.97 ±0.02 


18.1 ±0.33 


23.8 ±0.09 


17.24 ±0.01 






0.01 


-450 ±0.1 


-452.1 ±0.13 






0.05 


-103.7 ±0.03 


-103.8 ±0.04 


-115.4 ±0.045 




0.10 


-38.29 ±0.012 


-38.32 ±0.02 


-40.23 ±0.027 


-86.57 ±0.04 


0.15 


-16.54 ±0.007 


-16.55 ±0.015 


-17.01 ±0.02 


-17.28 ±0.008 


0.20 


-6.869 ± 0.004 


-6.868 ±0.009 


-7.098 ±0.012 


-6.237 ±0.005 



Table 6.6: The E clustering algorithm. 
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Vcut 


m/E = 


m/£; = 0.01 


m/E = 0.1 


m/E = 0.3 






0.01 


1337 ± 1.5 


1340 ± 1.7 


1135 ±0.8 


455.7 ±0.18 


0.05 


329.4 ± 0.3 


333.2 ±0.5 


291.2 ±0.4 


101.6 ±0.07 


0.10 


124.7 ±0.12 


126.3 ±0.3 


111.8 ±0.23 


29.24 ±0.04 


0.15 


53.69 ±0.06 


54.24 ±0.3 


48 ±0.17 


8.176 ±0.027 


0.20 


21.68 ±0.03 


22.02 ±0.2 


18.85 ±0.14 


1.52 ±0.024 






0.01 


-320.3 ±0.9 


-327 ±4 


-193.4 ± 1.5 


37.07 ±0.14 


0.05 


78.25 ±0.12 


74.5 ± 1.4 


75.8 ±0.6 


43.09 ±0.04 


0.10 


42.74 ±0.05 


39.6 ±0.8 


40.85 ±0.3 


20.79 ± 0.04 


0.15 


18.84 ±0.024 


18.5 ±0.5 


18.47 ±0.27 


9.353 ±0.015 


0.20 


7.15 ±0.014 


7.46 ±0.3 


6.45 ±0.5 


3.57± 0.011 




Bv^AVcut) 


0.01 


-454 ±0.1 


-453.4 ±0.14 


-391.2 ±0.13 


-157±0.07 


0.05 


-106.9 ±0.03 


-106.9 ±0.04 


-98.04 ±0.05 


-36.9 ±0.015 


0.10 


-40.73 ± 0.012 


-40.71 ± 0.025 


-37.79 ±0.03 


-12.24 ±0.01 


0.15 


-18.26 ±0.007 


-18.25 ±0.02 


-16. 7± 0.022 


-4.578 ±0.006 


0.20 


-7.976 ±0.004 


-7.97 ±0.01 


-7.125 ±0.013 


-1.53 ±0.005 




BTx:a iVcut) 


0.01 


1337 ± 1.5 


1339 ± 1.7 


1072 ± 0.9 


250.8 ±0.14 


0.05 


329.4 ± 0.3 


332.8 ±0.5 


277.2 ± 0.35 


58.88 ±0.07 


0.10 


124.7 ±0.12 


126.3 ±0.3 


107.1 ±0.26 


17.97 ±0.03 


0.15 


53.69 ±0.06 


54.2 ±0.2 


46.08 ±0.2 


5.185 ±0.02 


0.20 


21.68 ±0.03 


21.97±0.2 


17.9 ±0.25 


1.013 ±0.014 




Baxjf (^™t) 


0.01 


-320.3 ±0.9 


-326.5 ±4 


-164.4 ± 1.2 


43.13 ±0.07 


0.05 


78.25 ±0.12 


74.6 ± 1.4 


78.6 ±0.5 


33.78 ± 0.022 


0.10 


42.74 ±0.05 


39.5 ±0.8 


40.8 ±0.6 


17.18 ±0.013 


0.15 


18.84 ±0.024 


18.4±0.5 


18.7 ±0.2 


8.532 ±0.012 


0.20 


7.15 ±0.014 


7.36 ±0.35 


7.29 ±0.14 


3.529 ±0.007 




BA^piycut) 


0.01 


-454 ±0.1 


-453.1 ±0.14 


-370.7 ±0.13 


-88.5 ±0.04 


0.05 


-106.9 ±0.03 


-106.8 ±0.04 


-93.35 ±0.04 


-22.27 ±0.009 


0.10 


-40.73 ± 0.012 


-40.69 ± 0.024 


-36.04 ±0.03 


-8.066 ±0.006 


0.15 


-18.26 ±0.007 


-18.23 ±0.017 


-15.98 ±0.02 


-3.347 ±0.004 


0.20 


-7.976 ±0.004 


-7.965 ±0.01 


-6.83 ±0.013 


-1.247 ±0.0035 



Table 6.7: The EM clustering algorithm. 
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Ckapter 6. Numerical results 



Vent 


m/E = 


m/£; = 0.01 


m/E = 0.1 


m/E = 0.3 




^v,ca 


0.01 


1352 ± 1.1 


1352 ± 2 


1088 ± 1.6 


389.2 ±0.26 


0.05 


326 ± 0.25 


328.6 ±0.6 


275 ± 0.5 


82.34 ±0.09 


0.10 


122.7 ±0.1 


124.2 ±0.3 


106 ±0.26 


24.91 ± 0.05 


0.15 


53.55 ±0.05 


54.1 ±0.3 


46.84 ±0.26 


8.675 ±0.03 


0.20 


22.46 ±0.03 


22.86 ±0.2 


19.84 ± 0.14 


2.98 ±0.03 






0.01 


-790 ±0.7 


-787 ±5 


-423.2 ± 1.5 


24.27 ± 0.14 


0.05 


14.25 ±0.1 


11.7± 1.4 


25.2 ±0.7 


29.65 ±0.04 


0.10 


30.94 ± 0.04 


27.8 ± 0.8 


29.25 ±0.33 


15.33 ±0.027 


0.15 


17.81 ±0.024 


17.45 ±0.5 


16.64 ±0.24 


7.885 ±0.018 


0.20 


8.396 ±0.015 


8.7 ±0.3 


7.7 ±0.14 


3.796 ±0.012 






0.01 


-471.8 ±0.1 


-470 ±0.16 


-380.9 ±0.16 


-138.5 ±0.06 


0.05 


-111 ±0.03 


-110.8 ±0.05 


-96.21 ±0.05 


-32.06 ±0.015 


0.10 


-41.71 ±0.012 


-41.68 ±0.025 


-37.65 ±0.03 


-11.21 ±0.01 


0.15 


-18.35 ±0.007 


-18. 34± 0.017 


-16.84 ±0.025 


-4.688 ±0.007 


0.20 


-7.839 ±0.004 


-7.834 ±0.01 


-7.252 ±0.013 


-1.955 ±0.005 






0.01 


1352 ± 1.1 


1350 ±2 


1027 ± 1 


214.1 ±0.3 


0.05 


326 ± 0.25 


328.6 ±0.6 


261.9 ±0.4 


47.69 ±0.05 


0.10 


122.7 ±0.1 


124.1 ±0.3 


101.6 ±0.3 


15.35 ±0.03 


0.15 


53.55 ±0.05 


54.15 ±0.27 


45.25 ±0.2 


5.72 ±0.02 


0.20 


22.46 ±0.03 


23 ±0.2 


19.12 ±0.16 


2.113 ±0.02 






0.01 


-790 ±0.7 


-785 ± 5 


-384 ± 1.8 


33.72 ±0.1 


0.05 


14.25 ± 0.1 


11.7± 1.4 


30.6 ±0.5 


24.55 ±0.04 


0.10 


30.94 ±0.04 


27.9 ±0.8 


30 ± 0.36 


12.84 ±0.02 


0.15 


17.81 ±0.024 


17.3 ±0.5 


17.35 ±0.23 


6.91 ±0.01 


0.20 


8.396 ±0.015 


8.7 ±0.3 


8.05 ±0.17 


3.461 ±0.01 




B'/:^f(ycut) 


0.01 


-471.8 ±0.1 


-469.7 ±0.16 


-360.5 ±0.2 


-78.07 ±0.036 


0.05 


-111 ±0.03 


-110.8 ±0.05 


-91.5 ±0.045 


-19.17 ±0.009 


0.10 


-41.71 ±0.012 


-41.66 ±0.025 


-35.93 ±0.03 


-7.253 ±0.006 


0.15 


-18.35 ±0.007 


-18.33 ±0.017 


-16.1 ±0.02 


-3.29 ±0.004 


0.20 


-7.839 ±0.004 


-7.826 ±0.01 


-6.96 ±0.014 


-1.476 ±0.0035 



Table 6.8: The JADE clustering algorithm. 
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Ucut 


m/E = 


m/E = 0.01 


m/E = 0.1 


m/E = 0.3 






0.01 


442.5 ± 0.3 


443 ± 1.6 


329.9 ±0.5 


71.09 ±0.07 


0.05 


114.7 ±0.1 


116 ±0.3 


97.26 ±0.22 


12.98 ±0.05 


0.10 


45.05 ±0.06 


45.25 ±0.2 


39.53 ±0.2 


3.62 ±0.03 


0.15 


20.2 ±0.03 


20.4 ±0.18 


17.82 ±0.13 


1.16 ±0.026 


0.20 


8.593 ±0.02 


8.7 ±0.12 


7.61 ±0.1 


0.41 ±0.02 






0.01 


-126.5 ± 0.12 


-128.3 ±2 


-51.2 ±0.6 


23.01 ±0.04 


0.05 


6.29 ±0.04 


3.8 ± 1 


11.9 ±0.4 


10.58 ±0.024 


0.10 


6.625 ± 0.02 


6.1 ±0.5 


8.84 ± 0.2 


5.444 ±0.017 


0.15 


3.65 ±0.015 


4 ±0.3 


4.6 ±0.15 


2.824 ±0.01 


0.20 


1.658 ±0.01 


1.44 ±0.25 


2.09 ±0.1 


1.039 ±0.009 






0.01 


-163 ±0.03 


-162.5 ±0.07 


-127.3 ±0.06 


-31.63 ±0.014 


0.05 


-42.7 ±0.012 


-42.66 ±0.026 


-37.61 ± 0.034 


-7.768 ± 0.008 


0.10 


-17.39 ±0.007 


-17.37 ±0.015 


-15.8 ±0.02 


-2.996 ±0.006 


0.15 


-8.016 ±0.004 


-8.024 ±0.011 


-7.353 ±0.014 


-1.302 ±0.005 


0.20 


-3.502 ±0.0024 


-3.503 ±0.008 


-3.23 ±0.01 


-0.427 ±0.004 




B^^^Jycut) 


0.01 


442.5 ±0.3 


443.6 ± 1 


313 ±0.6 


41.43 ±0.06 


0.05 


114.7 ±0.1 


115.4 ±0.6 


93 ±0.3 


8.37 ±0.024 


0.10 


45.05 ±0.06 


45.07 ±0.22 


37.75 ±0.18 


2.554 ±0.02 


0.15 


20.2 ±0.03 


20.48 ±0.16 


17.2 ±0.14 


0.887± 0.016 


0.20 


8.593 ±0.02 


8.77 ±0.12 


7.3 ±0.1 


0.34 ±0.018 




BT^^ycut) 


0.01 


-126.5 ±0.12 


-127 ±2 


-42.3 ±0.7 


21.21 ±0.03 


0.05 


6.29 ±0.04 


3.9 ±0.9 


14.5 ±0.35 


9.74 ±0.017 


0.10 


6.625 ± 0.02 


6 ±0.5 


9.77 ±0.2 


5.213 ±0.011 


0.15 


3.65 ±0.015 


4.04 ±0.3 


4.84 ±0.2 


2.768 ±0.007 


0.20 


1.658 ±0.01 


1.47 ±0.25 


2.23 ±0.1 


1.012 ±0.006 




B^^Aycut) 


0.01 


-163 ±0.03 


-162.4 ±0.06 


-121.3 ±0.05 


-19.36 ±0.009 


0.05 


-42.7 ±0.012 


-42.65 ± 0.026 


-36 ±0.03 


-5.433 ±0.005 


0.10 


-17.39 ±0.007 


-17.36 ±0.015 


-15.15 ±0.02 


-2.345 ± 0.004 


0.15 


-8.016 ±0.004 


-8.026 ±0.012 


-7.06 ±0.013 


-1.097 ±0.0034 


0.20 


-3.502 ±0.0024 


-3.499 ±0.008 


-3.093 ±0.008 


-0.369 ±0.0025 



Table 6.9: The DURHAM clustering algorithm. 
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Conclusion 



In this thesis, we have described a next-to-leading-order calculation of the cross section 
for the heavy-flavour production in e^e~ collisions, including quark mass effects. The 
computation of the amplitudes was performed in a fully analytical way, and the 
different transition amplitudes are available as FORTRAN subroutines. 

The complete calculation is implemented in a program which behaves like a "par- 
tonic" Monte Carlo generator, in which pairs of weighted correlated events are pro- 



duced, according to the subtraction procedure described in Sec. pAl . The program can 
compute every kind of differential distribution for unoriented infrared-safe quantities, 
and the results can be plotted in histograms for future analysis. At the same time, it 
implements the calculation of the three-jet decay rate according to the E, EM, JADE 
and DURHAM jet-clustering algorithms. 

Instead of publishing figures illustrating the differential behavior of the shape 
variables, we have preferred to compute some moments of these distributions, because 
these results can be obtained with high accuracy and the comparison with other 
groups is easier. 

In this spirit, we have calculated the first five moments of the thrust, the heavy-jet 
mass, the c parameter and the energy-energy correlation, for the ratio m/E equal 
to 0.01, 0.1 and 0.3. We have also computed the three-jet decay rate with a cut 
parameter ranging from 0.01 to 0.2. We have calculated separately these quantities 
for a hypothetical vector boson with purely axial or purely vector couplings, and we 
have collected our results in the tables of Chapter |. In this way, we have succeeded 
to perform a partial comparison between our results for the jet-clustering algorithms 
with those of Ref. 0, and we have found satisfactory agreement. 

Some other applications of our calculation have appeared in the literature. 
In Ref. 1^, we have studied the momentum correlation in the decay of the Z/7 boson 
into a couple of heavy quarks, at order a^. We have found small corrections, of order 
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Conclusion 



of 1%, to this quantity, so that, non-perturbative effects, of order K/E (where A is a 
typical hadronic scale) may compete with the perturbative result. 

In Ref. |3^, we have given a verification of several ingredients that are used in 
the fragmentation-function formalism for heavy-quark production, such as the initial 
conditions, computed in Ref. |T^, and the NLO splitting functions in the time-like 
region pO| . This check has been a further test of the validity of our massive calcula- 
tion. In addition, by merging our fixed order cross section with the NLO resummed 
one, we have derived an improved differential cross section, which is accurate at NLO 
for E ^ m, but at next-to-leading-log for E ^ m. 



Appendix A 



Useful integrals 



In this appendix, we want to introduce some useful functions and integrals, with their 
properties, which have been extensively used throughout our computations. 

The Gamma function is defined by 



fOC 

r(z) 



poo 

/ dxe-^x'-^ Re^>0, (A.l) 
Jo 

and, from this definition, you can demonstrate 

fdxx'^{l-xf^2pdct> (Sin0f--^^(cos0f^^^^ ^^" + ^^;;^f . (A.2) 
Jo Jo i (a + p + zj 

A useful expansion is given by 

r(i + 6) = i-7^6+^^i±^6^ + o(6^) , (A.3) 

where 7^ = 0.5772157 ... is the Euler-Mascheroni constant. 
The dilogarithm function is given by 

log(l-^) 



(,)^_rrf,Mi^ ,<1, (A.4) 
Jo z 



Li 

/o 



and an immediate consequence of this definition is the following expansion in powers 
of e 

dxx-^-^'{l-axf' = £ dxx-^-'''[l + (3elog{l-ax) + 0{e^)' 

= -— -/3eLi2 + (e^) . (A.5) 
7 e ^ ^ 
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Appendix A. Useful integrals 



One of the most used properties is the analytic continuation of the dilogarithm func- 
tion 

Li2 {x ±irj) — — Li2 f — ^ — ^ log^ x + ^ ± ivr log x x > 1 , (A. 6) 

\xj 2 3 

that can be demonstrated with the help of 

log{— X ±ir]) — log X ± in x>0. (A. 7) 

In addition 

/ dxx'^'"^- —- = (1 + 0;)'/ x] a > -1 , 

where we have used the projective transformation 



(A.8) 



X ^1 ^ . (A.9) 

When calculating loop integrals, it is often useful to employ the Feynman parametriza- 
tion, which, in the most general form, reads 

i=l \k=l J 

where q are arbitrary complex numbers and 



c 



^Q. (A.ll) 



i=l 



For our future purpose, we need to compute the following Feynman integral 

d'^l 1 1 



f (TL 

J {2nY {I + pi)"^ — ml + it] {I + pi + p2)'^ — ml + it] 

^ (A.12) 



(/ + Pi + P2 H \- Pnf -ml + if] ' 

where {+irj), with rj > 0, gives the prescription of how the contour integral has to be 
deformed around the poles, and 

f^Pi-O. (A.13) 



1=1 
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We can look at / as being 



with Cj = 1 [i = 1,2, . . . ,n), c = n and Di being the single denominators in eq. ( |A.12|) . 
Using the identity ( |A.l(j| ) and integrating over I, we can rewrite / in the following 

way 



' r 



:-ir'T(„-i)f^, (A.15) 



(47r)(^) 
(47r 



2J \ ^ -I JJ 2 



where 

D = -^aiUj Sij + ^aiml - 17] , (A.16) 

i>j i 

and Sij is the square of the momentum flowing through the i-j cut of the diagram 
representing /. 
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Appendix B 



One-loop scalar integrals 



B.l Kinematical invariants 



We introduce now the following kinematical invariants (not independent) 



Sl 


= (g — p')'^ = (k + pY > rn^ 


(B.l) 


S2 


= (g — pY = {k + p'Y > 


(B.2) 


S3 


= {q - kf = {p + p'f > Arn^ 


(B.3) 


CTl 


= (g — p'Y — = g^(l — X2) 


(B.4) 


0-2 


= (g — pY — = g^(l — Xi) 


(B.5) 


0^3 


= iq-kf = q\l-Xg), 


(B.6) 



where xi and X2 are defined by ( p.6|) and Xg by (|3.38 



We can classify the different types of scalar integrals according to the number of 
massive propagators in the loop and according to the "shape" of the loop: box (B) 
and triangle (T). 
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Appendix B. One- loop scalar integrals 




Figure B.l: Diagram representing B2m- 



B.2 Box witk two massive propagators: 

We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type depicted in Fig. [B.l 



B. 



2m 



dH 1 



1 



(27r)^ P {l-kf {l+p-qy-m'^ {I + pf - m? 
Using eqs. ( |A.15| ) and ( [A.16| ), we can rewrite this integral as 

i . . X /■ [da] 



B2m = i^rcY r(2 + e) I 



D 



2+e 



d = A-2e 



(B.7) 



fB. 



where 
D 



tti q;2 m + ^3 [q — p) + ^2 ^3 g + ^2 «4 [p + k) + 03 0:4 
+ (0:2 + 03) ni^ -ir] . 
According to the definitions given in eqs. (|B.l|) -( pl6| ), D becomes 

D = — m^(ai a2 + — a2 — as) + S2 oi\ 0^3 + (^ol2 0:3 + s\ 0.2 + ir] 



(B.9) 



-m^(a2 + as)^ + (■§2 - m'^)ai 03 + g^a2 "3 + (si - m^)a2 «4 + iv 



.(B.IO) 



In order to perform the integration, it is useful to continue the denominator in a 
region where it is a negative-definite form. For this reason, we assume < 0, and 



B.2. Box witk two massive propagators: -62™ 
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we perform the integration of a function that now has no poles on the integration 
contour. At the end, we return to the physical region according to the irj prescription 

= -(m^ -iri) >0 . 



With this definition, we have 
D = 

and eq. ( [B.8|) becomes 



B 



2m 



{Any r(2 + e 



M^(a2 + "3)^ + o"2 Oil as + q^a2 03 + c^i 0^2 "4 

[da] \ 



(B.ll) 



fB.12) 



1 



-l-iriyj [-D] 



167r2 



{AtxY r(2 + e) e'^'I 



[B.13) 



We make the following change of variables 

ai = (1 - y)z 

a2 = xy 

"3 = (1 - x)y 

04 = (1 - y)il - z) 

where the Jacobian of the transformation is 

d{ai 020:3) 



(B.14) 



d{x y z) 

In this way, we can write I of eq. ( |B.13| ) as 

/ = / dx dy dz y (1 — y) 
Jo 

X 



1 



]\42y2 _|_ (X2yz{l — x){l — y) + q'^xy'^{l — x) + aixy{l — y){l — z) 

/ dx dy dz y~^~''{l — y) 
Jo 



2+e 



X 



1^ (72(1 — x){l — y) — (Jix{l — y) + M'^y + q'^xy{l — x) + aix{l — y)| 
The z integration is easily done 



2+e 



y 



-l-e 



1 



dx dy — 

1 + e ^2(1 - x) - aix I lM2y + ^23,^(1 _ 3,) + ^^3,(1 _ 



(T2(l — x)(l — + M^y + q^xyil — x) 



(B.15) 
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Appendix B. One- loop scalar integrals 



and, making the change of variable x — > 1 — a; in the second integral, we have 

-l-e 1 







1 -1 

dx dy 



y 



1 + e — (Ti(l — x) 



+ / dxdy- 



a2x{l — y) + JVPy + q^xy{l — x) 
1 



l+e 



l + e aix — (72(1 — X 

We then decompose / as the sum of the two integrals 

/ = CTa) + /2(o-i,cr2) 

which have the symmetry 

/2(cTi,CT2) = /i(cT2,Cri) , 

SO that / becomes 

/ = Ji (0-1,(72) + /i((T2,cri) 
and we have to perform only one integration 

1 1 



(Jix(l — y) + M'^y + q^xyil — x) 



l+e • 



(B.16) 



/•I -1 
h = dxdy- 
Jo 



-l-e 



l + e a2X-ai{l-x) {a2xy+' |i + ^ + + _ 3;)] | 

Using the identity ( [A.8| ) 

-1 1 1 



l+e • 







dx dy 



H (1 - X) 



l + e (T2X — (Ti(l — x) (cr2x)^+'' \0"2X (72 



X y 

-1 

X y 



-l-e 



-X + — + ^xfl - x) 



Mi + _ a;) 



2e 



(M^)-^ / dxdy 



M 



2e 



(72X — ai(l — x) (cr2x)^+^^ 

\ 2e 



+ g^x(l - x) 



-l-e 



1 -y + y 



X 



Ml 1 

(T2 (T2 



+ f-x(l -x) 



(l-y)^^ + (l-y)^^ , (B.17) 



where, in the curly braces, we have added and subtracted the same quantity. If we 
consider now the first two terms in the curly braces and expand them in e, we have 

\ 2e 



f{x,y) 



i-y + y 



X 



Ml A_ 1 
(T2 0'2 



2enog \-y\y 



X 



Mi + ^x(l 

(72 (72 ^ 



-log(l-y)| + 0(e2) . (B.18) 



B.2. Box witk two massive propagators: -62™ 
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Since we are keeping only terms till order (9(e°), f{x,y) gives contribution to the 
integral only if it is multiplied by factors at least of order 1/e. Such factors may arise 



only at singular points, which, for eq. (|B.17|) , are x = and y = 0: in fact f{x,y) is 
multiplied by In the limit y^O, we have 



l-2e„,-l-e I 



x-'-''y-'-'fix,y) = 2ex-'-'^y 



-y + y 



X 



0"2 172 



+ y + (xy^) 



2ex-^'y-' 



Mi + 1 

<T2 <T2 



+ ^-x{l-x) 



(B.19) 



that gives zero contribution in the limit e^O. In this way, f{x,y) gives no contribu- 
tions and can be disregarded. By integrating eq. ( p.l7|) in y, we obtain 



(M2)-^ 1 
1 + e 7 



X / dx 



lo (J2X — cri(l — x) a2X 

{M^y 1 F, 
1 + e e aia2 



1 + 



X dx 
Jo 



O"! + 0-2 



a2X — cri(l — x) X 



M2 



X 



-2e 



0-2 



(B.20) 



1 + 



(M2)-^ 1 F, 



dx 



(Ti + (72 



1 + e e (Ticr2 [-'o a2X — ai{l — x 
1 + elog h + -^x(l - x) J + O (e^ 



l + 2elog^ + 0(e2) 



— I dx \ ) x' 

Jo V 0^2 / 



l-2e 



1 + elogh +^x(l-x)j +0(e2 



(Xi + 0-2 



1 + e e (Ticr2 ["'0 cr2X — cri(l — X 

2 \ JVf2 



1 + 2e log + 2e2 log^ + ( e^ ) 

(T2 (72 ^ ^ 
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Appendix B. One- loop scalar integrals 



± + eJ^'dxx-'-'^\og(^l + ^x{l-x)^+0{e') |, (B.21) 



finite for e — + 



where 



r(l-e)r(l + 2e) 

r(i + e) 



l + —e' + Oe 



(B.22) 



The contribution given by Ai to the total integral I is zero. In fact, using eq. ( [B.16| ), 
we have 



dx 



CTi + (T2 



(T2X — ai{l — X 



-Ai + dx 
Jo 



0"2 + Ui 



and then 



2\-e 



(M^) 

1 + e cri(T2 



-2((Ti + (J2) / dx 
Jo 



(JlX — 0"2(1 — X 

Inp- £2£ 



-Ai = 



Cr2X — (Ji(l — x) 



-/ rfx -log 1 + ^a; 1-x + — + -log + log' . B.23 

Jo X \ j 2e^ e (72 cr2 J 



From eq. ( |B.16|) we have 

1 + e (Ti(T2 



X |-2(cri + (T2) dx 



log 



o-i(l-x) 



a2X — cri(l — x) 



1 1 / M2 M2 

+ ^ + - log + log 

e \ cTi 0-2 



M2 



+ log^ + log^ 

(Tl (72 



(B.24) 



In this way, we succeed in isolating all the divergent contributions, which have ap- 
peared as poles in e. The remaining integrals are finite. For the first integral in 
eq. (p.24 ) we have 



dx 



log 



<Tl(l-x) 



a2X-(7i(l-x) 2(0-1 + 0-2) 



, 2 ^1 , 2 
log h TT 

0-2 



that can be demonstrated using the contour integrals around the poles 0, 1 and 
(o-i + o"2), and the identity 



log 



02 X 

.o"! 1 — a; 



47r i 



(T2 X 
0-1 1 — X 



-log' 



(72 X 

ai 1 — X 



x+ie 



(B.25) 
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The last integral in eq. (|B.24|) is solved with 



/■I dx , 
/ —log 

JO X 



1 + ^(1-^)]^^ 



^ dx ^ ( X \ dx ^ ( X ^ 
— log 1-— + / — log 1- — 



Lio 



A4- 



1, 2-A- 



Lig 



(B.26) 



where we have defined 



Ah 



4M2\ 



2(1±A) 



(B.27) 



The integral / becomes 

1 + e 0-1(72 1 



A , + A_ = 1 . 



log^ ^ + TT^ 



log^ 1 



1 1 / , 

— + - log + log 



0"! 



0-2 



and, according to eq. ( [B.13| ) 



log h log 



(B.28) 



B2m = (47rrr(l + e)e-MM2)"^^U + - flog— + log— 1 



We return to the physical region with the analytic continuation of -82™ ? that is, 
using eq. ( [B.ll| ) 

(B.30) 



logM^ = log (^—m^ + i?]^ = logm^ + in 
1 ( I 1 , 



You can see that 



< A_ < 1/2 



1/2 < A+ < 1 



(B.31) 
(B.32) 

(B.33) 
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and then 

so that the logarithm acquires an imaginary part 



log (l-3^)=-log^ 
The final result is 

B2m = N{e) (m^)"'- 



1 ^+ ■ 

log ZTT . 

A_ 



1 f 1 1 / , 
W + ~ ^ log — + 2 log — log — 



-TT — lo£ 



where 



A_ 



N(e) 



2TTi 



1 , m , m , A+ 
- + log h log h log — 

e (Ti (72 A_ 



167r^ 



(47r)^r(l + e) . 



(B.34) 
(B.35) 



B.3 Box witk tkree massive propagators: B^m 




Fi gure B.2: Diagram representing B^m- 



We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type represented in Fig. [B.2 

d'^l 1 1 1 1 



B 



3m 



[B.36) 



B.3. Box witk tkree massive propagators: B^m 



105 



Performing the same steps done for the calculation of B2m, we have 



Bsm = (47r)^ r(2 + 6) I 



= 4 - 2e 



[B.37) 



where 
D 



ai a2 m? + ai (g — pY + ai + a2 (f + 0.2 0:4, {q — ky 



+ {a2 + as + a4)m - i-r] . 



(B.38) 



Using the definitions of the kinematical invariants of eqs. ( |B.1| )-( [BI6| ) and continuing 
the integral into the unphysical region of negative mass, according to eq. (|B.11|) , we 
have 



D 



«i(«2 + «4)^ + «! as S2 + 02 «3 ? + a2 ^4 S3 — (1 — cti) m + iri 



-m^(l - aiY + 0-2 ai 03 + 02 + S3 a2 + it] 
M^(l - aiY + (72 ai 0:3 + 0:2 ^3 + 02 "4 



The integral (p.37| ) then becomes 



5. 



3m 



167r2 



(47r)^ r(2 + e) e" 



2+e 



;i + e)N(e) e^"^/ 



(B.39) 



[B.40) 



By making the following change of variables 

ai = 1 — X 
a2 = x(l — z) 
0:3 = xyz 

= x(l — y)z , 

and taking into account the Jacobian of the transformation 

d{ai a2 03) 



(B.41) 



d{xy z) 



x^z 



we have 



dx dy dz 



dx dy dz 



2 

X z 



M^x"^ + 0-2(1 — x)xyz + s^x'^il — y)z{l — z) + q^x'^yz(l — z) 

zx'" 



2+e 



M^x + 0-2(1 — x)yz + S3x(l — y)z{l — z) + q^xyz{l — z) 



2+e 
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/■I -1 

/ dx dz 

io 1 + e 0-2(1 



X 



X 



0-2(1 — x) — S3x(l - ^) + q^x{l — z) 

1 1 



[M'^x + 0-2(1 - x)z + q'^xz{l - z) 

1 



l+e 



M'^X + SsXz{l — z) 



— r 

l + e Jo 



X 



0-2(1 - x) + (g2 - S3)a;(l - z) 



X 



X 



-1-26 



[m^x + 0-2(1 - x)^ + q^xz{l - z)\ [M2 + 53^(1 - z 



l+e 



(B.42) 



When x^O, the second integral is divergent: for this reason, we add and subtract 
the asymptotic expression of the divergent function when x goes to zero 



l + e 



Jo 



dx dz 



X 



0-2(1 - x) + (g2 - S2)x{l - z) 



X 



X 



-l-2e 



W^x + 0-2(1 - x)z + q^xz{l - z) M2 + S32;(l - z) 



l+e 



l+e 



L 



+ / dxdz 



X 



-l-2e 



(J2 



M2 + 53^(1 - z) 



l+e 



— / dxdz 



X 



-l-2e 



(J2 



M2 + S3^(l - 



l+e 



(B.43) 



We can make the following considerations: 



1 . the first term is finite when e — > 



2. the sum of the second and the third term is finite when a; — > by construction, 
and then we can again put e = 



3. the integral of the last term is 



It. = 



dx dz — 
0-2 

(M2)-^ /•! 



X 



-l-2e 



2o-2e Jo 



M2 + 53^(1 - Z) 

dz 



l+e 



(B.44) 



2o-2e 



/ 

^0 



1 + ^^(1-^) 
1 



l+e 



dz ; r H / dz ■ 

M2 + S3Z(1-Z) 2(72 



1 rl ^ log 1 + ^^(1-^) 



/ 

^0 



M2 + S3^(l - z) 
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By collecting a global factor (M^) % eq. (|B.43| ) becomes 



(M 



rl 



1 + e \Jo 



dx dz 



1 



(T2(l - x) + (g2 - s^)x{l - z) 



1 



+ 



W^x + ^2(1 - x)z + q^xz{l -z) X |M2 + S-iZ{l - z) 

1 



CT2 X 



+ /" dz 



M2 + S3Z(1 - Z) 

1 



1 /■! 

dz ■ 



2a2eJo M'^ + ssz{l-z) 2a2Jo + S3z{l - z 



l + jl,z{l-z) 



(M 



1 + e \Jo 



dx dz ■ 



- a^z + q^z{\ - z) 



0-2 



+ I dz 



2 0-2 e Jo M2 + 53^(1 -z) 2 0-2 70 
The integration in x gives 



M2 + 53^(1 - ^) M2x + 0-2(1 - + q^xz{l - z 

. (B.45) 



1 , log l + j^z{l-z 
dz 



dz 



1 



M2 + S3Z(1 - Z) 

+ ^2^(1 - z) 



1 + e ^2 M2 + 53^(1 - z) 



log 



0-2-2 



1 1, 

H log 

2e 2 ^ 



and decomposing the denominator 

1 



1 



M2 + S3;z(l - z) S3A' 



where 



and 



we have 



A' = , 1 + 



4M2 
S3 



2\-e 



1 + e 0-2S3A' Jo 



dz 



.2 1 



1 + log^-^log^ 



- log 2; + log(z - A_) + log(A+ -z)-]^ \og{z log(e+ - ^) I , (B.46) 
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where X± are given by eq. ( p.27|) . 

In order to compute these integrals, it is worth noticing that 

< A- < < 1 < e+ < A+ , 

where we have used the fact that > S3 (see eq. (p.3|) ). In addition 

<- 



dz 



21og- 



^4 



/. 



^ , \ogz J. ( ^ 

dz — = L12 — 

z-^± 



i^^log(.-e- 



dz 



dz 



dz 



dz 



dz 



z — 




\og{z 




z — 




log(e4 




z — 


e+ 


log(2 


-A-) 


z — 


e- 


\og{z 


-A-) 


z — 




log(A_ 


,-z) 


z — 





A' 



z-i- 



logA'log^^^ Lis f^rj + 



^log^(e+)-^log^(-e-) 



2log'(-e-)-2log'(e+) 



log A+ log — - log(-A_) log — Li2 

?7„ 



Li-; 



-A, 



log?7+log— ^ + Lia ( — I - Lis 



log^+ log(-A_) -log(-^^)logAH 
log^+ log(-A_) - log(-^_)logAH 



Z — \j 



log?7+log— ^+ Lia ( — ) - Lis (— ^ 
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/■I log(A+ - z) B.p. , / c M / X x , c ^ \ f\ log(C+ - ^) 
/ dz = log -C- log -A_ -logC+logA+ - / dz 

Jo Z — ^+ Jo Z — 

= log(-C_)log(-A_) -log^+logA+ 



1 T- /-^H 

log ?7_ log -; h Ll2 



where 



A+ V 



r/±^-(A'±A), 



and B.P. means integration "By Parts". The integral / then becomes 
7 = - + 21og^ + log— log^ 

_U.(1).L.(1) + L.(it)_u.(z^) 

+ 21ogA+logr;_ -21og(-A_)logr;_ +2Li2 (^^^ -2 Lis 

+ 21og(-e-) log(-A_) - 21oge+log A+ - loge+log(-A_) + log(-e-) log A+ 



+ Li2 — - Li2 I 

\V+J V V+ 



1 1 

+ 2 ^og^ ~ 2 + ^^S?7+ log 

We can use the following algebraic relations 
C+ 1 + A' l + A' + A-A A 



£^ A^ 



(B.47) 



77+ A' + A 



A' + A 



+ 1 



V+ V+ 



and some properties of the dilogarithm function 
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Li2 ( — I - Li2 ( - — 



-2 Li, 



log log h — 

V+ V+ 6 



Lio 



Li2 



2 Lis 



log' 



to obtain 
/ 



2\-e 



1 + e (JsSsA' [ e 
■ 1 



1 J 1 1 
; -log- 



^4 



2 log h log log — — 

0^2 S3 / ^+ 



+ 2 Lis 



2 Lis 



2 Lis 



'A. 



2 Lis 



A. 



TT 



'A_^ 
-A_ 



-2L„^-^J+2U.^-g,J---21og,_log 

+ 2 log 7]+ log^^ + log2^+ -21ogA' log(-^_) + log2 A' 



(B.48) 



To return to the physical region, we must continue analytically the solution, using 
eqs. ( [B30|) -(p3^ ) and 



1 / / 4m2 \ 1 
e. = - 1±^1-^ ±.,= -(l±A')±^,. 



(B.49) 



You can see that 



0<A_,e-<l/2, 



l/2<A+,e+<l 



(B.50) 



and then 
1 



> 1 



+ irj 



1 - A' 
1 + A' 

1 - A' 



< 



1 - A' 



+ 



1 + A' 1 + A' Vl- A' 1 + A 



IT] 



-— + ir] 



>o 
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In a similar way, you can demonstrate that 



and 
A_ 



1 + A 



ri+ A' + A 

1 + A 



A' + A 

>o 



> 1 



1 + ir] 



1 + A + ^^tT^ 



A' + A + .,(^i + |i.) 
1 



1 + 2!A^ 
1 + A A' + A 



tr) 



<o 



?7_ T]- 



?7_ — > r]_ + if] . 

Using eqs. (|A.6| ) and (|A.7|) we can write B^m as 

1 f 1 , ^- 



Bsm = N{e) {m') 



2\-e 



0-2 S3 A' I e ^+ 



0-2 



log^ 





log^ 


g — 

S3 ; 


v+\ 


1 2 

-log — 

V+ 




-2Li2 


(-!)■ 







7r2 
2 



+ 2 log ( I log ^ + log' ^+ - 2 log A' log ^_ + log' A' 



+ in 



1 a' e 

- + 2 log — + 4 log 7]+ - 2 log A' + 2 log — 

e ^2 ^4 



(B.51) 



B.4 First triangle witk two massive propagators: 

We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type illustrated in Fig. |B.3 



'™ J {2ttY P {I - j/y -m? {l+p + kf -m? ' ^ ' ' 
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Figure B.3: Diagram representing Tg, 



2m- 



Using eqs. (|A.15| ) and (|A.16|) we can write 

i f [da] 



TL = (4vr)^ r(l + e)(-l) / 1^ , d = A - 2e 



where, with the help of eqs. (|B.lD -( p^ , 



D 



[B.53) 



«! a2 + ai a^^p + kY + as + [0.2 + 0:3) — ir] 



«! + «! as Si + ^2 0^3 — (02 + 03) + irj 



-m^(l — ai)^ + 0"! «! 03 + g^a2 03 + 



(B.54) 



Continuing the integral into the unphysical region, according to eq. (|B.11D , we have 



D 



M^(l — + CTi «! as + q^a2 03 



The integral ( |B.53| ) then becomes 



rpq 

-'-2m 



167r2 

The following change of variables 



{AtiY r(l + e) 



[da] 



-D] 



l+e 



N(e) e'"^J 



(B.55) 



(B.56) 



ai = 1 — X 
a2 = x(l - y) 
a3 = xy , 
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together with the Jacobian of the transformation 



d{xy) 



X , 



gives 



dx dy 



dx dy 



X 



M^x^ + crix(l — x)y + q'^x'^y{l — y) 



M^x + 0-1(1 — x)y + q^xyil — y) 



l+e ■ 



(B.57) 



The integral / is finite when e — > 0, as can be seen also by inspecting eq. ( [B.52D , that 
has no soft or coUinear divergences. For this reason we put e = and integrate in x 
to obtain 

M2 + g2^(l-y) 



dy- 



1 



log 



/o M"^ - oiy + q^yil - y) aiy 
By decomposing the denominator in the following way 
1 1 



(B.58) 



M2 - aiy + g2^(l - y) _g2 



1 _ ^) _ 



q\^4 + W\y-p+ y-p- 



(B.59) 



where 



«! = 1 — 



0-1 



I Am? 



<ai<l 



we have 
/ = 



dy- 



tti ± Waf + 



log 



4M2 



° q'^|<^+^\y- p+ y-p-, 



-q' {y-X+){y-X-) 
0^1 y 



It is easy to see that 
and / becomes 

1 



p_ < A_ < < p+ < 1 < A 



+ ) 



(B.60) 



(B.61) 



dy 



q'^o[+^[-^o -y-p+ 



1 K-y , . y-X- 

log h log ■ 



K - P+ p+- A, 



1 y 

log — 

p+ 



dy 

y- 



, K-y , , y- X- , y 

log h log ^ log 



A+ -p- 



A„-p_ 



-P- 



(B.62) 
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We summarize here the results of the integration in y 

1 ^ ^ K-y T- f p+ \ J- ( ^- P+ 



Jo y-p+ A+-P+ \ A+-p_y \X+-p+^ 



I dy log — — = Li 

Jo y - p-i- p-i- - A_ 



dy log^ = -Li2 1 + — 

y-p+ p+ \ p+J 6 



Cy^log^ = Li. ( ^) - Li. 

y-p- K-p- \K-p-J \K-P- 



1 1 , y-X- 1 - A_ \ , . / A_ 

dy log = Li2 - Li2 

y-p- X--P- V ^--P~ \X--P- 



1 - p_ , 1 - A_ , -p_ , -A_ 
+ log log log log ■ 



A_ — p- A_ — p_ A- — p- A_ — p- 



I dy — ^- — log = Lis ( — ) - log(-p-) log 

io y- P- -P- \P- 



y-p- -p- \p-j -p- 

We can then write 

1 f_!l!+Li.fl-i-ULi/ 



- log(-p_) log ^—^ + Li2 { \ - Li2 ^ ^ 



-p- \K-p-J V K-p-) 

T- f ^- \ T- f 1 - A_ \ , 1 - P- , 1 - A 
+ Li2 - Li2 -T - log log 



X- — p- J \ X- — p- J X- — p- X- — p- 

-P- 1 -A- 
A_ — p- X- — p- 



+ log-^^log-^— I . (B.63) 
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According to eqs. (|R30D -( F32|) and noticing that 

P± 



If , 4m2\ 



(B.64) 



we see that the ordered sequence ( [B.61| ) becomes 

< A_ < p_ < p+ < A+ < 1 . 

This fact, together with the values of the signs of the imaginary parts, helps us to 
continue the solution into the physical region. In fact 

1- p+ 



1 


- P+ 


^+ 


- P+ 


1 


- P- 


A+ 


- P- 


1 


- A_ 


P- 


- A_ 




P+ 


P+ 


- A_ 


1 


> 1 - 




P- 




1 


- P- 


A_ 


- P- 


1 


- A_ 



> 1 



> 1 



> 1 



> 1 



K-P+ 
i-p- 

X+-P- 
1 - A_ 

p+ 

p+ - 



17] 



IT] 



17] 



A_-p_ 



< 



< 



< 



^ 1-p- 

A_ — p_ 
1 - A_ 

" a7^ 



— %7] 



17] 



-P^ -P- 

The final expression is 



TL = iv(6)K) 



2\-t 



17] 



iL^ Lijl-i. 



Li2 



P+ 



A+ - Ph 
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+ Li2 
- Li2 
+ Li2 



. (K-P+ 



1 


- P+ 


A+ 


- P- 




- P- 


P- 


- A_ 



1- A_ 



- Lio 



- Lio 



' P+ - 
. P+ . 

P- 



— Lio 



1 
2 



\+-p- 
\+- p+ 1 



+ Li. 



'i^p±: 

/ A_ 



- Li2 (p-) 



1-P+ 



+ ^ log^ - ; log^ - 1 log^ p_ 



1 ~ P- 1 1 2 A+ — P- 1 

- log p_ log - - log2 7 + - 

p_ 2 1 — P- 2 



-log 



+ i7r 



A_ 



1 1 A_-l , 
— log T h log 



A_ 



log 



1- A- 
-A_ 



A_ — p- A_ — p_ 



+ 



TT 



O 1 A+ - p_ 1 - Ph 

2 log T h log ■ 



A+ -p+ 



l-p_ 



(B.65) 



where we have factorized out the term (m^) ^, that gives zero contribution in the 
limit 0. 



B.5 SeconJ triangle witk two massive propagators: 




Figure B.4: Diagram representing T2m- 
We are referring, as far as the denominator structure is concerned, to a Feynman 



B.5. Second triangle witk two massive propagators: 
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diagram of the type depicted in Fig. 

dH 1 1 



2m 



1 



i f [da] 



167r2 



(47rrr(l + e)(-l) / 



d = A-2e 



(B.66) 



where we have performed the usual steps, aheady done for the previous integrations, 
and where 



D 



ai az + oti as + k)^ 

2 



+ (ai + 0:2) m - %r] 



(«! + a2 — «i «3) fn + 0:2 as si + ir] 
vr?{\ - as)^ + cTi a2 as + ir] 



(B.67) 



Continuing the integral into the unphysical region, according to eq. ( |B.11| ), we have 



D 



M^(l - as)^ + (Ti a2 as 



(B.68) 



The integral (|B.66| ) becomes 



2m 



167r2 



(47r)^ r(l + e) e^ 



\dQ\ 



N(e) 



(B.69) 



and, with the following change of variables 



ai = x(l - y) 
Ci2 = xy 
as = 1 — X , 



with Jacobian 



d{ai a2) 
dixy) 



X 



we get 



dx dy ■ 



X 



X 



dx dy 1 , 

[M^x + ai{l - x)y]^^' 



10 [M^x^ + aix{l - x)y]^^" 

This integral has no soft or coUinear divergences, so we can put e = 



(B.70) 



dx dy 



1 



M'^x + (ri{l — x)y Jo ai{l — x 



dx 



1 



log 



M'^x + ai{l -x) 



(B.71) 
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With the change of variable z — 1 — x we obtain 



ai Jo z 



log(l-z)| 



1 



TT 



(B.72) 



The analytic continuation into the physical region (M^ = — + irj) is straightfor- 
ward. In fact we have 

which gives 

= N{e) (m^) -'^[u,(-^)+log^log(l + ^)-inlog(l+ 



(B.73) 

where we have factorized out the term (m^)~^, that gives zero contribution in the 
limit e— > 0. 



B.6 Triangle witk one massive propagator: Tj 



Im 




Figure B.5: Diagram representing Ti^. 
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We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type represented in Fig. B]5 

d'^l 1 1 1 



Tlr 



{2ny P {l + kf {l + p + kf-m^ 
i , , ..I' [da] 



167r2 

where, with the definition ( |B.11| ) 



(4vrrr(l + e)(-l) / 



d = A-2e 



D 



di as (p + kY + a2 as + as — irj 



«! as si + a2 as — as + irj 



-m^al + o"i ai as + ir] 



M'^al + 0"! ai as 



The integral ( |B.74| ) becomes 



(AtcY r(l + e) e' 



[da] 



167r2 ' - V- ■ - - J 

so that, with the following change of variables 

ai = x(l - y) 
a2 = I — X 
a3 = xy , 

and the Jacobian 



N(e) e^'^J 



d{ai a2) 
dixy) 



X , 



we obtain 



X 



dx dv 

[MH^y^ + aix^il - y)y]^^' 



X 



-l-2e 



(B.74) 



(B.75) 



[B.76) 



dx dy 1 . . 



Integrating in x 



'2eIo ^i+e 



-l-e 



l+e ' 



(B.77) 
(B.78) 



and using the identity ( |A.8| ), we have 



2e 



2e 



dyy 



-l-e 



1 2e 



1 /. 1 fM^V" r 1 



— (M^) — 



CTi \ (Ji 



--2eLi2 (l-^) 



(B.79) 
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Expanding ("^^ in powers of e 



2e 



1 



M2 



2? + 7^°^^+^H^-]f^J+^°^V 



0-1 



M2 



and continuing analytically the solution 



Ti^ = iV(e) (m 



m 



m 



+ - log h log — log 1 + 



m 



Li. ( 



-7r + i TT 
6 



- + log(l + ^ +21og — 
m"'/ (Ti 



(B.80) 



(B.81) 



B.7 First ckeck triangle: Tf^^ 




Figure B.6: Diagram representing T. 



q—k 
2m ■ 



This appendix and the following one are deserved to compute two diagrams useful 
for a partial check of the calculation. The denominator structure is illustrated in 

Fig. El 

dH I I 1 



^^"^ / {2tiY P {I - V'Y -m? {1 + pY - m? 

- ' (4.).r(l+e)(-l)/l^"l 



d = 4 - 2e 



(B.82) 



B.7. First ckeck triangle: T2m^ 



where 



D 



ai a2 + ai 0^3 + ^2 as {q — + (0^2 + as) — 
-m^(l - aif + S3 0:2 03 + it] 
M^(l - ai)^ + 53020:3 • 



The integral ( B.82 ) becomes 



q-k ^ 



2m 



167r2 



(47r)^ r(l + e) e^ 



[da] 



-D] 



N(e) e'^'I 



and, by changing the variables of integration, 



ai = 1 — X 

02 = 

03 = x(l — z) 



with Jacobian 



d{ai 02) 
dixy) 



X 



we have 



[M2x2 + X^Z{1 - Z)ss] 



1+e 



dx dz 



X 



-l-2e 



[M2 + z(l - z)ss] 



1+e 



Performing the x integration, we obtain 



M2 



M2 



which is equal to the integral computed in eq. ( p.44|) , so that we have 



rpq—k 



Nie) e^- (m'Y'J-S^ log ^ + log — log 



+ -log^e+-logA'log^+ Li2 ( 1 - Lis I g; 



i- 1, 2 
log 

'^4 



The analytic continuation of this solution, using 



122 



Appendix B. One- loop scalar integrals 



IS 



pq-k 
'2m 



N{e) (m^y 



— r7 i - log ^ + log — log ^ - log ^- + o log ^4 



log A' log — + Li2 



2 ^+ 



«4 



-TT 



+ i TT 



- + log 2 log A' 

e S3 



(B.88) 



B.8 Second ckeck triangle: T^m 




Figure B.7: Diagram representing T^m- 



This appendix (as the previous one) is deserved to calculate a diagram used for a 



partial check of the calculation. The denominator structure is depicted in Fig. B.7 



- 3m 



{2'nY {I — p'Y — irfi {I + p — qY — m? [l + pY — m? 

[da] 



167r2 



(4vrrr(l + e)(-l) / 



d = 4 - 2e 



(B.89) 



B.8. Second check triangle: T^rn 



where 



D 



ai a2 (p + p) + ^2 as I? +m - ir] 



+ S3 ai as + q^a2 as 



The integral (p.89| ) becomes 



T- 



Sm 



167r2 



{A-kY r(l + e) e^ 



[da] 



-D] 



N(e) e^"^/ 



and, with the following change of variables, 



«! = xy 

012 = x{l - y) 

as = 1 — a; , 



we have a Jacobian 



so that 



d{ai 0:2) 
d{xy) 



dx dy ■ 



X 



X 



ll+e 



[ss x(l — x)y + q^x{l — x){l — y) + M^] 
This diagram has no soft or collinear divergences, and then we can put e 
can integrate in y to obtain 

1 



1 1 M"^ + xil - x)s^ 
dx log ■ 



s?j — JQ \ — X + a;(l — x)q'^ 

By splitting log (a/6) = log a — log 6 and using 



1 dx 



X 



log 



1 + x(l — x) 



S3 

M2 



1 dx 



X 



log 1 



^4 



+ 



1 dx 



X 



log 1 



we have 



Ts^ = iV(e) i ^ 



2 ss - g2 
The analytic continuation gives 



-Li2 — - Li, — =-log' 1 



3m 



'v(.)('"r1;7^{v(i-i)-v(^ 



2i7r 



log — - 1 - log — - 1 
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B.9 Partial ckeck of tke virtual integrals 
B.9.1 Ckeck of 

A partial check of the correctness of the above formulae can be performed in the 
following way. We consider first a check of i?2m- To this purpose, we define I as 



d'^l 1 + A{1- kf + B[{l+p-qf - m2] +C[{1+ pf 



m 



{2nY P(l - kf [{l+p- qf - m2] [(/ + pf - w?] 

1 + g[g2 -2p-q]+2l- [-Ak + B{j) - q) + Cp] + f{A + B + C) 
{2tiY P{1 - kf [{l + p- qY - m2] [(/ + pf - 

If we impose that J has no soft and collinear divergences, we have 

\ + B\q^ -2p-q\ =0 
k- \-Ak + B{p-q) + Cp\ = . 
The solution of this system is 

0-2 

with the term A undefined. We can then rewrite I as 

d'l 2/ ■ [-Ak - - g) - J-p] + /2 _ J_ _ J_ 



(B.96) 



{27rY P{1 - kf [{l+p- g)2 - m?] [{I + pf - m?] 



B2m + AT^ Tim 



(72 



(B.97) 



where the second line is obtained directly from eq. ([B.96| ), and the primed quantities 
are obtained with the substitution p p', that is ai ^ 02- The integral J is now 
convergent and the cancellation of the divergent part in the right-hand side can be 
checked directly (both in the real part and in the absorptive one). In fact: 
cancellation of the terms proportional to — 



1 1 



1 1 



a\a2 02 2cri o\ 2o2 



B2„ 



T 

-'in 
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cancellation of the terms proportional to 



m 



m 



1 1 



m 



1 1 



ra 



log V log h 27ri > < log h iTT > < log V iir 



o-i 



(72 (Ji 



ai (72 



0-2 



B2r, 



T 



cancellation of the finite terms: the finite integral / of eq. ( [B.QTp can be written 

as 

1 1 

0"! 0-2. 



I + [A- — - — ]T.'^ 



2m 



dH 



21- 



(B.98) 



{2ttY P{1 - ky [{l+p- qY - m2] [(/ + pf - w?] ' 

In order to compute /, we introduce the Feynman parameters (see eq. ( |A.1CI| )) to 
rewrite the denominator in the following form 

4 



den = ^ait^ + ai{l — kf' + + P ^ q) 



,2 2 

— q) — m 



+ a2 



{l + pf - } 



(ai + 02 + 03 + 04) f' + 21- [-k 0:4 + (p - g) 03 + p 02] + • • • , (B.99) 



where the are illustrated in Fig. [B.l| . With the following change of variable 

/' = / + [-fca4 + (p-g) 03+^02] , (B.lOO) 
we have a denominator that contains only /'^, so that I becomes 



r(4) 
2 



{2Tiy 



[da] 



I' -[-kai + {p~q)a^ + pa2] ■ - Ak ~ Mp - q) - j-p 



r + •••]' 

-2 [-/ca4 + {p-q)a^+p 02] • \-Ak - ^{p - q) - ^P 



= N{e){l + e)e- [da^ r 2, ^ ,2 ^ ^2 ^ v ■ 

J [-m^(a;2 + a^Y + a2 ai 0:3 + g 0:2 03 + cti 0:2 0:4] 

where we have put e = in the integral because it is finite by construction. Expanding 
the dot-product we get 

/ = A^(e)(l + e)e [da\-^ ■ — ^ ■ 

= 7V(e)(l + e)e^'^^(-2)[i7/3 + K/2] , (B.lOl) 
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where 



CTi (72 / 2(7 



'1 



if = A - ^ + 



2 VcTi (T2/ 2(72 

In the Euchdean region used to compute the one-loop integrals and with the change 
of variables of eq. (p3.14| ), we have 

U2 31 = / dy dz ^ ^ ^^^^ 7) . 

^'^ h [W^y + a2z{l - x){l - y) + q^xy{l - x) + aix{l - y)(l - z)Y 

Performing the z and y integration, we obtain 

/•If 1 , M"^ + q^x{l-x) 

/{2,3} = Ldx\-——j— . ^J og 



_l_ q2x(^l — x) — aix aix 

1 ^^^M^ + q^xil - x)\ {x, (1-x)} 



_l_ q'ix{l — x) — (72(1 — x) (72(1 — x) J (72(1 — x) — (Jix 

This integral can be expressed in terms of dilogarithm and logarithm functions; 
nevertheless we have performed a numerical integration in order to check the iden- 

tity i wm - 



B.9.2 CkeckofB 



3m 



With a reasoning similar to the previous one, we can check B^rn- We introduce the 
following integral 



dH l + A [{I + pf - m2] + 5 [(/ + p - qf -m^] + C [(/ - p' 



m 



i2n) 



P [(/ + py - m2] [{l+p- qf - m?] [{I - p'f 



dH l + B [g2 -2p-q] + 2l- [Ap + B{p - q) - Cp'] + f{A + B + C) 



(27r)' 



P [{I + pY - m2] [{l+p- qf - nP] [{I - p'f 



This integral has only soft singularities, but no collinear ones, which can be removed 
if we require that 



1 + 5 g - 2p ■ g 







that is 



B 



1 

0"2 
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with A and C undefined. We can rewrite I as 



21- 



Ap + -Mp-q)- Cp'] + l\A -^ + C) 



{2nY P [{I + pf - m?] [{l + p- qf - m?] [(/ - p'f 



Bsm + AT2^ T. 

(72 



q—k 



(B.102) 



where the primed quantities are obtained with the substitution p ^ p', that is 
o"! ^ (T2- The integral / is now convergent and the cancellation of the divergent 
part in the right-hand side can be checked directly (both in the real part and in the 
absorptive one). In fact: 



cancellation of the terms proportional to — 

e 



1 



(T2 S3A' 



log ^ +iiT 



rpq — k 
-'2m 



cancellation of the finite terms: the finite integral / of eq. ( [B.102D can be written 
as 

2l-[Ap~^^ip-q)-Cp'] 

27rY P [{I + p)2 - m2] [{l+p- g)2 - m^] [{I - p'f - rrP] ' ^ ' ' 
With the same procedure used to check i?2m in the previous appendix, we can write 

-2 [a2P + az{p-q)- a^p'] ■ \Ap - ^^{p - q) - Cp' 



[-m^(l - ai)^ + (72 «! tta + q^a2 as + S3 ^2 ^4] 
where we have put e = in the integral because it is now finite, and the are 
illustrated in Fig. |B.2| . Expanding the dot-product, we have 

—2 (02 J + ctsL + a^R) 



N{e) (1 



e e^ 



[da] 



[-m2(l - aiy + (72 ai a^ + q'^a2 a^ + S3 a2 a^y 



Nie) (1 + e) e'"^(-2) [Jh + L/3 + Rh] 



(B.104) 



where 



J=(A- — + C]m' + 

\ 0-2 

L = (a - — + C]m^ 

\ (72 

R = (a - — + C]m^ 

(72 



1 g2 



(72 



(72 



(72 



A 



- CTi - (72 



(T2 



c 



(7i - (72 



A 



1 

2 
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In the Euclidean region used to compute the one-loop integrals and with the change 
of variables of eq. ([B.41| ), we have 

r f^j J J xz{{l- z),yz,{l-y)z} 

-'{2 3 41 = / dxayaz- — — — — - . 

Jo [M^x + (72(1 — x)yz + S3x(l — y)z[l — z) + q^xyz[l — z)Y 



From the identity 



^ dx ^ - 
[ax + 6]^ a? 



. a + b b 
log + r - 1 



(B.105) 



b a + b 

we can integrate in x, where a and b are given by 

a = M^ - a2yz + 53(1 - y)z{l - z) + q^yz{l - z) 
b = (72yz . 

At this stage of the integration, /{2,3,4} are two-variable integrals, which can be ex- 
pressed in terms of dilogarithm and logarithm functions; nevertheless we have per- 
formed a numerical integration in order to check the identity (p.l02|) . 



B.IO Self energy witk one massive propagator: S 



Im 



^3 



1 



h+1 



Figure B.8: Diagram representing S 



Im ■ 



We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type illustrated in Fig. |B.8 



'S'lr 



d'^l 1 



1 



where /i^ > and 



D = —(X\ a2 h'^ + 0.2 ^ — if] ■ 



d = A-2e, 

(B.106) 

(B.107) 



B.IO. Self energy witk one massive propagator: Sim 
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We assume h"^ < and, at the end, we return to the physical region according to the 
ir) prescription. We then define 



so that 



D — H'^ai a2 + m^a2 , 
and, with the following change of variables, 

ai — 1 — X 
a2 — X , 



(B.108) 



we obtain (the Jacobian is equal to 1) 



^0 



167r2 



e 

1 



m?x + H'^x{l — x) 

x~'^ 



(B.109) 



The integral is finite, and then, expanding in e, we get 

Sim — N{e) (jn^^ - ll — e j dx \ogx — e j dx log 



= iV(.)(m^)-|l + 2-(l + ^Vogfl + ^ 



(B.llO) 



For }i? > rri^, we have to continue our solution according to 

U + I ' log U ~ — ~ 1 = log I — - 1 ) - ZTT 

and we obtain 

Sim = N{e) ( 



m 



1 2 _ ( 1 _ ] lo ( — 



1 + iTT 1 



m 



(B.lll) 



Two particular expressions are interesting: Sim and its derivate computed at h? = m?. 
The first expression is easily done 



'lm|/i2^^2 — Simm 



N{e) (m^) ' + 2 ) . 



(B.112) 
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For the derivate, we have to start from eq. (|B.109| ) 



dS 



Im 



h'^=m? 



N{e) (m^ 



1 



dx 



t^x ''{1 — x) 



N{e) (m 



e Jo 




1 














m? 


1 












/ dx 


x~ 




m? J 


'0 




















v"2^ 






)■ 



1+e 



l~x) 



(B.113) 



B.ll Self energy witk two massive propagators: S2m 




S. 



Figure B.9: Diagram representing S2m- 
We are referring, as far as the denominator structure is concerned, to a Feynman 

[da] 



diagram of the type represented in Fig. p.S 
d'^l 1 1 



2m 



{27rY P -m? {I- hf - m? W 
where h? > 4m^ and 



{Any r(e) 



= 4 - 2e , 
(B.114) 



D = -h'^aia2 + m'^ - if] . (B.115) 
Again we assume < and define 

= -(h^ + i^) > . 

The Jacobian of the following change of variables is equal to 1, 



B.ll. Self energy witk two massive propagators: S2m 
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so that 



5. 



2m 



167r2 



(4vr) 



r(i + 6) 



dx ■ 



e Jo [m^ + H'^x{l — x)Y 
1 



(B.116) 



l + ^xil-x] 



The integral is finite, and we can expand in e 



2m 



N{e) 'ijl 
N{e) M"'!- 



e dx lo£ 



1 H ^x(l - xj 



+ 2 + (r+-r_) log 
e r_ 



(B.117) 



where we have used the fact that 



^2 0-2 

n n 

1 H rx(l - x) = — T(r+ - - r_) 



with 



4m2 



If we are interested in the value of ]^ such that /i^ > 4m^, we continue analytically 
our solution 



r± = - l±Wl- 



4m^ 
7^ 



and we obtain 



52„^ = iV(e) (m^) ' |i + 2 + (r+ - r^) log ^ + 



2 TT r_i — r_ 



(B.118) 



(B.119) 



Two particular cases are interesting: the one with m = and the one with ^ 0. 
For both cases, we start from eq. (p.ll6|) . If m = we have 



S2mL=o = N(e) {h^) ' - /' dxx-'{l - x)"^ = N(e) [h' 



2 . (B.120) 



When h"^ ^ 0, we can expand eq. ( p.ll6| ) in e and h"^ 



5. 



1 



dx 



1 — e log 1 X (1 — X 



N(e) (m^) 

^ ' e Jo 

N{e) [m^y^l^dx 1 + e f^x(l-x)) + O (/i^) 



O 



'1 Ih^' 
e 6 



(B.121) 
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In addition 



(5: 



2m 



h 



2m 



J {2nY (/2 - m2) [(/ - /i)2 - m2] 2 
that can be demonstrated by contracting both sides with h^, and expanding 

h-l = -^-{[il-h) 
In the same way, we can write 



\2 2 

^ m 



;2 2 

/ — m 



{S: 



2m 



with 



A 



B 



d-1 
1 

d- 1 



-l\L,„,+ {-h^-m^]S; 



2m 



(B.122) 



;B.123) 



(B.124) 



(B.125) 



that can be obtained by contracting eq. (p.l24]) with h'^ and with g'^'^. The expression 
of Lim is given in eq. (|B.132| ). 



B.12 Vacuum polarization witk two massive propagators: 




Figure B.IO: Diagram representing P2m- 



This appendix and the following one are to be regarded as particular cases of 
Appendix |B.11| . Nevertheless, we summarize here the results. The structure of the 



B.13. Vacuum polarization witk no massive propagator: Pom 
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denominators of this Feynman diagram is depicted in Fig. B.IC 
dH 1 1 



2m 



with 



{27vY P - m2 (/ - ky - m? 167r 



/c^ = , D = («! + — ir] = 



The result is simply given by 



IGvr^ e Jo [rn?) 

that agrees with eq. (p.l21| ) for /i^ = 0. 



d = A-2e, 
(B.126) 

(B.127) 



(B.128) 



B.13 Vacuum polarization witk no massive propagator: Pom 



1-k 



Figure B.ll: Diagram representing P( 



Orra • 



We are referring, as far as the denominator structure is concerned, to a Feynman 
diagram of the type represented in Fig. [B.ll 



Om 



rf'^/ 1 



1 



ci = 4 - 2e . 



Since no invariant scalars are present, because fc^ = 0, we have 



(B.129) 



-Pom — . 



(B.130) 
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B.14 Massive loop: Li^ 




Figure B.12: Diagram representing L 



Im ■ 



The single loop is illustrated in Fig. |B.12 

J _ [ d'^l 1 

and can be easily integarted 
L 



d = 4:-2e, 



i.. = -7T^(4vrr^^^ = iV(e) H 



(An) 



m" I — hi 



(B.131) 



(B.132) 



A particular case is the massless loop, where, again, we do not have invariants to 
build the solution 



dH 1 



777=0 



(B.133) 
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Renormalization 



C.l Radiative corrections to external keavy-quark knes 

1 

# % 

^ i k k j 

Figure C.l: Diagram representing the quark self-energy T,ij{h). 

In this appendix we describe how to treat loop corrections to external heavy- 
quark hues and how to compute the mass counterterm. 

The one loop correction to a quark propagator reads 

where the sum over repeated indexes is meant and /i is the mass parameter of the 
dimensional regularization, introduced in order to keep gs dimensionless. With a little 
algebra we have 



{2nY P[{h + iy-m^] 



{2%^ P [{h + If - m? 
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+ {2-d) J 



In the last integral we can replace 

h„, 1 



l,~^l-h-f,=l[[{h + l) 



2 2 

— m 



(C.2) 



(C.3) 



and obtain 



m 



-1+e- (1-e)— 1 |(+(4-2e)m 



S 



Im 



+ ^(l-6)Li^ , 



(C.4) 



where we have used the definitions given by eqs. ( [B.106| ), ( [B.132| ) and ( |B.133| ) and we 
have specified (i = 4 — 2e. We are interested in the expansion of around ^ = m 



Using the identity 



+ 0[{li-m) 



(C.5) 



so that 



2m 



d 



(C.6) 



(C.7) 



and eqs. ( |B.112[ ) and ( |B.113| ), with simple algebraic passages, we can rewrite ( |C5| ) as 

■ ,2\^ .rQ n r 3 



m 



2 "'J 



5 + 4 

e 



m+ (J/i — m) 



4 

e 



+ 



m] 



= 5. 



—i 6m — i — m) zq + O (J^J^ 



— m] 



where 



Sm = -^g^^CFN{e){^ 







[5 + 4] 









m 





e 


[^4] 









(C.9) 
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The full quark propagator at first order reads 



+ - T,ki[h) - — '— + 0[a^j 



m Iji — m 



i 5 



i 5. 



ik 



-i5m) - — + 

f — m f — m 



(C.IO) 



If we want that the pole of the propagator is not displaced by radiative corrections, so 
that m corresponds to the pole mass definition, we have to add a mass counterterm 
to cancel the second term of the above expression. For this reason, we define the 
Feynman rule for the mass counterterm as the insertion, in the fermion propagator, 
of the vertex —i mc, where 



rric = -5m = -glCp [ 



5 + 4 

e 



m 



(C.ll) 



and is defined in eq. ( ^.38|) , that is 



N, = -iN(e) 



(47r)^r(l + e) 



In this way, slightly off-shell, the quark propagator behaves like 



— m 



with 



l + ZQ = l~glCFN, 



+ 4 



(C.12) 



(C.13) 



C.2 Radiative corrections to external gluon lines 



Gluon, light-fermion and ghost self-energy corrections to external gluon lines vanish 
in dimensional regularization, so we need to consider only the correction coming from 



the heavy-fiavour loop, represented in Fig. |C.2| . This contribution is given by 



n 



dH 

(2^ 



Tr 



/- 



m 



(C.14) 
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I 




Figure C.2: Diagram representing the gluon self-energy contribution coming 
from a massive-quark loop. 



n 



ab 



where the minus sign takes care of the fermion loop. With simple algebra 



d'^l Tr 


7^ (/ + m) 7^ [(/ - ^) + m] 


{27rY {12 _ ^2) 


(/ — Kf — m? 





{2'kY 



— m?) {I — h) — m? 



, (C.15) 



where Tr [t^t''j = Tp 5^^, Tp = 1/2. Using eqs. (|BTT|), ( [Bl22D and (|RT24| ), we can 
write 



2..2enn xab u2 



where 



^,,2x C?-2 1 , 



d-1 

In = 4 — 2e dimensions, we eet 



d-1 



'^-1 + 2^ 
2 



2m ■ 



(C.16) 



(C.17) 



n 



ab 



-gtfx'^TpS^'h'lg,, ^2 



e / 2 

1 - - + 2— 1 + -e 
3 \ 3 



3j 



Im 



2m 



(C.18) 



We are interested in the form of the propagator for small gluon virtuality: using 
eqs. ( IB.121I) and (|B.132| ) we obtain 



(C.19) 



C.3. Ckarge renormalization 
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It is worth noticing that, as previously stated, for small gluon virtuality and for the 
massless-quark loop, n(/i^) is zero, because Li^ = and 5*2™— *-Pom = 0. The full 
gluon propagator at first order reads 

= ^^-^'^^^^['^^) 37 + ^^-^"^ 

The part of the full propagator proportional to h^h,^ gives no contribution, since it is 
contracted with a conserved current. For this reason we can write, for the slightly 
off-shell gluon propagator 



G,{h) = ^^Z, , (C.21) 



with 



Zg = l + z, = l-N,TFgl(^\ ^. (C.22) 



C.3 Ckarge renormalization 



We carry out the charge renormalization in the mixed scheme of Ref. [^], in which 
the nif light flavours are subtracted in the MS scheme, while the heavy-flavour loop 
is subtracted at zero momentum. 

In this scheme, the heavy flavour decouples at low energies. The prescription for 
charge renormalization is 

tts *■ 'Zeis Q^s = ^vertex ^qulrk ^gluon '^s (C.23) 



^vertex quark gluon 

1 + = 1 + glN, ^ ^ 



(C.24) 



where Ca = A^^c = 3 for an SU(3) gauge theory. 

We try to justify the Tp structure of this expression without entering in the 
details of the calculation (albeit straightforward). Three renormalization constants 
determine the charge renormalization Za/. ^vertex, that takes account of the two 
Feynman diagrams describing quark-quark-gluon interaction; Zquark? which contains 
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the contribution coming from the quark self-energy and ^giuon, that includes the 
effects of vacuum polarization. 

The contributions proportional to Tp come only from the gluon self-energy part 
that takes account of the corrections due to quark loops. In fact: 

- the vertex correction contains the colour factors Ca and Cp, and, in the MS 
scheme, gives 



Z vertex — 1 ~" 5's^e~ 



1 \ 3 ■ 
Cp--Ca]+-Ca 



l-glN^{Cp + CA) (C.25) 



the radiative corrections to quark propagator brings (see eq. (|C.4|) ) 



^quark= l-^^X^C*^ (C.26) 

the gluon- and ghost-loop corrections to gluon propagator contribute to the 
renormalization constant with 

^g^on = 1 + 9lN^ {^Ca + ^Ca)=1 + <7.^iV,| Ca (C.27) 

the contribution to the gluon self-energy coming from a massless fermion loop 
is given by (see eq. (|C.18|) ) 



where we have used eq. (P.132| ) and the value of 5*2^ given by eq. (|B.120| ). In 



this way, in the MS scheme, the massless loop contributes to the renormalization 
constant of the gluon propagator by an amount 

Zl,^,^ = l-N,Tpn,,g',^^, (C.29) 

where we have summed up the contributions of riu light quarks. 

if we subtract the heavy-flavour loop at zero momentum, we perform the same 
steps we have done to obtain the expression of Zg (see eq. ( |U.22| )) and we get 

zSnon = Zg = l- N^Tp gl 1 . (C.30) 
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Combining eqs. (102^ ), ( 10291 ) and (10301 ), we have 



^gluon 



^gluon ^ ^gluon ^ ^gluon 



(031) 

From the definition (|023|) and the eqs. (|025|) , ( p.26|) and (|031|) , we obtain the 
expression of Z^^ given in eq. ( p.24|) . 

Observe that, in this scheme, the term corresponding to the heavy-fiavour loop 
(eq. (|O30|) ) compensates exactly the self-energy correction to the external gluon line 
coming from the heavy-fiavour loop (eq. (|022|) ): in fact, the final-state gluon is on 
the mass shell, so it is effectively renormalized at zero momentum by the heavy-quark 
loop, and thus decoupling applies. 
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Soft and coUiner amplitudes 



D.l Collinear limit for g — > gg splitting 




K+1 ' 1 

\ 

V 



Figure D.l: Gluon splitting into a gg couple. 

In this appendix we derive the singular part of the square of the invariant amph- 
tude when two colhnear gluons are produced. In the collinear limit, the amplitude 
for the emission of two gluons in the final state can be decomposed into two parts: 
the first one contains the graphs where the two gluons are emitted by a single virtual 
one (see Fig. U.l ), and the other one contains all the other graphs 



A^' = i^A'^Al + k) '^[l^l^y {-9s) r'^T'^'^i-k, -I, k + l)+ e,{k) i^il) , 

(D.l) 

where a and b are the colour indexes of the final gluons, P is the spin projector of the 
gluon propagator, gs is the strong coupling constant, /"'"^ are the structure constants 
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of the SU(3) gauge group, e and e are the polarization vectors of the final gluons, and 
T'^'^^ is the Lorentz part of the three-gluon vertex 

r^'^^'i-k, -I, k + l) = {-k + lygf"" + {-21 - k^g"'' + {2k + lYg"^ . (D.2) 

Only the first term of ( p.lD is singular in the collinear limit. We want to stress the 
fact that this term is singular in the soft limit too. Therefore one has to be careful, 
when considering the soft and collinear limit of the square amplitude, not to include 
this contribution twice. 

We introduce two light-like vectors t and r], that, in the centre-of-mass system, 
have components 

t= (|k + I|, k + I ) 

(D.3) 



T] = C X 



k + 1 



v|k + i|' |k + I|V 
and choose c = 1/4, so that 2t ■ rj = 1. We then decompose, in the collinear limit, 

+ = t^' + ^7]" , (DA) 

where 

^={l + ky = q'y . (D.5) 

We will work in the light-cone gauge, characterized by the light-like vector t], because, 
in this gauge, as we will see, the interference of the divergent term of ( p.l|) and of 
the finite term TZ does not contribute to the singular part. 

In this gauge, characterized by 

V^ = 0, rj-e = 0, r]-i = 0, (D.6) 
the gluon spin projector becomes 

r] ■ p 

and the transversality of gluon polarization gives 

k-e{k)=0, /■€(/) =0. (D.8) 
In addition, for every momentum p 

V.P'^^p) = V^P^'^P) = . (D.9) 
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We write / and k as (Sudakov decomposition) 

= (1 - v) t'^ + Cv^ - k^± , 

with k± such that 

t-k± = , r]-k± = . 

By imposing that k"^ = P = and that (/ + k)"^ = q'^y, we have 

kl = ~v{l - v) q'^y , ^' = {l-v)q^y, ^" = vq^y. 



(D.IO) 



:d.ii) 



fD.12) 



We can exphcitly build the components of k±: in fact, in the centre-of-mass system, 
the space-hke four-vector /c^ can be written 



(0,kj 



(D.13) 



with belonging to the plane spanned by k and 1, and perpendicular to k + 1, so 
that 



ak + 61 



ki • (k + 1) = 



k^-q = Q 



(D.14) 



In this way, the conditions ( p.ll|) are fulfilled, by virtue of eqs. ( p.3|) . It can be easily 
shown that eqs. ( |D.14| ) are satisfied, in the collinear limit, if 



so that 



a (g ■ A;) + 6 (g ■ /) = , 
N , N 



(D.15) 



(D.16) 



q - k q ■ I 

with N a normalization constant, determined by the first formula of eqs. ( p.l2|) , but 
whose value is irrelevant to us, because we are interested in normalized quantities, 

like 

fc'f 1 



^ V^2k~l 
From eqs. (|D.4| ) and (p.lOj) we obtain 

1 



Iq-l 
q ■ k 



Iq ■ k 



k^ 



l-v 
1 



+ {l-2v)q^y7]^ + k^^ 



V 



(1 -v)k^' - (1 -2v)q^y7]^' - k'^ 



(D.17) 



:D.18) 



(D.19) 
(D.20) 
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We are now in a position to simplify the expression of the invariant amphtude (p.l|) 
In fact, the first term of T^'^'^ of eq. ( D.2 ) becomes, with the help of eq. (p.20|) . 



-k + iyP^^{k + l)e^{k)i,{l) = [{l-2v){l + ky -2kl]P^^{k + l)e^{k)i,{l) 

1 - 2 v)2 q^yr]^ - 2 klP^^ik + I)] e^{k) e,{l) , 



where we have used eq. ( p.9| ) and, in the last line, 

{k + iyP^^ik + l) = 2q%r]^ 



(D.21) 



Using eqs. ( p.6| ), ( p.8D , ( p. 181 ) and ( p. 191 ) we can write the second and third term of 
eq. (p.2|) in the following way 



{-21 - k fe^ik) ~e,{l) = -2 l^^e,{k) e,(/) = - k^^e^{k) e,(/) 
{2k + ire^{k) e,(/) = 2 A;%(A;) l,{l) = ^ k'ie^{k) i,{l) . (D.22) 
The amplitude ( p.l| ) becomes 

[ q y 

X -2 kl + ^ A;l ^'^^ + k'i g^^ + O {y)] + nz\e,{k) ~e,{l) . (D.23) 

V 1 — V J J 

Observe that the first term is of order 1 / , so that a singularity with strength 1/y 
can arise only from the square of the first term, and the interference term does not 
contribute. For this reason we can now substitute 



Al{l + k) — > Al{t) = tree-level amplitude 

P'^^k + 1) — > P'^^t) = -g^^ + ^"^^ + = _g'n 

7] ■ t 



(D.24) 



(D.25) 



where we have neglected terms of order k±, which give no contributions to the di- 
vergent part of the amplitude. Keeping only terms that contribute to the coUinear 
singular part, we can write, with the help of eq. (p.6| ). 



q^y 



-2klg'r + -Kgr + T^klg'i" 
V 1 — V 



e^{k)~e,{l). (D.26) 
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By squaring the amplitude and summing over the colours and spins of the final 
gluons, we obtain, for the collinear singular part, 

gl 4CW r _ 1 1 



M 



col 

gg 



r y 

-2v(l 



1 

-2 + - 

V 



V){1 



1-v 



+ v{l-- 



9o 



^cW^r'W, (D-27) 



kl 2 - 2 e 

where we have used the gauge invariance t^A"{t) = to write the following identity 
and 



k^ k±^,y 



d-2 = 2-2e . 



fD.28) 



The first term of eq. ( p.27 ) is recognized to be the Altarelli-Parisi splitting function 
for the gluon-gluon process, in = 4 — 2e dimensions. The second term vanishes after 
azimuthal average in 4 — 2e dimensions, as can be seen from eq. ( [5.28[ ). 

Coming now to our problem, we can further specify the structure of A'^{t)Al°'' (t). 
In fact, by using eq. ( [4.33|) , we can write eq. (p.27|) in the following form 



M 



col 



2 2e 
9 s 2 



Q y I 

-2v{l-v){l-e) 



1 1 

V I — V 

k \ rrk \ rr' 



+ v{l-v) 



2 - 2e 



X M 



aa 
b 



[D.29) 



D.2 Collinear limit for g — > qq splitting 

In this appendix, we derive the singular part of the square of the amplitude for the 
emission of a collinear couple of massless quark-antiquark. The invariant amplitude 
of the process, represented in Fig. p.2| , is 



A, = A'^Ak + I) '-^0^ m {-Wsl%) vil) , (D.30) 

where P is given by eq. (p.7| ) and are the generators of SU(3) gauge symmetry. By 
squaring this amplitude and summing over the spins and colours of the final quarks, 
we obtain 

-^q? = ^^c(^ + 1)^7' (k + l)Pp,{k + l)Pp>Ak + /) Tr {mY') , (D-31) 
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k 



Figure D.2: Gluon splitting into a qq couple. 



srafe 



where are normalized such that Tr 

Considering now eqs. (p.lO|) , we see that, in the colhnear hmit, the trace is of the 
order of y, so that the singular part can be obtained by putting y = in the rest of 
the numerator 



M 



col 

qq 



■A:{t)Ar'{t)g^p,g^,.^,i:T{mY' 



;D.32) 



where we have used the definition of t given in eq. ( p.4| ). Evaluating the trace and 
keeping in the numerator only the terms proportional to j/, we obtain 



M 



col 

qq 



that is 



M 



col 
qq 



A:it)AT it) 



(D.33) 



g2 y 



2 - 2e 



+ 2v(l-v) 



+ {1 -vf 

gi.ua' 



>A:{t)Ar'{t) 



kl 2 - 2e 

Here again we can recognize the Altarelh-Parisi kernel for g ^qq splitting. 



(D.34) 



As previously done for eq. ( |D.27| ) , we can specify this formula to the problem we 
are studying: with the help of eq. (|4.33| ), we can write 



M 



col 

qq 



2 2e 



4Tp 

q^y 



1 



2 - 2e 



k', 



g±aa' 

2 - 2e 



X M 



(7(7 
b 



[D.35) 



D.3. Soft limit for the invariant amplitude QQgg 
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D.3 Soft limit for tke invariant amplitude QQgg 

In this appendix we derive the divergent part of the invariant amphtude for the 
process 

Z/liq) -> Q{p) + W) + 9{k) + g{l) , (D.36) 

in the hmit when the momentum / of the gluon is soft. A soft singularity appears 
only if the soft gluon is emitted from one of the external legs. 

If the emitting external particle is the gluon, the amplitude of the process, in the 
Feynman gauge, is 

Af^^ = A^Jil + k) i-g,) r'^T^^'^i-k, -I, k + I) e,{k) ^1) , (D.37) 

where we have added to eq. ( p.l|) the colour indexes i,i of the produced quarks. 
As / goes to zero, this term develops a singularity: in fact, by using the gauge condition 
k" A^^^ik) = and eq. ( p.8|) , we can write this amplitude as 

^g'fe) = gj'^'^'^ A- A^J{k) ^Al) + non-singular terms. (D.38) 

Similarly, if we consider the emission of a soft gluon of colour index h from an 
external quark leg with momentum p and colour index i, that is 

Qn{v + l)^Qi{v) + gb{l) , 

we can write the invariant amplitude 

^y^^ = u{p){-igsYtl) ^ZiP + e,ik) ~e.il) , (D.39) 

where A refers to the rest of the process from which the quark external line takes 
origin. In the limit of I going to zero, we can rewrite this amplitude as 



^afe(Q) ^ ^^/^^p) eJk) e^(/) + non-singular terms, (D.40) 

p ■ I ■' 

where we have defined 



AZip) = uipHTAp) ■ (D-41) 

In the same way, we can obtain the limit of the amplitude for the soft emission 
from an antiquark with momentum p' and colour index j 

Af^^ = -g, 4^ AZ{p) 4 e.(0 + non-singular terms. (D.42) 
p ■ L 
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Considering that A^J of eq. ( D.3S ) can be written as 



(D.43) 



where A'^ does not contain any colour element, and the similar ones for eqs. ( |D.40| ) 
and (p.42|) , we can sum the three amplitudes to obtain 



At' = 9s < if 



k-l 



where we have disregarded the non-singular terms. 

By squaring the amplitude and summing over the spins and colours of the final 
gluons and quarks, we have 

p ■ k ^ p' ■ k 



(p-O ik-i) + (p'-o (k-i) 



Ca 
2 



p ■ p 



m 



+ 



m 



{p-lf (p'-l) 



X M^^ (D.45) 



{p ■ I) ip' ■ I) 

where we have made use of eq. ( [4.33[ ). 

The same result applies in the case of k soft, once the interchange / <-> /c is made. 



Appendix E 



List of integrals for tke four-jet singular 
soft contributions 



In this appendix, we summarize the values of the integrals required to isolate the 
singular terms of the four-jet cross section, in the soft-gluon hmit. The form of the 
integrals we are going to compute is the following 

rdcj>ism4>)-''^ Tdy ['dvf{y,v,<P). (E.l) 
iV<^ Jo Jo Jo 

We fix our attention to the y and v variables and we split the integration range into 

two triangles 

rx rl rx fv/x rx rl 

dy dvf{y,v,(j))= dy dvf{y,v,(p)+ dy dv f{y,v,(j)) . (E.2) 

Jo Jo Jo Jo Jo Jy/x 

We perform the following change of variables: v — yt in the first integral and v — y/t 
in the last one. In this way, we obtain 

'l/x 



px pi px t-i/x px y f y 

I dy dvf{y,v,4>) = dy dty f{y,yt,4>) + dy dt — f[y,-, 

Jo Jo Jo Jo Jy V \ t 

pi/x px r y f y \ 

= dt dyy f{y,yt,(^) + dt dy — f [y, -, (f)\ , 

where, in the second line, we have exchanged the order of integration. 
In this way we succeed in writing I as the sum of two contributions 

I = 7(i) + /(2) (E.3) 

/W = ^ r d(t>{sin(t>)-^^ j"^ dt r dyyf{y,ytA) (E.4) 
iV^ Jo Jo Jo 
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Ii-type integrals 



The first type of integrals to compute is 

h{x,h) = r d<p{sm<f))-'^ r dy C dv[v{l-v)Y'y-' \ 

IS^ Jo Jo Jo y[y + n v\ 



(E.6) 



As f{y,v,(j)) = l/y[y + hv] is 0- independent, the integration is straightforward 
and, according to eq. ( |E.4| ), we have 



= / ^dt Tdyy-'-'^il-ytr- 
Jo Jo 1 

dyy-'-'^ 



l/x f-e 

dt- , 

1 + ht Jo 



+ ht 

I- e\og{l-yt) + (e^ 



l/x f-e 

dt — 

1 + ht 



It 



-X 



-2e 



X 



y 



(E.7) 



The integral in y is finite so it gives contribution of order e that can be neglected. In 
this way we have 



I? = -j-x-'^l^'^^dt{l-e\ogt + 0{e'))^ 



ht 



2hi 



X 



-2e 



logU + ^j +e 



log X log 1 + - - Ll2 



(E.8) 



where terms of order e are not computed. The integral of eq. ( [E.5|) gives rise to 



(2) 



X rt 



JO 



dt dyy 



-l-2e 



1 - - 



ft 

t + hJo 



r(2) , r(2) 
-'la "T -'ife ! 



tj t+h 



l-elog(l-| 



O e 



(E.9) 



where, for simplicity, we have separated out the integration into two terms. Perform- 

(2) 

ing the y integration in /^^ , we get 

1 



t(2) 
-'la 



1 , t"^-' 
— / dt- 



2e Jo t + h 
2he \ e 



dt { - 



1 



2 he Jo \t t + h 
1 -elogt + 0(e2)l 



dt 



t + h 



2hi 



X 



log ( 1 + - ) + e 



log X log ( (1 + + Li2 (^-^ 



(E.IO) 
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where we have neglected terms of order e. To compute /^^ , we perform the change 
of variable z = y/t to obtain 



7-(2) 
'l6 



■ , t-^-' /•! , log(l - Z) 

-e I dt / dz ' 



t + h 



yl+2e 



TT^l 

^ 6 h 



. (E.ll) 



Collecting all terms together, according to eqs. (|E.3|) , ( |E.8| )- (|ErTl| ), we have 

,2 



/,(.,.)^±U-llog.-log^^ + ilog^.-!^-2L.(-^ 



0{e) . 
(E.12) 



l2-type integrals 

The second type of integrals has the following expression 

l2{x,h) = r d<P{sm<f))~'^ rdy C dv[v{l-v)Y'y-' ^ . (E.13) 
i\tf) Jo Jo Jo [y + av\'^ 

Here again, the integration is straightforward, and we have, according to eq. 

i-l/x PX 



, , fl/X fX 

= / dt dyy-'-'^il-yt) 
Jo Jo 



dt- 



(l + ht) Jo 

l/x 



{i + hty 

dyy-'-^^ \l-e\og{l-yt) + 0(e^ 



1 /"i/^ 1 

— / rft(l-elogt) o+0(e 

2e Jo ^ ^ ^ (l + ht? ^ 



<! 1 - 



X 



2 he [ X + h 
From eq. ([E.5|) we obtain 



+ e 



{i + hty 

\ (. h\ h ' 

log 1 + - H -^loga; 

V X X + h 



+ 0{e) .(E.14) 



/f - fdtfdyy-^-^^U-y'-' 
Jo Jo 



t) {t + hy 



x rt 

dt o / dyy 

{t + hfJo 

1 , t-' 
-— dt- 



-l-2e 



y 



l-elog(l-|)+0(e2 



2eJo {t + hy 
1 { X 



-7^ Tc/t (l -elogt + Ofe 
It Jo ^ ^ 



2he \ x + h 



h X 
log — — + — — logx 

X + h X + h 



(t + hy 

+ 0{e) . (E.15) 
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Adding together these two contributions, we have 

h{x, h) = J« + = -L |_i _ 2 log + ^) + \ogh\ + 0{e) . (E.16) 



Is-type integrals 

The last type of integrals has the following form 

/3 = -^ r#(sin0)-2V3(x) , (E.17) 
N4, Jo 

where 

hix) ^ r dy f dv [v{l - v)]-^y-^^ . V ' (^'l^) 

JO Jo y + hv y — ccoscp^yjv + gv 

According to eq. ( [E.4|) and ( [E.5|) , we have 

= 1^''' dt r dyy~^~^'{l-yt)-' — 

•^0 ^0 (1 + /it) (1 - ccos^Vt + 



— e 



-±. a;-2e / ^ + O (e) 

2e Jo {1 + ht) [1 - c cos (t)Vt + gtj 

(1 - 2elogx) / dt ^ ^ ^ + O (e) 

2e' 'Jo {l + ht)(l-ccos(PVi + gt) 

1 /-V^ , 1 
/ dt- 



2eio + (l-ccos0v^ + ^t 
+ f'" it V + O (.) (E.19) 

■''0 + (l-CCOS0v^ + ^t) 

Jo Jo 



ty {t + h) it - c cos 4)^/1+ g 



-T^ / 7 F ^+0{t) 



I dt 7^ ''-^^ r + 0(e] 



-— I dt 



it + h)(t- 


- c cos (f)y/i + g^ 


1 


— e log t 


it + h)(t' 


- ccos0\/t + g^ 




1 



2e^o {t + h) (t- ccos(pVi+ g) 
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rx 

+ / dt- 



hogt 



{t + h) (t — c cos + g 
From eqs. ( |E.17|) and ( [E.18| ) we have 



+ 0(e) . 



h = 



1 

N^Jo 
1 m 



r(2), 



[1 — 2e log (sin i 



'l!,'\x) + lP{x) 



N^Jo 

U-^/. + ^[/3. + /3.]+/. + 0(e)} 



N, 



with 



^3a 



roo rn 

dt ^ 
Jo Jo 



rdt / 

JO Jo 



dt 

"'O 



JO JO 



(l + ht) (l - 


- c cos (j)^/t -+ 


-gt) 


2 log 


X + log t 




(l + /it) (l - 


- CCOS0v^ + 


gt) 


lo 






{h + t) (t — c cos (pVi + g 


) 


log 


; (sin 0) 




(l + ht) (l - 


- C COS (f)^/t -+ 


-gt) 



We can now integrate in the first three expressions, using the identity 

/ d(f) ; = : , 

JO a — b cos (J) V — 6^ 

and in t the last one, to obtain 

■ 1 71 

dt 



(l + ht) ^(l + ^i)2_^2^ 

h (c2 -2g + 2h + 2^{g - hy + c^h 



ha 



'36 



^^^=^^^= log 

(^ - + c^h h^-{g- ^{g - hf + c^h 

k , 2 log a; + log t n 
dt 



(1 + ^^) J(l + ^tf-c2t 



r „ logt 

dt 



(/i + t) ^(t + gf _ 



(E.20) 



(E.21) 

(E.22) 
(E.23) 
(E.24) 
(E.25) 



(E.26) 

(E.27) 
(E.28) 
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log (sin (p) 2 

{g — hy + h (ccos0) 



f2c{g + h)cos(j) CCOS0 g^ 

[ y/Ag — & cos^ v4(7 — & cos^ /i ^ 



i^^iogv/r^ 



v^r^ {g-h)^ + hcH^ 

{2 c(q + h) t ct , , X , 9 1 N 

where we have used the symmetry of the integral around the point it/ 2 and made 
the change of variable t = cos 0. 

We are now in a position of further reduce the expression of I^, to account for 
the actual values of the parameters: in fact, with the help of eqs. ( |5.15| )-( |5.17| ), ( [B.(j| ) 
and ( |B.49| ), we have 
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